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Abstract 



This thesis is dedicated to the holographic study of flavor dynamics. The technique 
employed is a D7-brane probing of various D3-brane backgrounds. The first topic 
covered studies the influence of an external magnetic field on a fiavored large N Yang- 
Mills theory. The theory exhibits spontaneous chiral symmetry breaking. A discrete 
self-similar structure of the spontaneous symmetry breaking mechanism is observed, 
this structure is reflected in the meson spectrum. The meson spectrum exhibits Zee- 
man splitting and characteristic GMOR relation. The second topic examines thermal 
properties of the dual gauge theory. The study reveals a flrst order phase transition 
associated to the melting of mesons. The critical ratio of the bare quark mass and 
temperature at which the transition happens is computed. The third topic studies the 
phase structure of the flnite temperature dual gauge theory in the presence of magnetic 
fleld. A phase diagram of the theory is obtained. The temperature restores the chiral 
symmetry while the magnetic fleld has a freezing effect on the meson melting. The 
meson spectrum exhibits Zeeman splitting and characteristic GMOR relation. Ther- 
modynamic quantities such as free energy, entropy, and magnetization are computed. 
The fourth topic studies the addition of an external electric fleld. The observed effect 
is dissociation of the bound quarks, favoring the meson melting. For sufficiently strong 
electric fields a global electric current is induced. Thus, the dissociation of mesons 
corresponds to an insulator /conductor phase transition. This transition persists at 
vanishing temperature. The fifth topic studies the addition of an R-charge chemical 

XV 



potential via brane probing of the spinning D3-brane geometry. The corresponding 
phase diagram is obtained. The chemical potential favors the dissociation of mesons. 
For high chemical potential, a finite phase difference between the bare quark mass and 
the quark condensate is induced. The last topic explores universal properties of gauge 
theories dual to the Dp/Dq system. A universal discrete self-similar behavior asso- 
ciated to the insulator /conductor phase transition is observed and the corresponding 
scaling exponents are computed. A similarity between the electric field and R-charge 
chemical potential cases is discussed. 
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Chapter 1: Introduction 



1.1 Introduction to the AdS/CFT correspondence 

The goal of theoretical physics is to provide an economical and self-consistent descrip- 
tion of physical reality, by means of physical laws and first principles. This implies 
borrowing the structure of mathematics, while using the concepts of philosophy. How- 
ever, the interplay between these subjects has always been mutual, as the discovery of 
new physical phenomena required the development of both novel philosophical concepts 
and mathematical tools. Indeed, it is now diflacult to draw a solid line between math- 
ematical and theoretical physics, while the philosophical aspects of general relativity 
and quantum mechanics remain challenging. 

Despite the broad meaning of the term theoretical physics, a considerable number 
of theoretical physicists are interested in the study of the fundamental interactions of 
matter. It is well established that there are four basic interactions of matter, namely, 
the electromagnetic, weak nuclear, strong nuclear, and the gravitational interactions. 
It is somewhat ironic that albeit the concept of gravitational interaction was the first 
to be developed, it remains a challenge to come up with a consistent quantum descrip- 
tion of gravity. It is believed that string theory, which has the graviton present in 
the spectrum of the fundamental string, is the best candidate of quantum theory of 
gravity. Furthermore, string theory provides a natural framework for a unified theory 
of fundamental interactions [39, 49, 75]. 
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Historically, string theory emerged as an attempt to delineate the strong interac- 
tions by what was called dual resonance models. However, shortly after its discovery 
Quantum Chromo Dynamics (QCD) which is a SU{3) Yang-Mills gauge theory, super- 
seded it. The matter degrees of QCD consist of quarks which are in the fundamental 
representation of the gauge group, while the interaction between the fundamental fields 
is being mediated by the gluons which are the gauge fields of the theory thus trans- 
forming in the adjoint representation of SU{3). A remarkable property of QCD is the 
fact that it is asymptotically free, meaning that at large energy scales, or equivalently 
at short distances, it has a vanishing coupling constant. This makes QCD perturba- 
tively accessible at ultraviolet. However, the low energy regime of the theory is quite 
different. At low energy QCD is strongly coupled, the interaction force between the 
quarks grows immensely and they are bound together, they form hadrons. This phe- 
nomenon is called confinement. Additional property of the low energy dynamics of 
QCD is the formation of a quark condensate which mixes the left and right degrees 
of the fundamental matter and leads to a breaking of their chiral symmetry. It is 
extremely hard to examine the properties of the strongly coupled low energy regime of 
QCD, since the usual perturbative techniques are not applicable. We need to coin new 
non-perturbative tools. One such approach is lattice QCD which has had a tremen- 
dous success in describing the static properties of the theory, with the use of large scale 
computer simulations. However, some of the most interesting and important problems 
of QCD, such as the description of the low energy mechanisms for confinement and 
spontaneous chiral symmetry breaking remain unsolved. This is why it is of extreme 
domination to develop novel analytic techniques applicable in the strongly coupled 
regime of Yang-Mills theories. 

The AdS/CFT correspondence, as we shall demonstrate in this section, is a pow- 
erful analytic tool providing a non-perturbative dual description of non-abelian gauge 
theories, in terms of string theory defined on a proper gravitational background. Let 
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us proceed by tracing the sequence of the ideas that lead to the development of this 
gauge/string correspondence. 

A crucial milestone was the large N limit proposed by t'Hooft [84]. Instead of using 
SU{3) as a gauge group, t'Hooft proposed to consider SU (N) Yang-Mills theory and 
take the limit N ^ oo, while keeping the so-called t'Hooft coupling fixed A = Ngyj^- 
t'Hooft proved that in this limit only planar diagrams contribute to the partition 
function which makes the theory more tractable. On the other side, the expansion in 
corrections of the QCD partition function and the genus expansion of the string 
partition function exhibit the same qualitative behavior, suggesting that perhaps a dual 
description of the large N limit of non-abelian gauge theories might be attainable in 
the frame work of string theory. 

These early hints about a possible gauge/string duality came very close to real- 
ity with the development of the concept of Dp-branes and their identification as the 
sources of the well-known black p-brane solutions of type HB supergravity. The key 
observation was that the low energy dynamics of a stack of coincident Dp-branes 
can be equally well described by a SU{N) supersymmetric Yang-Mills theory in p + 1 
dimensions and an appropriate limit of a p-brane gravitational background. The first 
gauge theory studied in this context is the large N J\f = 4, SU{N) supersymmetric 
Yang-Mills theory in 1 + 3 dimensions which is a maximally supersymmetric con- 
formal field theory. The corresponding gravitational background is, as proposed by 
Maldacena [66], the near horizon limit of the extremal 3-brane solution of type HB 
supergravity. This was the original formulation of the standard (by now) AdS/CFT 
correspondence. 

On the other side, by the time when the AdS/CFT correspondence was proposed, 
there was enough evidence that a consistent theory of quantum gravity should ex- 
hibit holographic properties [83], more precisely the dynamics of the physical degrees 
of freedom contained in a certain volume of space-time should be captured by the 



3 



physics of its boundary. As we shall demonstrate, the framework of the AdS/CFT 
correspondence provides a natural description of string theory on certain gravitational 
background, in terms of a dual field theory defined at a proper asymptotic boundary. 
This is how the AdS/CFT correspondence satisfies the holographic principle. 

1.1.1 Low energy dynamics of D3— branes 

Let us consider a stack of N coincident D3-branes. This system has two different kinds 
of perturbative type IIB string theory excitations, namely open strings that begin and 
end on the stack of branes and closed strings which are the excitations of empty space. 
Let us focus on the low energy massless sector of the theory. 

The first type of excitations corresponds to zero length strings that begin and end 
on the D3-branes. The orientation of these strings is determined by the D3-brane 
that they start from and the D3-brane that they end on. Thus, the states describing 
the spectrum of such strings are labeled by Xij, where i,j = l,...,N. It can be shown 
that in the case of oriented strings [49] Xij transform in the adjoint representation of 
U{N). On the other side, the massless spectrum of the theory should form a A/^ = 4 
supermultiplet in 1 + 3 dimensions. The possible form of the interacting theory (if 
we take into account only interactions among the open strings) is thus completely 
fixed by the large amount of supcrsymmctry that we have and is the J\f = 4 U (N) 
supersymmetric Yang-Mills theory. Note that U{N) came from the transformation 
properties of Xij. On the other side, U{N) can be thought of as a direct product of 
L'"(l) and SU{N), geometrically the U{1) corresponds to the collective coordinates of 
the stack of D3-branes. We will restrict ourselves to the case when those modes were 
not excited, we refer the reader to ref. [1] for further discussion on this point. 

The second kind of excitations is that of type IIB closed strings in flat space. The 
low energy massless sector is thus a type IIB supergravity in 1 + 9 dimensions. 
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The complete action of the system is a sum of the actions of those two different 
sectors plus an additional interaction term. This term can be arrived at by covari- 
antizing the brane action after introducing the background metric for the brane [64]. 
It can be shown [1] that in the a' — > limit the interaction term vanishes and the two 
sectors of the theory decouple. 

To summarize: the low energy massless perturbative excitations of the stack of 
D3-branes are given by two decoupled sectors, namely A/" = 4 SU (N) supersymmetric 
Yang-Mills theory and supergravity in flat 1+9 space-time. Our next step is to consider 
an equivalent description of this system in terms of effective supergravity solution. 

1.1.2 The decoupling limit 

Let us consider the extremal black 3-brane solution of type IIB supergravity. The 
appropriate gravitational background is given by [49] : 

ds'^ = H~^^\^dx^'dx'' + Hl^'^dx'dx' , (1) 

2$ 2 

e = gs , 

C(4) = H^^g'^dx^ A • • • A dx^ , 

where /x = 0, . . . , 3, i = 4, . . . , 9, and the harmonic (in six dimensions) function is 
given by: 

^3 = 1 + 1^25^. (2) 

The integer number N quantizes the flux of the five- form field strength dC^^y It can 
also be interpreted as the number of D3-branes sourcing the geometry. Taking the 
near horizon limit of the geometry corresponds to sending a' — > while keeping the 
quantity u = r/a' fixed. Such a limit serves two goals: first it enables one to zoom 
in the geometry near the extremal horizon and second it corresponds to a low energy 



5 



limit in the string theory defined on this background. After leaving only the leading 
terms in a', one can obtain the following metric [66]: 

ds^ = ^{-dxl + dxl + dxl + dxj) + + R^dnl , (3) 

1 

C(4) = — -^dx^ A dx^ A dx'^ A dx^ , 
9 s 

e = 9s , 

The background in equation (3) is that of an AdSs x space-time of radius R. Note 
that from a point of view of an observer at r — > oo, the type IIB string theory excita- 
tions living in the near horizon area, namely superstring theory on the background (3), 
will be redshifted by an infinite factor of = H^^^*. Therefore, we conclude that 
type IIB superstring theory on the background of AdSs x should contribute to the 
low energy massless spectrum of the theory seen by an observer at infinity. However, 
an observer at infinity has another type of low energy massless excitations of type IIB 
string theory, namely type IIB supergravity on flat 1 + 9 space or gravitational waves. 
Those two different types of excitations can be shown to decouple form each other. 
To verify this one can consider the scattering amplitudes of incident gravitons of the 
core of the geometry (the near horizon area). It can be shown that at low energies 
{lo <^ ^/R) the absorption cross-section of such a scattering (Tabs goes like [42, 56] 
(Tabs c)c LO^R^. Therefore, one has that (Tabs and those two types of excitations 
decouple in the low energy limit u ^ 0. 

Similar to the description from the previous subsection, we learned that the mass- 
less low energy spectrum of the theory contains two types of perturbative excitations, 
namely type IIB string theory on the background of AdSs x S^, and type IIB su- 
pergravity in flat 1 + 9 space-time. Now we are ready to conjecture the AdS/CFT 
correspondence. 
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1.1.3 The AdS/CFT correspondence 

As we learned from the previous subsections, the massless sector of the low energy 
dynamics of N coincident D3-branes allows two possible descriptions. Conjecturing 
that these descriptions are equivalent is the core of Maldacena's AdS/CFT correspon- 
dence [66] . Notice that in both descriptions one part of the decoupled sectors is a type 
IIB supergravity in flat space. Thus, we are naturally led to the conclusion that type 
IIB superstring theory on the background of AdS^ x background is dual to Af = 4 
SU{N) supersymmetric Yang-Mills theory in 1 + 3 dimensions. We have presented 
this statement in a diagrammatic way in Figure 1. 



type IIB superstring theory 
on the background of AdSs x 



+ 



type IIB supergravity 
in fiat 1 + 9 space-time 



AT = 4 SU{Nc) SUSY YM + 



type IIB supergravity 
in fiat 1-1-9 space-time 



Figure 1: A diagrammatic statement of the AdS/CFT correspondence. 

A further hint supporting the AdS/CFT correspondence is that the global sym- 
metries of the proposed dual theories match. Indeed, an AdSs space-time of radius R 
can be embedded in a M^'^ flat space. It can then be naturally described as a hyper- 
boloid of radius R and thus has a group of isometry SO{2,Q) which is also the group 
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of rotations in 2 + 4 dimensions. On the other side, the part of the geometry has a 
group of isometry SO (6) (the group of rotations in 6 dimensions). This leads us to the 
conclusion that the total global symmetry of string theory on AdS5 x gravitational 
background is 50(2,4) x 5*0(6). It is satisfying that the corresponding gauge theory 
has the same global symmetry. Indeed, it is a well-known fact that the = 4 super- 
symmetric Yang-Mills theory in 1 + 3 dimensions is a conformal field theory. As such 
it should has the global symmetry of the conformal group in 1 + 3 dimensions, but this 
is precisely 50(2, 4). Actually, since the theory is supersymmetric, the full global sym- 
metry group is the superconformal group which includes an SU (4) global R-symmetry. 
In particular this group rotates the gauginos of the super Yang-Mills theory. However, 
it is well known that 5C/(4) = 50(6), and therefore the global symmetry of the gauge 
theory is indeed 50(2,4) x 50(6)! 

An important aspect of the correspondence is the regime of the validity of the 
dual description. Depending on the precise way in which we are taking the a' ^ 
limit, there are two basic forms of the AdS/CFT correspondence. The strongest form 
is that the string/gauge correspondence holds for any A''. Unfortunately, this strong 
form of the conjecture cannot be tested directly since it is not known how to quantize 
superstring theory on a curved background in the presence of Ramond-Ramond fluxes 
[26]. The second form of the conjecture holds in the t'Hooft limit when N ^ oo 
and the t'Hooft coupling A oc ggN is kept fixed. In this way on gauge side of the 
correspondence only planar diagrams contribute to the partition function while on 
string side the required gs ^ limit suggests semiclassical limit of the superstring 
theory on AdSs x 5^. An important observation is that large A 3> 1 suggests large 
AdS radius R oc A^/^ and hence small curvature of the AdS background. This implies 
that the supergravity description is perturbative and thus provides an analytic tool for 
perturbative studies of non-perturbative field theory phenomena. On the other side, if 
we are at weak t'Hooft coupling (A <S 1) we can perform perturbative studies on gauge 



8 



side of the correspondence and transfer the result to the pert ur bat ivcly inaccessible 
regime of the supergravity description. Therefore, we conclude that the AdS/CFT 
correspondence is a strong/weak duality. In this work we will concentrate solely on 
the study of strongly coupled A ^ 1 Yang-Mills theories. Hence, we will perform the 
analysis on the supergravity side of the AdS/CFT correspondence. 

1.1.4 The AdS/CFT dictionary 

Let us focus on the precise way that the AdS/CFT correspondence is implemented. 
After a closer look at the geometry of the AdSg x background, we conclude that 
it has five non-compact directions. Four of them are parallel to the D3-branes world 
volume and correspond to the 1 + 3 directions of the dual gauge theory. The fifth 
non-compact direction is the radial direction u (radial in the transverse, to the D3- 
branes, space) and its interpretation in the dual gauge theory is not obvious. To 
shed more light on it, let us consider the action of a free massless scalar field in 1 + 3 
dimensions [26]: 



The corresponding field theory is conformal and thus has a global symmetry 5'0(2,4) 
which is the conformal group in 1 + 3 dimensions. Therefore, we can consider the trans- 
formation properties of the scalar field under the action of the dilatation operator. 
One can verify that the transformation: 



leaves the action (4) invariant. Furthermore, we learn that the scalar field (f) has a 
scaling dimension one. On the other side, the 50(2,4) group is the group of isometry 




(4) 



X e"x; — s- e "(j); 



(5) 
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of AdS5 and one can verify from equation (3) that the transformation x — > e^x, 
suggests: 

u e~'^u . (6) 

Therefore, we learn that the coordinate u scales as energy under dilations and thus 
has a natural interpretation as an energy scale of the dual gauge theory. 

The further development of the AdS/CFT correspondence resulted in a map be- 
tween gauge invariant operators in A/" = 4 super Yang -Mills in a particular irreducible 
representation of the SU (4) R-symmetry group and supergravity fields in the isomor- 
phic representation of the 50(6) global symmetry. These representations are obtained 
after Kaluza-Klein reduction of the supergravity fields on the internal sphere. Let 
us consider for simplicity the case of a scalar field of mass m, propagating on the 
AdSd+i background. The relevant action is [26]: 

S = J d'^xduy/^{g''''da(t>di,<l) - m^^) • (7) 

The solution of the corresponding equation of motion have the following asymptotic 
behavior at large u: 

</)(n,x)= (^i) <t^o{x)+{^ (0(x)), (8) 



where 



A = f + y^ + i?W. (9) 

Note that the supergravity field ^{u) is a scalar field and is thus invariant under the 
action of the dilatation operator because the latter is one of the generators of the 

global symmetry S'0(2,4). Therefore, we conclude that and {0{x)) carry scaling 
dimensions 4 — A and A, respectively. In [41] it was suggested that and {0{x)) 
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correspond to the source and the vacuum expectation value of the gauge invariant 
operator 0{x). It is also worth noting that the expression: 

J d'^x^Q{x){0{x)) = inv. (10) 

is invariant under the 5*0(2,4) global symmetry. It was suggested that the exact form 
of the map is given by the relation [41, 86]: 

^g/d<^.0o(x)(O(x)>^^p^ = 2sugra[0o(x)] , (11) 

where 

^Sugra[0o(a;)] = lim2:sugra[95'o(l/e,a;) = (^'^~^(t)Q{x)\ . (12) 

i.e. the generating functional of the conformal field theory coincide with the generating 
functional for tree level diagrams in supergravity. We refer the reader to the extensive 
review ref. [1] for more subtleties on the precise way of taking the e — limit in 
equation (12). 

Formula (11) has been tested by comparing correlation functions of the = 4 
quantum field theory with classical correlation functions in AdS^. Note that the tree 
level approximation on supergravity side is valid only at strong t'Hooft coupling and 
therefore the corresponding conformal field theory is strongly coupled. This is why 
the correspondence was tested in this way only for correlation functions which satisfy 
non-renormalization theorems and hence are independent on the coupling [26]. In 
particular it applies for the two- and three- point functions of 1/2 BPS operators 
[31, 63]. 

Further checks of the correspondence beyond the 1/2 BPS sector was started with 
the so-called plane-wave string/gauge theory duality, where one takes appropriate 
plane-wave limit of the AdS5 x background [15]. Key point of this limit is that 
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superstring theory on this background can be exactly quantized. Recently a significant 
progress towards quantizing superstring theory on AdS^ x has been achieved using 
integrable spin chain models. We refer the reader to refs. [14, 74] for more details on 
these vast subjects. 

1.2 Introducing flavor to the correspondence 

Direct consequence of the confining property of QCD is the fact that the low energy 
dymanics of the theory is governed by color singlets, such as mesons, baryons and 
glueballs. Mesons and baryons are bound states of quarks, the latter transform in the 
fundamental representation of SU{3). The fact that at low energy QCD is strongly 
coupled suggests that it is not accessible for perturbative studies. This is why it is 
important to come up with an alternative non-perturbative techniques describing the 
strongly coupled regime of Yang-Mills theories and in particular Yang-Mills theories 
containing matter in the fundamental representation of the gauge group, such as QCD. 

Further need of alternative non-perturbative techniques applicable to the prop- 
erties of the fundamental fields in the strongly coupled regime of non-abelian gauge 
theories is required by the very recent discoveries of the properties of matter obtained in 
heavy ion collision experiments. More precisely the fact that the quark-gluon plasma 
which is the phase of matter of the fireballs obtained in such experiments, is not the 
expected weakly coupled quark-gluon plasma predicted by the standard perturbative 
QCD but is classified as a strongly coupled quark-gluon plasma. A novel phase of 
matter that provides challenge for the society of theoretical physicists. 

One of the purpose of the study of the AdS/CFT correspondence is to develop the 
above mentioned analytic tools for the study of strongly coupled Yang-Mills theories. 
The original form of the conjecture that we described in the previous section, focuses 
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on a gauge theory with a huge amount of symmetry, namely the = 4 super Yang- 
Mills theory. On way to make the correspondence more applicable to realistic gauge 
theories, such as QCD, is to reduce the amount of the supersymmetry of the theory 
by introducing additional gauge invariant operators. This approach though fruitful 
still has the weakness that the matter content of the dual gauge theory, more precisely 
the fermionic degrees of freedom, transform in the adjoint representation of the gauge 
group. In other words there are no fundamental fields in the theory. The reason is 
that both ends of the strings, producing the field content of the gauge theory, are 
attached to the same stack of D3-branes and the corresponding states transform in 
the adjoined representation of the gauge group. In order to introduce fundamental 
matter, one needs to consider separate stack of D-branes. 

The easiest way to introduce fundamental fields in the context of the AdS/CFT 
correspondence is to consider an additional stack of Nf D7-branes [50]. (Prom now on 
we will use AT^ as a notation for the number of the D3-branes sourcing the gravitational 
background.) Since the D7 branes' world volume is higher dimensional and non- 
compact in the transverse to the D3-branes dimensions, the D7-branes have infinite 
"internal" volume and thus their gauge coupling vanish making their gauge symmetry 
group a global symmetry. In this way we introduce family of fundamental matter 
with global flavor symmetry SU {Nf). To be more precise let us consider two stacks of 
parallel Nc D3-branes and Nf D7-branes embedded in the following way: 








1 


2 


3 
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5 
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D3 






















D7 























Table 1: Embedding of the flavor D7-branes. 

The low energy spectrum of the strings stretched between the D3- and D7-branes 
directions give rise to the M = 2 hypermultiplet containing two Weyl fermions of 
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opposite chirality coming from the light-cone modes of strings stretched along the NN 
and DD directions (2,3,8,9) and two complex scalars coming form strings stretched 
along the ND directions, namely 4,5,6,7. Now if we consider Nf <^ Nc and take the 
large Nc limit. We can substitute the stack of D3-branes with a AdS5 x space and 
study the Nf D-branes in the probe limit using their Dirac-Born-Infeld action. On 
gauge side this corresponds to working in the quenched approximation (Nf <C Nc) and 
taking the large Nc t'Hooft limit. If the D3- and D7-branes are separated in their 
transverse (8,9)-plane, then the strings stretched between them has a final length and 
hence final energy. It can be shown that [49] the mass of the hypermultiplet is given 
by the energy of the string or the distance of separation L multiplied by the string 
tension nig = L/(27ra'). 

Let us study closer the symmetry of the set up. If the D3- and the D7-branes 
overlap the hypermultiplet is massless (m^ = 0). In this case the SO (6) rotational 
symmetry of the transverse space is broken to the product 50(4) x SO{2), corre- 
sponding to rotations along the ND directions (4,5,6,7) and the DD directions (8,9), 
respectively. This is equivalent to a SU{2)l x SU{2)ii x U{1)r global symmetry, and 
suggests that the gauge theory has a R-symmetry group SO{2)ji x U{l)jt [59], which 
is indeed the case, when the hypermultiplet is massless. If the D3- and D7-branes are 
separated it is known that the R-symmetry is just SU{2)fi, which again fits that fact 
that the SO{2) rotational symmetry in the (8,9)-plane is broken. 
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1.2.1 The dictionary of the probe brane 

Let us now focus on the precise way that the AdS/CFT dictionary is implemented. 

The dynamics of the D7-brane probe is described by the Dirac-Born-Infeld action 
including the Chern-Simons term [49]: 



^ = -/^7 J e-^d^C^-det{P[Gab] + + ^^^A^r j P[C(4)] A A , 



iVV 

Ms, Ms, 

(13) 



where {2'Ka')Tah = -P[-Ba6] + (27rQ;')Fa6 and /X7 = [(27r)^a'^] ^. It is convenient to 
introduce the following coordinates: 

p = ucos6; L = usinO; (14) 

and consider the ansatz: 

L = L{p); (j) = const . (15) 
Then the lagrangian describing the D7-brane embedding is: 

C oc pV1 + ^'(/>)' (16) 
leading to the equation of motion: 



At large p the solution has the behavior: 



c 



L{p) = m+-^ + ... . (18) 
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Now if we introduce the field: 



XN= 2^t] \2 = -m+\c+...; u^ = p^ + L{pf, (19) 

we can see that x(u) has the same behavior as the field <p{u,x) from equation (8). 
This is quite suggesting. The asymptotic value of L(oo) = m is exactly the separation 
of the D3- and D7-branes and is thus related to the mass of the hypermultiplet via 
niq = m/(27ra'). Since the hypermultiplet chiral fields are our quarks we will call rriq 
the bare quark mass. Therefore, the coefficient c should be proportional to the vev 
of the operator that couples to the bare quark mass but this is the quark condensate! 
This is an example of the how the generalized AdS/CFT dictionary works at the level 
of a D7-brane probing. Let us provide the exact relation between the quark condensate 
(Oq) and the coefficient c: 

We refer the reader to the appendix of Chapter 2 for more details on the last calculation 
and to ref. [53] for an elegant presentation of the holographic renormalization of probe 
D-branes in AdS/CFT. 

Now let us go back to equation (17) and note that in order for the D7-brane to 
close smoothly in the bulk of the geometry, we need to impose L'{0) = 0. This is 
possible only for c = and thus we conclude that the condensate of the theory vanish. 
But the dual gauge theory is supersymmetric this is why it is not surprising that the 
quark condensate is zero. Furthermore, since there is no potential between the D3-^ 
and D7-branes (because of the unbroken supersymmetry) , the D7-brane should not 
bend at infinity and this is why the solution for the D7-brane embedding should be 
simply L = m, as it is. 

Note that the analysis that led to equation (20) requires that the gravitational 
background be only asymptotically AdSs x S^. In fact, in all cases that we are going 
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to consider in this work, there will be some sort of the deformation of the bulk physics, 
coming either from the gravitational background or from the introduction of external 
fields. This will break the supersymmetry and will capacitate the dual gauge theory 
to develop a quark condensate. In Chapter 2, we will use this approach to provide a 
holographic description of magnetic catalysis of chiral symmetry breaking. Through 
the rest of the thesis, the study of the quark condensate as a function of the bare 
quark mass will enable us to explore the phase structure of the dual gauge theory and 
uncover a first order phase transition associated to the melting of the light mesons 
of the theory. Let us conclude by redirecting the reader to the very recent extensive 
review on the subject ref. [26]. 

1.3 Outline 

This work is dedicated to the holographic study of fundamental flavor. We focus on the 
D7-brane probing of asymptotically AdS^ x gravitational backgrounds. In particular 
we concentrate on meson spectroscopy, description of non-perturbative phenomena 
such as chiral symmetry breaking, and study of the thermal properties of the theory. 

In Chapter 2, we consider a D7-brane probe of AdS5 x in the presence of a 
pure gauge S-field. In the dual gauge theory, the B-field couples to the fundamental 
matter introduced by the D7-brane and acts as an external magnetic field. The B- 
field supports a 6-form Ramond-Ramond potential on the D7-branes world volume 
that breaks the supersymmetry and enables the dual gauge theory to develop a non- 
zero quark condensate. We explore the dependence of the quark condensate on the 
bare quark mass niq and show that at zero bare quark mass a chiral symmetry is 
spontaneously broken. We also study the discrete self-similar behavior of the theory 
near the origin of the parameter space given by the bare quark mass and the condensate 
of the theory. We calculate the critical exponents of the bare quark mass and the 
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quark condensate. A study of the meson spectrum supports the expectation based on 
thermodynamic considerations that at zero bare quark mass the stable phase of the 
theory is a chiral symmetry breaking one. Our study reveals a self-similar structure 
of the spectrum near the critical phase of the theory, characterized by zero quark 
condensate. We calculate the corresponding critical exponent of the meson spectrum. 
A further study of the meson spectrum reveals a coupling between the vector and 
scalar modes, and in the limit of weak magnetic field we observe a Zeeman splitting of 
the energy levels. We also observe the characteristic ^rriq dependence of the ground 
state corresponding to the Goldstone boson of spontaneously broken chiral symmetry. 

Chapter 3 is studying the finite temperature dynamics of flavored large AT^, SU{Nc) 
gauge theory with fundamental quark flavors in the quenched approximation. A quark 
condensate forms at finite quark mass, and the value of the condensate varies smoothly 
with the quark mass for generic regions in parameter space. At a particular value of 
the quark mass, there is a finite discontinuity in the vacuum expectation value of the 
condensate, corresponding to a first order phase transition. We study the gauge theory 
via its string dual formulation using the AdS/CFT conjecture, the string dual being 
the near-horizon geometry of Nc D3-branes at finite temperature, AdSs-Schwarzschild 
xS^, probed by a D7-brane. The D7-brane has topology x x and allowed 
solutions correspond to either the or the shrinking away in the interior of the 
geometry. The phase transition represents a jump between branches of solutions having 
these two distinct D-brane topologies. 

In Chapter 4, using a ten dimensional dual string background, we study aspects 
of the physics of finite temperature large Nc four-dimensional SU{Nc) gauge theory. 
We focus on the dynamics of fundamental quarks in the presence of a background 
magnetic field. At vanishing temperature and magnetic field, the theory has Af = 2 
supersymmetry and the quarks are in hypermultiplet representations. In Chapter 2, 
similar techniques were used to show that the quark dynamics exhibit spontaneous 
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chiral symmetry breaking. In this chapter we begin by estabhshing the non-trivial 
phase structure that results from finite temperature. We observe, for example, that 
above the critical value of the field that generates a chiral condensate spontaneously, 
the meson melting transition disappears, leaving only a discrete spectrum of mesons at 
any temperature. We also compute several thermodynamic properties of the plasma 
such as the free energy, entropy and magnetization. The study of the meson spectrum 
allows us to examine the stability of the theory. Our study shows that for sufficiently 
strong magnetic field, bellow the critical value, there is a metastable chiral symmetry- 
breaking phase that eventually becomes the true stable phase of the theory. We also 
obtain the meson spectrum of the coupled vector and scalar modes at strong magnetic 
field and verify the Zeeman splitting of the spectrum. 

In Chapter 5 we perform a study analogous to the one in Chapter 4 but for the 
case of external electric field. At zero temperature, we observe that the electric field 
induces a phase transition associated with the dissociation of the mesons into their 
constituent quarks. This is an analogue of an insulator-metal transition, since the 
system goes from being an insulator with zero current (in the applied field) to a con- 
ductor with free charge carriers (the quarks). At finite temperature this phenomenon 
persists, with the dissociation transition become subsumed into the more familiar me- 
son melting transition. Here, the dissociation phenomenon reduces the critical melting 
temperature. We also focus on the geometric aspect of the transition by performing 
an analogous T-dual study showing that the nature of the instability leading to the 
observed insulator/conductor phase transition is related to the over spinning of the 
T-dual D6-branes. We conclude the chapter with the identification of a peculiar class 
of embeddings that have a conical singularity above the horizon. We propose that 
those are fixed by stringy corrections. 

Chapter 6 is dedicated to the effect of an R-charge chemical potential on the 
meson melting phase transition. We begin by introducing the relevant gravitational 
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background which is that of spinning D3-branes. We comment on the deformation 
of the internal symmetries due to the three angular momentums of the geometry and 
identify ansatz relevant for the introduction of a probe D7-brane. We also dedicate 
a special section to the detailed calculation of the temperature and elaborate on the 
existence of extremal horizons in the case of multiple R-charges that would enable us 
to study the theory at zero temperature and non-zero R-charge chemical potential. 
In the second part of the chapter we study the properties of the D7-brane embeddings 
for two different cases. First we consider the case of only one non-zero R-charge. We 
calculate the equation of state in the quark condensate versus bare quark mass plane. 
We use this to study the dependence of the critical mass on the R-charge and develop 
analytic expression for the bare quark mass at large R-charge. Finally we obtain 
the phase diagram of the theory using appropriate dimensionless parameters which 
at fixed bare quark mass and t'Hooft coupling are proportional to the temperature 
and the R-charge chemical potential of the theory. After that we consider the more 
interesting case of three equal charges. Here we perform analysis analogous to the one 
for the one charge case and again obtain the phase diagram of the theory in coordinates 
corresponding to the temperature and the chemical potential. Our study shows that 
the phase diagram is similar to the phase diagram for the Hawking- Page transition of 
the adjoint fields. In the last section of the chapter we discuss the physical meaning of 
the constant of integration obtained while regularizing the behavior of the D7-brane 
probes at the ergosphere of the background. We conclude that this constant is related 
to the phase difference between the bare quark mass and the quark condensate. 

In Chapter 7 we show how two important types of phase transition in large Nc gauge 
theory with fundamental flavors can be cast into the same classifying framework as the 
meson-melting phase transition. These are quantum fluctuation induced transitions 
in the presence of an external electric field, or a chemical potential for R-charge. The 
classifying framework involves the study of the local geometry of a special D-brane 
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embedding which seeds a self-similar spiral structure in the space of embeddings. The 
properties of this spiral, characterized by a pair of numbers, capture some key universal 
features of the transition. Computing these numbers for these non-thermal cases, we 
find that these transitions are in the same universality class as each other, but have 
different universal features from the thermal case. The phase transitions that we 
consider are the thermal studied in Chapter 3, electrically driven studied in Chapter 5 
and the R-charge driven phase transition considered in Chapter 6. We present a 
natural generalization that yields new universality classes that may pertain to other 
types of transition. 

Chapter 8 is the concluding chapter. It summarizes our results and outlines possible 
directions for future studies. 
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Chapter 2: Flavored large Nc gauge theory in an 
external magnetic field at zero temperature 

2.1 Introductory remarks 

In recent years progress has been made towards the study of matter in fundamental 
representation in the context of AdS/CFT correspondence. As we discussed in Chap- 
ter 1, one way to achieve this is by introducing space filling fiavor D7-branes in the 
probe limit [50] and in order to keep the probe limit valid the condition A^^ <C Nc 
is imposed. The fundamental strings stretched between the stack of Nc D3-branes 
and the flavor Nj D7-branes give rise to an Af=2 hypermultiplct, the separation of 
the D3- and D7- branes in the transverse directions corresponds to the mass of the 
hypermultiplct, the classical shape of the D7-brane encodes the value of the fermionic 
condensate, and its quantum fluctuations describe the light meson spectrum of the the- 
ory [59] . This technique for introducing fundamental matter has been widely employed 
in diflFerent backgrounds. Of particular interest was the study of non-supersymmetric 
backgrounds and phenomena such as spontaneous chiral symmetry breaking. These 
phenomena were first studied in this context in ref. [12], where the authors developed 
an appropriate numerical technique. In recent years this approach received further de- 
velopment, and has proven itself as powerful tool for the exploration of confining gauge 
theories. In particular, for the description of their thermodynamic properties and for 
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the building of phenomenological models relevant to QCD. The chapter is organized 
as follows: 

In the second section we describe the method of introducing magnetic field to the 
theory, employed in ref. [29]. We describe the basic properties of the D7-brane em- 
bedding and the thermodynamic properties of the dual gauge theory, in particular the 
dependence of the quark condensate on the bare quark mass. We describe the spon- 
taneous chiral symmetry breaking caused by the external magnetic field and comment 
on the spiral structure in the condensate versus bare quark mass diagram. We perform 
analysis similar to the one considered in ref. [32] for the study of merger transitions 
and calculate the scaling exponents of the bare quark mass and the quark condensate 
[28]. We also describe the discrete self-similarity of the spiral and calculate the scaling 
factor characterizing it. 

In the third section we study the light meson spectrum of the dual gauge theory. 
First we derive the relevant equations of motion for the scalar and vector meson spec- 
trum. The study of the fluctuations along the axial scalar reveals a Zeeman slitting of 
the energy levels at weak magnetic field and a characteristic Gell-Mann-Oakes-Renner 
relation [34] for the pion of the softly broken chiral symmetry. 

Next we consider the meson spectrum of the states corresponding to the spiral. We 
study the critical embedding corresponding to the center of the spiral and reveal an 
infinite tower of tachyonic states organized in a decreasing geometrical series. After 
that we consider the dependence of the meson spectrum on the bare quark mass and 
confirm the expectations based on thermodynamic considerations that only the lowest 
branch of the spiral is stable. We observe that at each turn of the spiral there is 
one new tachyonic state. We comment on the self-similar structure of the spectrum 
and calculate the scaling exponent of the meson mass. We also consider the spectrum 
corresponding to the lowest branch of the spiral and for a large bare quark mass 
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reproduce the result for pure M = 2 Supersymmetric Yang Mills Theory obtained in 
ref. [59]. 

We end with a short discussion of our results and the possible directions of future 
study. 

2.2 Fundamental matter in an external magnetic field 
2.2.1 Basic configuration 

Let us consider the AdSs x geometry describing the near-horizon physics of a col- 
lection of Nc extremal D3-branes. 

2 ^2 

ds"^ = "^{-dxl + dxj + dxl + dxD + R^^ + R^dnl , (21) 

«4 



9sC{4) = ^ dx^ ^ dx^ ^ dx^ 5 



3 = 9s , 



Where dCl'l is the unit metric on a round . In order to introduce fundamental matter 
we first rewrite the metric in the following form, with dft^ the metric on a unit S^: 

ds^ = ^—^2 — [-dxl + dxj + dxl + dxl\ + 2 ^ ]^2 l-^P^ + p^d^l + dL^ + L^d(j?\ , 
dnl = dV'^ + cos^ + sin^ V'd7^ (22) 



where p, ■0, /3, 7 and L, 4> are polar coordinates in the transverse and respectively. 
Note that: u'^ = + L^. We use xo,xi, X2, x^, p,ip, P,'y to parameterize the world 
volume of the D7-brane and consider the following ansatz [59] for its embedding: 

(f) = const, L = L{p), 
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leading to the following form of the induced metric on its world-volume: 

ds = ^^J'^^\ -dxl + dxi + dxl + dxl] + _^^^^^^^ [{l + L'{pf)dp^ + p^dnl] . (23) 

Now let us consider the general DBI action (to simplify notations we temporary set 
the number of flavor branes to one Nf = 1): 

Sdbi = -fi7 J d^Ce-^i-detiGab + Bab + 27Ta'Fab)]^^^ ■ (24) 
Ms 

Here //y = [(27r)^a'*^]~^ is the D7-brane tension, Gab and Bab are the induced metric 
and S-field on the D7-brane's world volume, while Fab is its world-volume gauge field. 
A simple way to introduce magnetic field would be to consider a pure gauge 5-field 
along parts of the D3-branes' world volume, e.g.: 

5(2) = Hdx2 A dx2 . (25) 



Since Bab can be mixed with the gauge field strength Fab, this is equivalent to a 
magnetic field on the world-volume [29]. Recently a similar approach was used to 
study drag force in SYM plasma [69]. Note that since the B-field is pure gauge, 
dB = 0, the corresponding background is still a solution to the supergravity equations 
of motion. On the other hand, the gauge field Fab comes at next order in the a' 
expansion compared to the metric and the S-field components. Therefore to study 
the classical embedding of the D-brane one can study only the {Gab + Bab) part of the 
DBI-action. However, because of the presence of the S-field, there will be terms of 
first order in a' in the full action linear in the gauge field Fab- Hence integrating out 
Fab will result in a constraint for the classical embedding of the D7-brane. 
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Since for our configuration, we have that: 

5(2) ^ ^(2) ^ Q ^ ^(2) ^^^^ ^ Q ^ 

and at first order in a' the only contribution to the Wess-Zummino is 

27raV7 J F^2) A C(6) . (26) 
By using the following expansion in the DBI action: 

[-det{Eab + 27ra'F„6)] V2 = + na'^EE^^'Fab + 0{F^y, E = -detEab , (27) 

where we have introduced Eab = Gab + Bab as a notation for the generalized induced 
metric, we obtain the following action to first order in a': 

Sf = Tra'^ y" d^^VEE^-'^^F^.^ + 27ra' in J F(2) A C(6) • (28) 
^ Ms 

The resulting equation of motion does not contain Aa and sets the following constraint 
for the C(6) potential induced by the gauge 5-field. 

|^a6Ai...A6 5^C'/i,...^e = -daiVEE^'^^); a,b,h, ■ ■ ■ h & Ms ■ (29) 

Note that C(q) has a dynamical term proportional to 1 / Kq in the supergravity action, 
and that the D7-brane action is proportional to /X7 = Itt/kq. Therefore they are at 
the same order in a' [49]. We must solve for C(e) using the action: 

Sc,,, =1^7 J 5(2) A C(6) - ^ / c/^'^xV^ldqe)!' . (30) 
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The solution obtained from equation (30) has to satisfy the constraint given in equa- 
tion (29). Our next goal will be to find a consistent ansatz for C(g). To do this let us 
consider the classical contribution to the DBI action: 



Sns = -- d^'^VE . (31) 



Prom equation (31) one can solve for the classical embedding of the D7-brane, which 
amounts to second order differential equation for L{p) with some appropriate solution 
Lq{p). After substituting Lo(p) in (31) we can extract the form of the C(6) potential 
induced by the field. However, one still has to satisfy the constraint (29). It can be 
verified that with the choice (25) for the S-field and the ansatz of equation (23) for the 
induced metric, the right-hand side of equation (29) is zero. Then equation (29) and 
the effective action (30) boil down to finding a consistent ansatz for C^g) satisfying: 

5^(V=GdC6^°i'"^"^) =-^^H5{L-Lo{p)) , (32) 
d.i^dCf''^^''^) =-L',{p)^H5{L-Lo{p)) . (33) 
^afeAi-./iea^C^^ = 0; a, 6, /ii, . . . /ie £ • (34) 

One can verify that the choice: 



is a consistent ansatz and the solution for the C^g) field strength can be found to be: 

^^(6)^o« = ^^Z(^^^©(^-^o(p))sinV^cos^ . (36) 

It is this potential which breaks the supersymmetry. It is important to note that 
there is no contradiction between the fact that the i?-field that we have chosen does 
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not break the supersymmctry of the AdSs x supergravity background, on the one 
hand, and the fact that the physics of the D7-brane probing that background does 
have supersymmetry broken by the 5-field, on the other. This is because the physics 
of the probe does not back-react on the geometry. 

In what follows, we will study the physics of the D7-branes and the resulting dual 
gauge theory physics. Among the solutions for the D7-brane embedding, there will 
be a class with non-trivial profile having zero asymptotic separation between the D3- 
and D7-branes. This corresponds to a non-zero quark condensate at zero bare quark 
mass. Therefore the non-zero background magnetic field will spontaneously break the 
chiral symmetry. Geometrically this corresponds to breaking of the SO (2) rotational 
symmetry in the (L, -plane [59]. 

2.2.2 Properties of the solution 

We now proceed with the exploration of the properties of the classical D7-brane em- 
bedding. If we consider the action (31) at leading order in a', we get the following 
effective lagrangian: 



C = -f^p^sm^Pcos2l;Vl + L'\ 1+ . f'^' ■ (37) 

gs V (p^ + L^r 



The equation of motion for the profile Lo{p) of the D7-brane is given by: 



As expected for large (Lq + p^) — oo or — 0, we get the equation for the pure 
AdSs X background [50]: 
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Therefore the solutions to equation (38) have the following behavior at infinity: 



Lo{p) = m+ + .... (39) 

where the parameters m (the asymptotic separation of the D7- and D3- branes) and 
c (the degree of bending of the D7-brane) are related to the bare quark mass niq = 
m/2'Ka' and the quark condensate — c respectively [60]. As we shall see below, 

the presence of the external magnetic field and its effect on the dual SYM provide a 
non vanishing value for the quark condensate. Furthermore, the theory exhibits chiral 
symmetry breaking. 

Now notice that H enters in (37) only through the combination H^R^. The other 
natural scale is the asymptotic separation m. It turns out that different physical 
configurations can be studied in terms of the ratio = m? / {HR^): Once the fh 
dependence of our solutions are known, the m and H dependence follows. Indeed, let 
us introduce dimensionless variables via: 

p = rVHp , Lo = rVhL , L'oip) = L'{p) . (40) 

The equation of motion (38) then takes the form: 



d3\p^^^^=^l + ^ — + / ^ — = (41) 



The solutions for L{p) can be expanded again to: 



Hp) = m + ^ + . . . , (42) 
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and using the transformation (40) we can get: 



c = cR^H^/^ . (43) 

It is instructive to study first the properties of (41) for m S> 1, which corresponds 
to weak magnetic field H <C m?/B?, or equivalently large quark mass m S> Ry/H. 

2.2.3 Weak magnetic field 

In order to analyze the case of weak magnetic field let us expand L(p) = m + r/(p) and 
linearize equation (41) while leaving only the leading terms in (p^ + m^)"^. The result 
is: 

which has the general solution: 

r?(p) = . -T^ -2^ + • (45) 

From the definition of ri{p) and equation (42) we can see that C\ = c and since 
jyl/o-^oo = we have C2 = 0. Now if we consider m large enough, equation (45) should 
be valid for all p. It turns out that if we require that our solution be finite as p — >^ 
we can determine the large m behavior of c. Indeed, the second term in (45) has the 
expansion: 

1 1 1 9^ / ^ 

- = -7^^ + 7^ +0 P^) . 46) 

4p2(m2 + p2) 4m p^ 4m'^ 
Therefore we deduce that: 
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and finally, we get for the profile of the D7-brane for m ^ 1: 

" 11 TTh 

L{p) = rh+ —— - . (48) 

If we go back to dimensionful parameters we can see, using equations (43) and (47) 
that for weak magnetic field H the theory has developed a quark condensate: 

M) oc-c=-—H^ . (49) 
4m 

However, this formula is valid only for sufficiently large m and we cannot make any 
prediction for the value of the quark condensate at zero quark mass. To go further, 
the involved form of equation (41) suggests the use of numerical techniques. 

2.2.4 Numerical results 

In this subsection wc solve numerically equation (41) for the embedding of the D7- 
brane, using Mathematica. It is convenient to use initial conditions in the IR. We 
use the boundary condition L'{p)\p=o = 0. We used shooting techniques to generate 
the embedding of the D7-brane for a wide range of rh. Having done so we expanded 
numerically the solutions for L(p) as in equation (42) and generated the points in 
the (m, — c) plane corresponding to the solutions. The resulting plot is presented in 
Figure 2. 

As one can see there is a non zero quark condensate for zero bare quark mass, 
the corresponding value of the condensate is Ccr = 0.226. It is also evident that the 
analytical expression for the condensate (47) that we got in the previous section is 
valid for large m, as expected. Now using equation (43) we can deduce the dependence 
of Ccr on H: 

= CcrR^H^/^ = 0.226R^H^/'^ . (50) 
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-c(m) 




Figure 2: The black line corresponds to (47), one can observe that the analytic result 
is valid for large rh. It is also evident that for m = {ipip) ^ 0. The corresponding 
value of the condensate is Ccr = 0.226. 

It is interesting to check the consistency of our numerical analysis by solving equa- 
tion (38) numerically and extracting the value of Ccr for wide range of R^H, the re- 
sulting plot fitted with equation (50) is presented in Figure 3. 

Another interesting feature of our phase diagram is the spiral behavior near the 
origin of the (m, — c)-plane which can be seen in Figure 5. Note that the spiral presented 
in Figure 5 has two arms, we have used the fact that any two points in the (m, — c) 
plane related by reflection with respect to the origin describe the same physical state. 
A similar spiraling feature has been observed in ref. [3] , where the authors have argued 
that only the lowest branch of the spiral corresponding to positive values of m is 
the stable one (corresponding to the lowest energy state). The spiral behavior near 
the origin signals instability of the embedding corresponding to Lq = 0. If we trace 
the curve of the diagram in Figure 5 starting from large m, as we go to smaller 
values of m we will reach zero bare quark mass for some large negative value of the 
quark condensate Ccr- Now if we continue tracing along the diagram one can verify 
numerically that all other points correspond to embeddings of the D7-brane which 
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Figure 3: A plot of the magnitude of the quark condensate at zero bare quark mass 
Ccr as function of R^H, the black curve represents equation (50). 

intersect the origin of the transverse plane at least once. After further study of the 
right arm of the spiral, one finds that the part of the diagram corresponding to negative 
values of m represents solutions for the D7-brane embedding which intersect the origin 
of the transverse plane odd number of times, while the positive part of the spiral 
represents solutions which intersect the origin of the transverse plane even number 
of times. The lowest positive branch corresponds to solutions which don't intersect 
the origin of the transverse plane and is the stable one, while the upper branches 
have correspondingly 2,4, etc., intersection points and are ruled out after evaluation 
of the free energy. Indeed, let us explore the stability of the spiral by calculating the 
regularized free energy of the system. We identify the free energy of the dual gauge 
theory [4, 25] with the wick rotated and regularized on-shell action of the D7-brane: 



(51) 




(52) 
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The second term under the sign of the integral in (52), corresponds to the subtracted 
free energy of the L{p) = embedding and serves as a regulator. Now we can eval- 
uate numerically the integral in (52) for the first several branches of the spiral. The 
corresponding plot is presented in Figure 4. Note that we have plotted Idi versus |m|, 
since the bare quark mass depends only on the absolute value of the parameter m. 
The lowest curve on the plot corresponds to the lowest positive branch of the spiral, 
as one can see it has the lowest energy and thus corresponds to the stable phase of the 
theory. 




Figure 4: The lowest lying curve correspond to the positive fh part of the lowest branch 
of the spiral, suggesting that this is the stable phase of the theory. 

In the next section we will provide more detailed analysis of the spiral structure 
from Figure 5 and explore the discrete self-similarity associated to it. 

2.2.5 Criticality and spontaneous chiral symmetry breaking 

In the following section we analyze the spiral structure described in the previous sub- 
section. The technique that we employ is similar to the one used in ref. [32] and [67] , 
where the authors studied merger transitions in brane/black hole systems. 
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Figure 5: A magnification of Figure 2 to sliow the spiral behavior near the origin of 
the (— c, m)-plane. We have added the second (left) arm of the spiral representing the 
(m, — c) {—rn, c) symmetry of the diagram. 

Let us explore the asymptotic form of the equation of motion of the D7-brane 
probe (41) in the near horizon limit + !?■ ^ 0. To this end we change coordinates 
to: 



p —>■ \p\ L —> XL; 



(53) 



and consider the limit A ^ 0. The resulting equation of motion is: 





= . 



(54) 



+ L2 yrri^ 



(/52 + L2)2 



Equation (54) enjoys the scaling symmetry: 



p pp; L pL . 



(55) 
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In the sense that if L = f{p) is a solution to the E.O.M. then j^f{fip) is also a 
solution. Next we focus on the region of the parameter space, close to the trivial L = 
embedding, by considering the expansion: 



L = + (27ra')x 



(56) 



and linearizing the E.O.M. . The resulting equation of motion is: 



pdpipdpx) + 2x = 



(57) 



and has the solution : 



X = A cos( V 2 lnp) + B sin( V 2 In p) . 



(58) 



Now under the scaling symmetry p ^ p,p the constants of integration A and B trans- 
form as: 



V 



(59) 



cos In ij, sin -\/2 In p 
■ sin \/2 In p cos \/2 hipj 

The above transformation defines a class of solutions represented by a logarithmic 
spiral in the parameter space {A,B) generated by some (Aj„,Sj„), the fact that we 
have a discrete symmetry x ~^ ~X suggests that {—A^n, —Bin) is also a solution and 
therefore the curve of solutions in the parameter space is a double spiral symmetric 
with respect to the origin. Actually as we are going to show there is a linear map from 
the parameter space {A, B) to the plane (rn, — c), which explains the spiral structure, a 
subject of our study. To show this let us consider the linearized E.O.M. before taking 
the A — ^ limit : 



p^/l + f%{p^/l + f%x) + 2x = 0; X = Ax , 



(60) 
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with the solution: 



X = cos \/2 In 



Expanding at infinity: 



+ Bsm.V2ln 



(61) 



X 



rh + ^ + 
P 



n B 1 
V2P^ 



(62) 



we get: 



m 



V" 



A 

-B/V2 



(63) 



Now if we match our solution (61) with the solution in the p — limit (58) we 
should identify {A,B) with the parameters {A,B). Combining the rescaling property 

of (A, B) with the linear map to (m, — c) wc get that the embeddings close to the trivial 
embedding L = are represented in the {in, — c) plane by a double spiral defined via 
the transformation: 







f rh 









— \/2sin v^ln|U 
-4= sin \/2 In /i cos \/2 In ji 





(64) 



Note that the spiral is double because we have the symmetry (m, — c) — (— m, c). This 
implies that in order to have similar configurations at scales fii and //2 we should have: 



cos \/2 In |Lti = ± cos \/2 In /X2 



(65) 



and hence 



V21n^ 



-mr, 



(66) 



37 



which is equivalent to: 

Ml 

Therefore we obtain that the discrete self-similarity is described by a rescaling by a 
factor of: 

q = e-""!^ « 0.10845 . (68) 

This number will appear in the next subsection where we will study the meson spec- 
trum. As one may expect the meson spectrum also has a self-similar structure. 

It is interesting to confirm numerically the self-similar structure of the spiral and to 
calculate the scaling exponents of the bare quark mass and the quark condensate. It is 
convenient to use the separation of the D3- and D7-branes at p = 0, Lin = L{0) as an 
order parameter. There is a discrete set of initial separations Lj„, corresponding to the 
points Hq, Hi, H2, ... in Figure 5 for which the corresponding D7-brane embeddings 
asymptote to fh — — as p — > 00. The trivial L = embedding has Lin — and 
is the only one which has a zero quark condensate (c = 0), the rest of the states have 
a non zero c and hence a chiral symmetry is spontaneously broken. Each such point 
determines separate branch of the spiral where c = c(m) is a single valued function. 
On the other side, each such branch has both positive m and negative m parts. The 
symmetry of the double spiral from Figure 5, suggests that the states with negative rh 
are equivalent to positive rh states but with an opposite sign of c. This implies that 
the positive and negative rh parts of each branch correspond to two different phases 
of the theory, with opposite signs of the condensate. As we can see from Figure 4 the 
lowest positive branch of the spiral has the lowest free energy and thus corresponds to 
the stable phase of the theory. In the next subsection we will analyze the stability of 
the spiral further by studying the light meson spectrum of the theory near the critical 
L = embedding. 
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Here we are going to show that both the bare quark mass rh and the quark con- 
densate c have critical exponent one, as L^n 0. Indeed, let us consider the scaling 
property (59), (64). If we start from some and transform to Lin = ^-^in' 
solve for /x and using equation (64) we can verify that the bare quark mass and the 
quark condensate approach zero linearly as — > 0. To verify numerically our analysis 
we generated plots of m/Lj„ vs. \/21og Lj„/27r and c/Li„ vs. \/21og Lj„/27r presented 
in Figure 6. 

The red curves in Figure 6 represent a fit with trigonometric functions of a unit 
period, as one can see the fit is very good as Lj„ — > 0. On the other side, for large Lj„ 
we obtain the results for a pure AdS^ x space, namely L = const, c = 0. It is also 
evident from the plots that the scaling exponents of m and c are equal to one. 

2.3 Meson spectrum 
2.3.1 General properties 

We study the scalar meson spectrum. To do so we will consider quadratic fluctuations 
[59] of the embedding of the D7-brane in the transverse {L, ^)-plane. It can be shown 
that because of the diagonal form of the metric the fluctuation modes along the cf) 
coordinate decouple from the one along L. However, because of the non-commutativity 
introduced by the i?-field we may expect the scalar fluctuations to couple to the vector 
fluctuations. This has been observed in ref. [10], where the authors considered the 
geometric dual to non-commutative super Yang Mills. In our case the mixing will be 
even stronger, because of the non-trivial profile for the D7-brane embedding, resulting 
from the broken supersymmetry. 
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Figure 6: The red curves represent fit with trigonometric functions of unit period. For 
smah Lin the fit is very good, while for large Lin we get the results for pure AdS^ x 
space, namely L = const, c = 0. The plots also verify that the scaling exponents of rh 
and c are equal to one. 
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Let us proceed with obtaining the action for the fluctuations. To obtain the con- 
tribution from the DBI part of the action we consider the expansion: 

L = Lo{p) + 2Tra'x, = + 27ra' , (69) 

where Lo{p) is the classical embedding of the D7-brane solution to equation (38). To 
second order in a' we have the following expression: 

Eab = El^ + 2TTa'El^ + {i-Ka'fEl^ , (70) 

where E^, E^^E"^ are given by: 

E^^^ = Gab{p,Lo{p),^) + Bab, 

^ab = {daXS', + d,xS'a) + ^LoCaOX + F^b (71) 

^ab = {daxdbX + LldM^) - ^f!;'% {dax6', + d^x^O X + \dl,G,bX^ . 

P -\- Ijq {P -\- -tvQj z 

Here Gab and B^b are the induced metric and B field on the D7-brane's world volume. 
Now we can substitute equation (71) into equation (31) and expand to second order 
in a'. It is convenient [10] to introduce the following matrices: 

11^0,11-1 = 5 + J, (72) 
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where S is diagonal and J is antisymmetric: 



\\s-'\\ 

jab 
Gil 

Grvry 



= diag{-Gr/,Gri\ 
H 



G 



11 



Gf 1 + m 



t xct xb xb xct\ 



Gil 

^^pp jG,^,Gg,^,Gg0} , (73) 



G„ 



,2(1 + ^0 



-,2 j.ry . ri _ „^„2 



G 



^2 J- r-^ 



cos VG^^ ; G/3/3 = sin i\)G^^ . 



(74) 



(75) 



Now it is straightforward to get the effective action. At first order in a' the action 
for the scalar fluctuations is the first variation of the classical action (31) and is satisfied 
by the classical equations of motion. The equation of motion for the gauge field at first 
order was considered in Section 2.2 for the computation of the G(6) potential induced 
by the B- field. Therefore we focus on the second order contribution from the DBI 
action. 

After integrating by parts and taking advantage of the Bianchi identities for the 
gauge field, we end up with the following terms. For 



12 
2A 



(76) 



and for F: 



(77) 



and the mixed x~^ terms: 



sin2V' 



(78) 
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and for 

?2r 2 



(X l^/^J^^S-'da^d,^ , (79) 



where the function / in (78) is given by: 



f{p) = dp ( 5(p)-^i^ J23 ) + J^-'SlM + 25(/9) J^^^S^^a^oGn , (80) 



with 5(p) =^^^- =p^^i + Lq>^Ji + 



2\2 • 



sin •!/; COS V (p^ + Lq) 

As can be seen from equation (78) the A2, A3 components of the gauge field couple 
to the scalar field % via the function /. Note that since for p ^ 00 and L ^ 00, we 
see that J^^ — > 0, the mixing of the scalar and vector field decouples asymptoticly. In 
order to proceed with the analysis we need to take into account the contribution from 
the Wess-Zumino part of the action. The relevant terms to second order in a' are [10]: 

Swz = ^^^^M7 J i^(2) A F(2) A C(4) + {2^Ta')f^^ J F(2) A 5(2) A P[C(4)] , (81) 

where C(4) is the background R-R potential given in equation (21) and (7^4) is the pull 
back of its magnetic dual. One can show that: 

C4 = — -—^ v^L^ sin ih cos ibdib A da A d(3 A dd) . (82) 

gs [p^ + L^r 

Writing cf) = 27ra'$ we write for the pull back P[C'(4)]: 

P[C'(4)] = - — ^^!^K{p)da^diP A da A dp A (83) 



where we have defined: 
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Now note that the S-field has components only along x'^ and x^, therefore dx" in 
equation (83) can be only dp,dx^ or dx^. This will determine the components of the 
gauge field which can mix with However, after integrating by parts and using the 
Bianchi identities one can get the following simple expression for the mixing term: 

-(2W)2^ I d^e^Hd^K^F,, , (85) 

resulting in the following contribution to the complete lagrangian: 

oc ^-^^HdpK^Foi . (86) 

Note that this means that only the ^-nd Ai components of the gauge field couple 
to the scalar field Next the contribution from the first term in (81) is given by: 

(2W)2^ I d'^^-^±^F^,F^e<^'<^ , (87) 

where the indices take values along the p, ip, a, (3 directions of the world volume. This 
will contribute to the equation of motion for Ap, A^, A^ and Ap which do not couple to 
the scalar fluctuations. In this section we will be interested in analyzing the spectrum 
of the scalar modes, therefore we will not be interested in the components of the gauge 
field transverse to the D3-branes world volume. However, although there are no sources 
for these components from the scalar fluctuations, they still couple to the components 
along the D3-branes as a result setting them to zero will impose constraints on the 
^0 • • • ^3- Indeed, from the equation of motion for the gauge field along the transverse 
direction one gets: 

3 

^S^''dbdaAa = 0, b = p,^,a,P, (88) 

a=0 



44 



(Here, no summation on repeated indices is intended.) However, the non-zero in- 
field explicitly breaks the Lorentz symmetry along the D3-branes' world volume. In 
particular we have: 



which suggests that we should impose: 

-^0^0 + diAi = , ^2^2 + ^3^3 = . (90) 

We will see that these constraints are consistent with the equations of motion for 
^0 • • • ^3- Indeed, with this constraint the equations of motion for ^ and A^^, fj, = 
. . . 3 are, for %: 

' ^ ( 77^1 + ^ + 7^°X + (91) 



and for $: 



9 \l+L'^ (p2 + L^)2 g 
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and finally for A^. 




p2 (p2 + L^)2 g 



(93) 



1— + / 2 I r2N2 °-^l + ^^0^ = 0' 

= 0, 

9 

+ U2X = 0. 
9 



We have defined: 

n = -dl + dl+ . (94) 

As one can see the spectrum splits into two independent components, namely the 
vector modes ^1 couple to the scalar fluctuations along while the vector modes 
A2,As couple to the scalar modes along x- However, it is possible to further simplify 
the equations of motion for the gauge field. Focusing on the equations of motion for 
^0 and Ai in equation (93), it is possible to rewrite them as: 

9^' ) + + Ip^TW' " + d,)^-0 (95) 

1^ / gdpi-doAo + diA,) \ ^ An,{-doAo + diAi) 
9 1 + J 

+ (^2|^°M0^O + 5iAi) = 0. 

Note that the first constraint in (90) trivially satisfies the second equation in (95). In 
this way we are left with the first equation in (95) . Similarly one can show that using 
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the second constraint in (90) the equations of motion in (93) for A2 and boil down 
to a single equation for F23: 



1. / gdpF^s \ 



+ ,2nt"^^"^- +^-i9i + di)x = 0.(97) 



Now let us proceed with a study of the fluctuations along 

2.3.2 Fluctuations along $ for a weak magnetic field 

To proceed, we have to take into account the Fqi component of the gauge field strength 
and solve the coupled equations of motion. Since the classical solution for the embed- 
ding of the D7-brane is known only numerically we have to rely again on numerics to 
study the meson spectrum. However, if we look at equation (38) we can see that the 
terms responsible for the non-trivial parts of the equation of motion are of order H^. 
On the other hand, the mixing of the scalar and vector modes due to the term (86) ap- 
pear at first order in H. Therefore it is possible to extract some non-trivial properties 
of the meson spectrum even at linear order in H and as it turns out, we can observe 
a Zeeman-like eflFect: A splitting of states that is proportional to the magnitude of 
the magnetic field. To describe this, let us study the approximation of weak magnetic 
field. 

To first order in H the classical solution for the D7-brane profile is given by: 

Lo = m + 0{H^), (98) 
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where m is the asymptotic separation of the D3- and D7-branes and corresponds to 
the bare quark mass. In this approximation the expressions for g{p) and dpK{p), 
become: 



dip) = P , OpK{p) 



(p2 + m?)^ ' 

and the equations of motion for $ and FqIj equations (92) and (95), simphfy to: 

^(pWa,*) + ^4.+ j^;,^D#-4ff^^;,^F.. = 0, (99) 

where □ = -5^ + + a| + 5|, V'^ = -dl + 5^ . 

This system has become similar to the system studied in ref. [10] and in order to 
decouple it we can define the fields: 

^± = Foi ± mV^ , (100) 

where V = -sJ—dQ + di- The resulting equations of motion are: 

Note that is the Casimir operator in the (xq, x-\) plane only, while □ is the Casimir 
operator along the D3-branes' world volume. If we consider a plane wave e^^-^ then 
we can define: 

and we have the relation: 

= Mo\ - A;| - A;| . (102) 
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The corresponding spectrum of is continuous in k2,k3. However, if we restrict 
ourselves to motion in the (xq, a;i)-plane the spectrum is discrete. Indeed, let us 
consider the ansatz: 

</>± = ?7±(p)e-^^°'=o+^'=i^i _ (103) 

Then we can write: 



yMP%^^) + (p2 + 1,2)2 ^ ^(^^f^M±r?± = , (104) 



M± = Moi± . 



Let us analyze equation (104). It is convenient to introduce: 



y = -4; M± = — M±; P±(y) = (1 - y)"±v?±; (105) 



2a± = 1 + Jl + Ml- e = H- „ . 



With this change of variables equation (104) is equivalent to: 



y{l - y)P'l + 2(1 - (1 - a±)y)P' - a±(a±_i)P± ± e^^^|pP± = . (106) 



Next we can expand: 



P± = Po lb ePi + 0{e^) ; q± = oq ± ecti + 0{e^) ; (107) 

M± = Mo±eai^^°?>^ + 0(e') ; Mq = 2 Vao(ao + 1) • 
Mo 
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leading to the following equations for Pq and Pi: 



yil-y)P^ + 2{l-{l-ao)y)P^-aoiao-l)Po = 0, (108) 
y{l-y)P^ + 2{l-{l-ao)y)Pi-ao{ao-l)Pi = (ai(2ao - 1) 



The first equation in (108) is the hypergeometric equation and corresponds to the 
fluctuations in pure AdS^ x S^. It has the regular solution [59]: 

Poiy) = F{-ao,l-ao,2,y) . (109) 

Furthermore, regularity of the solution for ri{p) at infinity requires [59] that ao be 
discrete, and hence the spectrum of Mq: 

l-ao = -n, n = 0,l,... (110) 
Mo = 2V(n + l)(n + 2) . 

The second equation in (108) is an inhomogeneous hypergeometric equation. However, 
for the ground state, namely n = 0, Pq = -^(~1) 0, 2, y) = 1 and one can easily get the 
solution: 

P^{y) = ^ ln(l - y) + {6a, - Mo)(ln(-y) + ^) - . (Ill) 

On the other hand, using the definition of P±{y) in (105) to first order in e we can 
write: 

V± = (l T ln(l - ,)) (1 ± eP^iy)) , (112) 
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for the ground state oq = 1 and we end up with the following expression for r]±: 



1 Mo € , _ / 1 ln(l - y) 

V± = -, ±e—-^±- (6ai-Mo) ln(-y)+ ^ 



1-y 4(1 -y)2 1-y 



y 6 



(113) 



Now if we require that our solution is regular at y = and goes as 1/p^ oc 1/y at 
iufiiiity, the last term in (113) must vanish. Therefore we have: 



ai 



Mo 



(114) 



After substituting in (107) and (105) we end up with the following correction to the 
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Figure 7: Plot of M = MR^ /m vs. H/m for the first three states. The dashed black 
lines correspond to the spectrum given by equation (110). 



ground sate: 



M± = Mo ± 



H 



m 



(115) 



We observe how the introduction of an external magnetic field breaks the degeneracy 
of the spectrum given by equation (110) and results in Zeeman splitting of the energy 
states, proportional to the magnitude of H. Although equation (115) was derived 
using the ground state it is natural to expect that the same effect takes place for 
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higher excited states. To demonstrate this it is more convenient to employ numerical 
techniques for solving equation (104) and use the methods described in ref. [12] to 
extract the spectrum. The resulting plot is presented in Figure 7. As expected we 
observe Zeeman splitting of the higher excited states. It is interesting that equation 
(115) describes well not only the ground state, but also the first several excited states. 

It turns out that one can easily generalize equation (115) to the case of non-zero 
momentum in the (a;2, a;3)-plane. Indeed, if we start from equation (101) and proceed 
with the ansatz: 

ct>± = fi±ip)e-'--' , (116) 



we end up with: 



i,d,{p%fj±) + j^^^^Mlf,^ T H^^^^Moi^fj^ = , (117) 



Moi± = y M| + fc^3; k23 = yk^ + ki . 

After going through the steps described in equations (105)-(113), equation (114) gets 
modified to: 



Note that validity of the perturbative analysis suggests that ai is of the order of ao 
and therefore we can trust the above expression as long as /C23 is of the order of Mq. 
Now it is straightforward to obtain the correction to the spectrum: 



H L kj 



M, = Mo±^Wl + ^. (119) 



We see that the addition of momentum along the (0:2 — X3)-plane enhances the splitting 
of the states. Furthermore, the spectrum depends continuously on A;23. 
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2.3.3 Fluctuations along $ for a strong magnetic field 

For strong magnetic field we have to take into account terms of order H"^, which means 
that we no longer have an expression for Lq(p) in a closed form and we have to rely 
on numerical calculations only. Furthermore, there is no obvious way to decouple 
equations (92) and (95). However, it is still possible to extract information about the 
spectrum of the scalar modes if we restrict ourselves to fluctuations along the (x^, x^) 
plane. In this way there is no source term in equation (93), and we can consistently 
set Foi equal to zero. We consider time independent fluctuations satisfying the ansatz 
g-m23r-23^ (where r23 is the radial coordinate in the {x2 — X3)-plane). The damping 
factor in the exponent can be thought of as the mass of the scalar meson in 2 Euclidean 
dimensions. Indeed, let us consider the ansatz: 

$ = /i(p)e-*'=2^'-^*^3xV;(53) , (120) 

where Yi{S^) are the spherical harmonics on the sphere satisfying: Aq^YJ = + 
2)Yi. With this set-up the equation of motion for $, equation (92), reduces to equation 
for h{p): 

where we have defined: 

777.23 ~ ~ ^3 • (122) 

Before we proceed with the numerical analysis of equation (121) let us introduce di- 
mensionless variables by performing the transformation (40) and defining: 



R 

m23 = -j=m22. ■ (123) 
V H 
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The resulting equation is: 



1 + Z' (l + (^2 + L2)2)l/2 



1/2 



L2 



(p2 + ^2)2 + 1 



h{p) = . 



, , --d-nhip) + 



(124) 



In order to study the spectrum we look for normalizable solutions which have asymp- 
totic behavior h[p) oc l/p^ for large p and satisfy the following boundary conditions 
at /) = 0: 

h'{Q) = 0; /i(0) = 1 . (125) 

Let us consider first the lowest level of the spectrum. The spectrum that we get as a 
function of the bare quark mass is plotted in Figure 8. 

|S(e)| vs \og{ojt) 

|S(e)| 




log[<y/] 



Figure 8: Spectrum of 77123 vs. ifi. The dashed line represents the lowest level of the 
meson spectrum for pure AdSs x space. 

For large fh the spectrum asymptotes (the dashed line in Figure 8) to the one for 
pure AdS^ x space obtained by ref. [59] 



Mo = |^\/(n + / + l)(n + / + 2) 



(126) 
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Figure 9: Enlargement of part of the spectrum of 77123 vs m from Figure 8. The black 
solid curve shows the oc fit- 



with the substitution n = 0, ? = 0, to obtain our case. Therefore we are describing 
the lowest possible state of the meson spectrum. In figure 9 we have zoomed in the 
area near the origin of the (m, 'm23)-plane, one can see that for small values of m = 



mass rriq, which is to be expected since the chiral symmetry associated with the spinor 
representation of S0(2) is spontaneously broken [34]. 

It is interesting to look for modes corresponding to higher excited states (non- 
zero n). In Figure 10 we have presented a plot of some of these. Again, the dashed 



values of rh one can observe the qualitative difference of the behavior of the spectrum 
corresponding to the n = state from that of the higher excited states. Indeed, as 
m — > the n = states follow the \/m behavior plotted in Figure 9, while the excited 
states tend to some finite values at zero bare quark mass. The n = states merge into 
the Goldstone boson of the spontaneously broken chiral symmetry. 




2'Ka'mq/ R\/~H we observe oc y/fn^ dependence of the ground state on the bare quark 



line correspond to the pure AdS^ x 



spectrum given by (126) for I = 0. For small 
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Figure 10: Spectrum of 77123 vs fh for n = ... 4. The dashed hues represent the 
spectrum for AdS^ x space. 



2.3.4 The critical L = embedding 

In this section we study the L = embedding and in particular the spectrum of the 
fluctuations along the L coordinate. Let us go back to dimensionful coordinates and 
consider the following change of coordinates in the transverse space: 



p = ucose , (127) 
L = u sin 6 . 

In these coordinates the trivial embedding corresponds to ^ = and in order to study 
the quadratic fluctuations we perform the expansion: 

9 = + {27:a)6e{t,u) , (128) 
Se = e~'^^h{u) . (129) 



Note that in order to study the mass spectrum we restrict the D7~brane to fluctuate 
only in time. In a sense this corresponds to going to the rest frame. Note that due to 
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the presence of the magnetic field there is a couphng of the scalar spectrum to the vector 
one. However, for the fluctuations along 9 the coupling depends on the momenta in 
the (x2, X3) plane and this is why considering the rest frame is particularly convenient 

Our analysis follows closely the one considered in ref. [48], where the authors have 
calculated the quasinormal modes of the D7-brane embedding in the AdS-black hole 
background by imposing an in-going boundary condition at the horizon of the black 
hole. Our case is the T — > limit and the horizon is extremal. However, the ^ = 
embedding can still have quasinormal excitations with imaginary frequencies, corre- 
sponding to a real wave function so that there is no flux of particles falling into the 
zero temperature horizon. The resulting equation of motion is: 

It is convenient to introduce the following dimensionless quantities: 




(131) 



and make the substitution [48] 

M.) = ^ = ^ + ^(1^, (132) 

leading to the equation for the new variable f{z): 

f"{z) + {u'-V{z))f{z) = 0. (133) 
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Where the effective potential is equal to: 



V{z) 



3 (l + 3z4)(l 



4z2 



[l + z 



4^2 



(134) 



The potential in (134) goes as ^ for z ^ and as — foi' z ^ 00 and is presented in 
Figure 11. As it was discussed in ref. [48] if the potential gets negative the imaginary 
part of the frequency may become negative. Furthermore, the shape of the potential 
suggests that there might be bound states with a negative oj"^. To obtain the spectrum 
we look for regular solutions of (133) imposing an in-falling boundary condition at the 
horizon {z 00). 

V(z) 



0.2 



-0.2 



-0.4 



-0.6 




Figure 11: A plot of the effective potential V{z) given in equation (134). 

The asymptotic form of the equation of motion at z — > 00 is that of the harmonic 
oscillator: 

r(z)+a;2^(z)=0, (135) 



with the solutions e^*"^^, the in- falling boundary condition implies that we should 
choose the positive sign. In our case the corresponding spectrum turns out to be 
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tachyonic and hence the exponents are real. Therefore the in-faUing boundary condi- 
tion simply means that we have selected the regular solution at the horizon: 2; — ^ oo. 
We look for a solution of the form: 



f{z) = e+'^'S{z) . (136) 

The resulting equation of motion for S{z) is: 

(-3 - 6z^ + 9z^)S{z) + 4^2(1 + {2iujS'{z) + S"{z)) = . (137) 

Next we study numerically equation (137). After solving the asymptotic form of 
the equation at the Horizon, we impose the following boundary condition at z = 1/e, 
where e is a numerically small number typically e = 10~^ : 

S(l/e) = l-^; S'aA) = ^. (138) 

after that we explore the solution for a wide range of w = itoj. We look for regular 
solutions which have \S{e) \ « 0, this condition follows from the requirement that x oc 2;^ 
as z ^ 0. It turns out that regular solutions exist for a discrete set of positive uj <^1. 
The result for the first six modes that we obtained is presented in table 2. 

The data suggests that as — > the states organize in a decreasing geometrical 
series with a factor q ^ 0.1084. Up to four significant digits, this is the number from 
equation (68) which determines the period of the spiral. We can show this analytically. 
To this end let us consider the rescaling of the variables in equation (137) given by: 

z = Xz; Lu = uj/X; A — ^ 00 . (139) 
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Tl 


(n) 


(n) , (n-1) 


n 
u 






1 
1 


z.oyuz X lu 


u. iuyzo 


2 


3 1 348 X in~3 


10846 


3 


3.3995 X 10"^ 


0.108-15 


4 


3.6865 X 10-^ 


0.10844 


5 


3.9967 X 10-6 


0.10841 



Table 2: Numerical data for the first six quasi- normal modes of the critical embedding. 
The data suggests that as n ^ oo the states organize in a decreasing geometrical series. 

This is leading to: 

9S{z) + 4:Z^{2iLuS'{z) + S"{z)) + ©(A"^) = . (140) 

The solution consistent with the initial conditions at infinity (138) can be found to be: 

S{z) = ^^e~'^e-'^^V^H'.^{zuy, cb = itbi , (141) 

where H^^}- is the Hankel function of the first kind. Our next assumption is that in the 
ujj ^ limit, this asymptotic form of the equation describes well enough the spectrum. 
To quantize the spectrum we consider some zo = zq/X 1, where we have 1 ^ ^ A 
so that the simplified form of equation (140) is applicable and impose: 

^(^o) = . (142) 

Using that zu = izuj this boils down to: 

H^^iiu^izo) = . (143) 
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Now using that luiZq <C 1 for a sufficiently small coj, we can make the expansion: 



(144) 

where Ai and A2 are real numbers defined via: 



+ iA2 = --i(i/2)-*^r(iV2) . (145) 

TT 



This boils down to: 



cos(V21n(w7Zo) + 0) = 0; ^ = 7r/2 - arg(Ai + iA2) . (146) 
The first equation in (146) leads to: 

a;('^) = le-"^e-"^=a;JV, (147) 



suggesting that: 



g = e ^/2 « 0.10845 . (148) 



This is the number given in (68). Note that the value of zq is a free parameter that 
we can fix by matching equation (147) to the data in table 2. On the other side, S{z) 
given in equation (141) depends only on zu = iujz and therefore once we have fixed 
Zq we are left with a function of ooi which zeroes determine the spectrum, equation 
(142). It is interesting to compare it to the numerically obtained plot of \S{e)\ vs. a;/, 
that we have used to determine the spectrum numerically. The result is presented in 
Figure 12, where we have used the n = 3 entry from table 2 to fix zq. One can see the 
good agreement between the spectrum determined by equation (142), the red curve in 
Figure 12 and the numerically determined one, the dotted blue curve. 
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|S(e)| vs log(t^/) 

|S(e)| 




Figure 12: The dotted blue curve corresponds to the numerical solution of equation 
(137) while the thick red curve is the one determined by equation (142). The plots are 
scaled to match along the vertical axis. 

2.3.5 The spectrum near criticality 

In this section we study the light meson spectrum of the states forming the spiral 
structure in the (m, — c) plane, Figure 5. In particular we focus on the study of the 
fluctuations along L. The corresponding equation of motion is given in (97). The effect 
of the magnetic field H is to mix the vector and the meson parts of the spectrum. 
However, if we consider the rest frame by allowing the fluctuations to depend only 
on the time direction of the D3-branes' world volume, the equation of motion for the 
fluctuations along L decouple from the vector spectrum. To this end we expand: 

L = Lo(/o) + (2W)x(p,t) , (149) 
X = h[p) cos Mt . 
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Here Lq{p) is tlie profile of the D7-brane classical embedding. The resulting equation 
of motion for h{p) is: 



where g{p, Lq, L'^) = pW^ + Lo'^l + jff^^ . 



It is convenient to introduce the dimensionless variables: 



r h ~ Lq ^ p ~ MR 

h= — ; Lo= ^;p=^^;M = ^, 151) 

R^/H RVH RVH VH 



leading to: 

S,(SjiJy) + {Sj^^-g^+S,(fL^))h = 0, (152) 
with 9(p, Lo, it) = p'^lTIifyiT^^JZ^ . 

We study the normal modes of the D7-brane described by equation (152) by im- 
posing Neumann boundary conditions at p = 0. Since our analysis is numerical we 
solve the equation of motion (152) in terms of a power series for small p and impose 
the appropriate initial conditions for the numerical solution at p = e, where e is some 
very small number. In order to quantize the spectrum we look for numerical solutions 
which are normalizable and go as l/p^ at infinity. 

Let us study the dependence of the spectrum of M on the bare quark mass rh, for 
the states corresponding to the spiral structure from Figure 5. A plot of the spectrum 
of the first three excited states is presented in Figure 13. The classification of the states 
in terms of the quantum number n is justified because at large rh the equation of motion 
for the fluctuations asymptotes to the equation of motion for the pure AdS^ x space, 
considered in ref. [59], where the authors obtained the spectrum in a closed form. Note 
that the diagram has a left-right symmetry. This is because we plotted the spectrum 
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for both arms of the spiral in order to emphasize its self-similar structure, physically 
only one side of the diagram is sufficient. 



M 

n=2 ' n=2 




-0.30 

I 



Figure 13: A plot of the meson spectrum corresponding to the two arms of the spiral 
structure at the origin of the (m, — c) plane. The ground state (n = 0) becomes 
tachyonic for the inner branches of the spiral while only the lowest branch is a tachyon 
free one. The tachyon sector of the diagram reveals the self-similar structure of the 
spectrum. 

Let us trace the blue curve corresponding to the n = state starting from the 
right-hand side. As m decreases the mass of the meson decreases and at m = it 
has some non-zero value. This part of the diagram corresponds to the lowest positive 
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branch of the spiral from Figure 5 (the vicinity of point Hq). It is satisfying to see that 
the lowest positive m branch of the spiral is tachyon free and therefore stable under 
quantum fluctuations. Note that despite that the negative m part of the lowest branch 
has no tachyonic modes in its fluctuations along L, it has a higher free energy (as can 
be seen from Figure 4) and is thus at best metastable. 

One can also see that the spectrum drops to a zero and becomes tachyonic exactly 
at the point where we start exploring the upper branch of the spiral. This proves 
that all inner branches correspond to true instability of the theory and cannot be 
reached by super-cooling. As we go deeper into the spiral, the n = spectrum remains 
tachyonic and spirals to some critical value. The dashed line denoted by o;}'^'* in Figure 5 
corresponds to the first entry in table 2. As one can see this is the critical value 
approached by the spectrum. 

Now let us comment on the n = 1 , 2 levels of the spectrum represented by the red 
and green curves, respectively. As one can see the n = 1 spectrum becomes tachyonic 
when we reach the third branch of the spiral (the vicinity of point H2 in Figure 5) and 
after that follows the same pattern as the n = level, spiraling to the second entry Uj^^ 
from table 2. The n = 2 level has a similar behavior, but it becomes tachyonic at the 
next turn of the spiral and it approaches the next entry from table 2. Similar feature 
was reported in ref. [68] where the authors studied topology changing transitions. The 
above analysis suggests that at each turn of the spiral, there is one new tachyonic 
state appearing. It also suggests that the structure of the n-th level is similar to the 
structure of the n + 1-th level and in the n — > 00 limit this similarity becomes an exact 
discrete self-similarity. The last feature is apparent from the tachyonic sector of the 
diagram in the second plot in Figure 13, the blue, red and green curves are related by 
an approximate scaling symmetry, the analysis of the spectrum of the critical L = 
embedding suggests that this symmetry becomes exact in the n — > 00 limit with a 
scaling factor of q given in equation (68). 
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It is interesting to analyze the way the meson mass M approaches its critical value 
and compute the corresponding critical exponent. Let us denote the critical value of 
M by M* and consider the bare quark mass m as an order parameter, denoting its 
critical value by m*. We are interested in calculating the critical exponent a defined 
by: 

|M-M*| oc |m-m*|" . (153) 

We will provide a somewhat heuristic argument that a = 2 and will confirm this 
numerically. To begin with let us consider the energy density of the gauge theory E as 
a function of the bare quark mass m. Now let us consider a state close to the critical 
one, characterized by: 

M = + SM; m = rh^ + Srh- E = + 5E . (154) 

Next we assume that as we approach criticality the variation of E and M are pro- 
portional to the variation of the energy scale and hence 5E oc 5M. Therefore we 
have: 

^oc^occ, (155) 
om dm 

where c is the quark condensate. The second relation in (155) was argued in ref. [60]. 
In the previous section we argued that the critical exponent of the condensate is one 
and since the critical embedding has a zero condensate it follows that c (x\m — ?fi*|. 
Therefore we have: 

, , , 

- — oc a m — m* oc m — m* (156) 
om 

and hence a = 2. 

Now let us go back to Figure 13. As we discussed above, for each energy level n the 

(n) 

tachyonic spectrum spirals to the critical value ujj , corresponding to the center of the 
spiral. If we focus on the m = axis, we can see that for each level we have a tower 
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of tachyonic states at a zero bare quark mass, corresponding to the different branches 

~ (n) 

of the spiral. Let us denote by the imaginary part of the meson spectrum, 

corresponding to the A;-th tachyonic state of the n-th energy level, at a zero bare quark 
mass fh. As we go deeper into the spiral, k ^ oo and M^"^ m!^^\ the data in 
Figure 13 suggests that Mi"^ = Lo^p\ On the other side, if the meson spectrum has a 
critical exponent of two, one can show that for a large k: 



(n) 



(n) 



(n) 
k-1 



M, 



(n) 



(157) 



where q is given by equation (68). We can solve for M, 



in). 



Mi") = Mfc_i + ^ 



(158) 



Now assuming that for k = 1,2 the approximate geometrical series defined via 
(157) is already exact we calculate numerically m["'\M2^^ for the n = 0,1,2 levels 
and compare the value of M^"^ obtained by equation (158) to the first three entries in 
table 2. The results are presented in table 3. 



n 




Mi") 


Mi") 


..in) 
Uj 





2.7530 X 10-^ 


2.6460 X 10-1 


2.6447 X 10-1 


2.6448 X 10-1 


1 


3.0162 X 10-2 


2.8917 X 10-2 


2.8902 X 10-2 


2.8902 X 10-2 


2 


3.2715 X 10-^ 


3.1363 X 10-3 


3.1347 X 10-3 


3.1348 X 10-3 



~ (n) 

Table 3: Comparison between the critical value M* ' for the n-th energy level of 
the meson spectrum and the corresponding quasi-normal mode Wj") of the critical 
embedding. The data suggests that they match. 



One can see that up to four significant digits the critical value of the meson spec- 
trum is given by the imaginary part of the quasi normal modes presented in table 2. 
This supports the above argument that the meson spectrum has a critical exponent of 
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two. Another way to justify this, is to generate a plot of the meson spectrum similar to 
the one presented in Figure 6 for the bare quark mass rh and the fermionic spectrum 
c. Notice that M approaches criticality from above while the parameter rh oscillates 
around the critical value m^: = 0. This suggests to use M as an order parameter and 
to generate a plot of ifi/{M — M^,)"^ vs. \/21og \M — M*|/27r. Note that according to 
equation (157) the plot should represent periodic function of an unit period. The re- 
sulting plot for the n = level, using from table 3 as a critical value, is presented 
in Figure 14. 
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Figure 14: A plot of the bare quark mass meson vs. the meson spectrum, in an 
appropriate parameterization, determined by the scaling exponents of rh and M. The 
discrete self-similar structure of the spectrum is manifested by the periodicity of the 
plotted function. 



2.3.6 The stable branch of the spiral 

In this subsection we consider the spectrum corresponding to the states far from the 
origin of the (m, — c) which define the outermost branch of the spiral ending at point Hq 
from Figure 5. The fluctuations of the D7-brane corresponding to the massless scalar 
(j) were studied in Subsection 2.3.3 and some features consistent with the spontaneous 
chiral symmetry breaking, such as a characteristic ^/m behavior [34] were reported. 
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Here we complement the analysis by presenting the results for the fluctuations along 
the L coordinate. Since this is the massive field in the spontaneous chiral symmetry 
breaking scenario, we expect a \/ const + rh behavior of the meson spectrum for small 
values of m. Note that such a behavior simply means that the spectrum of the L 
fluctuations has a mass gap at zero bare quark mass and that the slope of the spectrum 
vs. the bare quark mass function is finite. It is satisfying that our results are in accord 
with this expectations. 

To obtain the spectrum, we solve numerically equation (152) imposing Neumann 
boundary conditions at p = 0. A plot of the first five energy levels is presented in 
Figure 15. As one can see at large m the spectrum approximates that of the pure 
M =2 Flavored Yang Mills theory given in equation (126) which we replicate here: 

Mo= --^V(" + ^ + !)(" + ^ + 2) . (159) 

Here I is the quantum number corresponding to the angular modes along the internal 
sphere wrapped by the D7-brane and is zero in our case. After introducing the 
dimensionless variables defined in (151), equation (159) boils down to: 

Mo = 2A/(ra + l)(ra + 2)m . (160) 

The black dashed lines in Figure 15 represent equation (160). The fact that the meson 
spectrum asymptotes to the one described by (160) justifies the use of the quantum 
number n to classify the meson spectrum. One can also see that as expected the 
spectrum at zero bare quark mass has a mass gap. 
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Figure 15: A plot of the meson spectrum corresponding to the stable branch of the 
spiral. The black dashed lines correspond to equation (160), one can see that for large 
fh the meson spectrum asymptotes to the result for pure AdS^ x space. One can 
also see that at zero bare quark mass rh there is a mass gap in the spectrum. 

2.4 Concluding remarks 

In this chapter we performed a detailed analysis of the spiral structure at the origin 
of the condensate versus bare quark mass diagram. We revealed the discrete self- 
similar behavior of the theory near criticality and calculated the corresponding critical 
exponents for the bare quark mass, the quark condensate, and the meson spectrum. 

Our study of the meson spectrum confirmed the expectations based on thermody- 
namic considerations that the lowest positive rh branch of the spiral corresponds to a 
stable phase of the theory and that the inner branches are real instabilities character- 
ized by a tachyonic ground state and cannot be reached by a supercooling. The lowest 
negative rh branch of the spiral is tachyon free and thus could be metastable. 

The supercooling mentioned above could be attempted by considering the finite 
temperature background, namely the AdS-Black hole geometry, in the presence of an 
external magnetic field. We could prepare the system in the phase corresponding to the 
trivial L = embedding and then take the T ^ limit. If some of the inner branches 
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of the spiral were metastable the theory could end up in the corresponding phase. The 
study of the finite temperature case is of a particular interest. Due to the additional 
scale introduced by the temperature, the theory has two dimensionless parameters and 
is described by a two dimensional phase diagram. The effect of the temperature is to 
restore the chiral symmetry and is competing with that of the external magnetic field. 
On the other side, the magnetic field affects the melting of the mesons [4]. We will 
come back to that in Chapter 4 of this work. 

2.5 Appendix: Calculating the condensate of the theory 

Let us consider the on-shell action: 



Pmax 



S = -27r\,Nf I dpp'Jl + -^^^V^Tl7W (161) 
* 

it diverges as pmax ^ oo. To rectify this we regularize the action by subtracting the 
action for the L = embedding: 

Pmax j 

Ssub f , L R^H"^ 1 



// H^H^ 1 
^WVI + -^^p'LaWpta. + R^H^+ (162) 



Pm 



This results to the following regularized action: 

Sreg\L-i Pmax] = S — Sgub ■ (163) 

On the other-side the hamiltonian density of the theory can be written as [60] : 

'H = J (femgQQ + Ho , (164) 
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where TCq is the mass independent part of the hamiltonian and Q, Q are the two chiral 
fields of the hypermultiplet in J\f = 1 notations. Now by making the identification: 



{H) = - Hm Sreg , (165) 

pmax 



and using equations (161)-(163) and the asymptotic of L{p) as pmax 



we can obtain: 



L(p) = m+^— + ... , (166) 

Pmax 



(i^) = (v;V) = -27ra' lim = (167) 

Oniq Pmax^OO dL 



hm ATr^a'jyNfpl^^^L'ipmax) = Sn^a'-nNfC . 

Pmax ^OO 



And hence we finally get: 



i^^) = - .o_.nL2 g • (168) 



2.6 Appendix: Numerical techniques 

We present the Mathematica 6 code that can be used to generate some of the main 
numerical results reported in Chapter 2. In fact variations of the packages bellow can 
generate most of the numerical results reported in this thesis and are quite represen- 
tative for our research. 

Let us begin by presenting the code used to solve numerically equation (38) and 
generate the — c versus rh plot in Figure 17. The first two lines clear the varibles that 
we use and define equation (38). The assignment for b, R, h = 1, is equivalent to the in- 
troduction of dimensionless variables. The variables mmin and mmax specify the range 
of initial conditions for L(0) that we explore. The rest of the code solves numerically 
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equation (41) and reads off the asymptotic behavior of L{p) at jj = Rhomax = 300 
which is our numerical approximation for oo. The main Mathematica's command in 
this code is ND Solve. The last few lines of the code generate the plot of — c versus rh 
presented in Figure 2. Modification of the code presented in Figure 17 can be used 
to generate virtually all of the equation of states c(m) reported in this thesis. One 
need only change the equation of motion and the way the initial conditions are im- 
posed. The last is sometimes non-trivial from physics point of view because of the 

Clear [b, R, h, L] 
Equation : = 

(2 r L[r] + (3 + + 3 L[r] ^ (3 + R* + 3 L[r]^ + L[r] ) L' [r] ) (l + L'[r]^) + 

r (r^ + L[r]^) (r* + R* + 2 L[r] ^ + L[r]*) L"[r] 

b = 1; 
R = 1; 
h = 1; 

mmln = . 0001; 

mmax = 1.6; 

n = 100; 

Array[einb, n] ; 

6 = 10*-8; 

Rhomax = 300; 

Array [p, {n, 2)]; 

Array [massO, {n, 2)]; 

For [i :=0;i = l, i<n + l, i + + , { 

mO = nnnin + (nmnax - nnnin) i / n; 

mas50[i, 1] = mO ; 

s = NDSolve[ {Equation == 0, L[6] == mO, L ' [6] == 10*-8}, 

L, {r, 6, Rhomax}, AccuracyGoal -> Infinity] ; 
l[r_] := L[r] /. First[s]; 
p[i, 1] =l[Rhomax]; 
massO [i, 2] = p[i, 1] ; 
p[i, 2] =1/21' [Rhomax] Rhomax*3; 

emb[i] =Plot[l[x], {x, e, 30), DisplayFunction -» Identity] ; } 

]; 

listl = Array [p, {n, 2}]; 

w = ListPlot[listl, BaseStyle -> {FontFamily -» "Times" , FontSize -> 24} , Joined -» True, AxesLabel -> 

["m" , "-c(m)"}, PlotStyle -» {Hue[0] , Thickness[0 . 003] } , DisplayFunction -» Identity ]; 
w2 = ParametricPlot[{l + h*2 / (41^3), -lh''2/(4 1)}, {1, (3h^2/4)^(l/4), 2.5}, 

PlotStyle -» { Thickness [ . 003] } , DisplayFunction -» Identity ] ; 
Show[w, DisplayFunction -> $DisplayFunction ] 

Figure 16: Mathematica 6 code used to generate the — c versus rh plot reported in 
Figure 2. 

variety of vanishing loci that we may have. However, numerically it results to the same 
procedure. 

Let us focus now on the analysis of the free energy presented in Figure 4. The 
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Equation = 2 r L[r] (l + L'[r]^) + 

(l + 3r*) L[r]= (3 (L'[r] +L'[r]') +rL"[r]) +3r^L[r]^ (3 (L'[r] +L'[r]') +rL"[r]) + 

L[r]« (3 (L'[r] +L'[r]^) +rL"[r]) + ((l + 3r*) L'[r] (l + L'[r]') +r (l + r^) L"[r]); 
n = 450; 

Array[mass, {2, n}]; 
Array [F, n] ; 
Rhomin = 10^-6; 
Rhomax = 100; 
Lmax = . 6 ; 
6 = 10^-6; 

Ssub = Integrate , {r, Rhomin, Rhomax}]; 

For ^i :=0; i = l, i<n + l, i + 
Lin = Lmax * i / n; 

si = NDSolve [ {Equation == 0, L[Rhomin] == Lin, L ' [Rhomin] == 0}, L, {r, Rhomin, Rhomax}] ; 

l[r_] : = L[r] / . First [si] ; 

mass[l, i] =l[Rhomax]; 

mass[2, i] = . 5 1 ' [Rhomax] * Rhomax^ 3; 



F[i] =NIntegrate[ r^ ll+ 1 + 

'[ A| (r=*l[r]^)' 

}]'• 



l'[r]^ 



{r, Rhomin, Rhomax} j - Ssub; 



ListPlot[Table[{abs[mass[l, i] ] , F[i]}, {i, 1, n}]. 
Joined -» True, BaseStyle -> {FontFamily ^ "Times" , FontSize -» 24} , 
AxesLabel ^ {" |m| ", "im"} , PlotStyle -» {Hue[ . 6] , Thickness [0 . 003] } ] 

Figure 17: Mathematica 6 code used to generate the free energy /d7 versus |m| plot 
reported in Figure 4. 

code has the same structure as the one presented in Figure 17. The main difference 
is the addition of an extra array F[i] which corresponds to the regularized on-shell 
action of the D7-brane for a particular value of the bare quark mass parameter rh. 
The two most important Mathematica commands are NDSolve and NIntegrate. Again 
variation of this code can be used to study the free energy for all different physical 
systems considered in this document. 

And finally let us provide the code used to demonstrate the Gell-Mann-Oaks- 
Renner relation 77123 y/m demonstrated in Figure 8. The code runs in two parts. 
The first part extracts the value of the initial condition L(0), corresponding to the D7- 
brane embedding which has zero separation at infinity {rh = 0) and does not intersect 
the p axis of the (p, L)-plane. This embedding corresponds to the critical quark con- 
densate Ccr ~ 0.226, given in equation (50), and breaks the rotational U{1) symmetry 
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in the transverse to the D7-brane plane. Thus, providing a holographic description of 
spontaneous symmetry breaking. This part of the code is reported in Figure 18. 

b = 1; 
R = 1; 

Equation : = 

(2 r R* L[r] + (3 + R^ + 3 L[r] ^ (3 + R* + 3 L[r]^ + L[r]^) ) L' [r] ) (l + L' [r]^) + 

r (r^ + L[r]^) (r* + R* + 2 L[r]^ + L[r]*) L"[r] 

mmin = 1 . 508 ; 
mmax = 1.55; 
n = 200; 
e = 10''-8; 
Rhomax = 150; 
h = 10; 

Array [p, {n, 2}]; 
Array [massO, {n, 2}]; 
flag = 1; 

For [i :=0; i = n, i>0, i = i- l, { 
mO = mmin + (mmax - mmin) (i/n)'^1.5; 
massO[i, 1] = mO; 

s = NDSolve[ {Equation == 0, L[e] ==mO, L ' [e] ==10^-8}, L, {r, e, Rhomax}]; 
l[r_] := L[r] /. First [s] ; 
massO[i, 2] =l[Rhomax]; 

If [ (massO[i, 2] > 0) && (flag == 1) , lin = i, {flag = 0; 

}];}] 

Print [ "ready Iin=", lin]; 
Print["mO=", massO[Iin, 1]]; 
lists = Array [massO, {n, 2}]; 
ListPlot[list3] ; 
lin 

massO[Iin, 2] 
massO[Iin-l, 2] 

Figure 18: Mathcmatica 6 code used to extract the initial condition 1/^(0) for the 
embedding corresponding to the critical quark condensate Ccr ~ 0.226. 

Here mmin and mmax specify the interval where we expect to find the initial 
condition Lcr(O). We need a very high precision because we are trying to extract the 
behavior of the meson spectrum for very small m. The main output of this code is the 
position lin in the array of initial conditions massO[i, 1] for which massO[Iin, 1] = Lcr{0)- 
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EquatlonPhl := — Sr ^ 9r f [r] 



^3 Vl + 1 ' [r] *2 / 1 



(r~2tl[rr2)-2 



l[r] *2 r*3 



/ . R'4 
-y + (r-2*l[r]' 



Vl + 1 ' [r] *2 



l[r] ^2 R^4 
(r*2 + l[r] *2) '-2 



ml 2 + 



(r»2*l[r]-2)-2 I 



f[r] 



R= 1; 
h = 1; 

Array [M, {n, 1}] ; 
M[Iin - 1] = .1; 

For [^i :=0; i = Iin, i<n + l, i + + , | 

mO = mmin + (mmax - nimin) (i / n) ^1.5; 

s = NDSolve[ {Equation - 0, L[6] == mO, L ' [e] ==10*-8}, L, {r, 6, Rhomax}]; 
l[r_] :=L[r] /. First[s]; 
ml = 0; 
K = 220; 

Clear [H] ; 
Array [H, {2, K}] ; 

Step = 20 (2) * .5 l[Rhomax] *10*-2; 
Mstart = M[l - 1] ; 

in2 = Mstart; 
flag = ; 

H[l, 1] = m2 - Step; 
H[2, 1] = 1; 

For [ j ; = 0; j = 2, flag == 0, j + + , 
{ 

H[l, j] = in2; 

nsol = NDSolve[ {EquatlonPhl = 0, f[6] =1, f ' [e] =0}, f, {r, e, Rhomax}]; 
F[r_] :=f[r] /. First [nsol] ; 
H[2, j] = Re [F [Rhomax] ] ; 

If[H[2, j-l]*H[2, j] >0, flag = 0, 

H[l, j] -H[l, j -1] 

j flag = 1; M[i] = H[l, j - 1] H[2, j - 1] ; Print["M=", M[i] ] ; Prlnt[ j] ; [ ; 



H[2, j] -H[2, j-1] 



■step;}]}] 



llstlZ = Table [{mass0[l, 2], M[l]}, {1, Iln + 1, n}]; 
Flt[llstlZ, {x*(l/2)}, X] 

graphl = ListPlot [llstlZ, Joined -» False^ AxesLabel -> {"m", "m23"|, 

PlotStyle -» {Hue[0 . 6] , Thickness[0 . 004] } , BaseStyle -» {FontFamlly -» "Times" , FontSlze -» 24} , 

Ticks -» {Automatic, Automatic) , DisplayFunction -» Identity] ; 
graph2 = Plot[Evaluate[Fit[listlZ, {x''{l/2)}, x] ] , {x, 0, massO[n, 2] + .01), 

PlotStyle -» {Thickness [0 . 0022] ) , DisplayFunction -» Identity] ; 
graph20 = Show[Graphics[Dashlng[ {0 . 03, 0.03)]], graph2, DisplayFunction -» Identity] 
Show[graphl, graph20, DisplayFunction -> $DlsplayFunctlon] 

Figure 19: Mathematica 6 code used to obtain the m23 cx: \/^ relation reported in 
Figure 8. 



In Figure 19 we have presented the second part of the code generating the 77123 
^ffli relation. The code uses the value of lin calculated in the part presented in Fig- 
ure 18 and the same array of initial conditions massO[i, 1] to scan only through em- 
beddings corresponding to positive ih. For each initial condition it solves the equation 
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of motion for the classical embedding (the variable, equation, defined in the code in 
Figure 18) of the D7-brane and uses to solve the equation for the fluctuations along 
(f) (given by the variable, EquationPhi, defined in the top line in Figure 19). The code 
then extracts the meson spectrum by requiring regularity of the numerical solution 
for the fluctuations at p = Rhomax. This is the purpose of the inner cycle which is 
varying the integer j. It is to be noted that the code presented in Figure 19, after minor 
modifications, can be applied to other studies of meson spectrum. The main part of 
the code is the j-cycle which verifies numerically the normalizability of the quadratic 
fluctuations and quantizes the spectrum. 
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Chapter 3: Thermal properties of large Nc flavored 
Yang— Mills theory 



3.1 Introductory remarks 

In this chapter we perform a holographic study of the thermal properties of large Nc 
flavored Yang-Mills theory. The presented material is based on a paper written in a 
very close collaboration with Tamecm Albash, Arnab Kundu and Clifford Johnson [6], 
similar result was also reported in ref. [67]. Furthermore, the main result was presented 
in the Ph.D. thesis of Ingo Kirsch [54]. This is why this chapter should be regarded 
as an introduction to the ideology of Chapters 4,5,6 and 7, where the reported results 
are more original. 

We study the geometry of AdSs-Schwarzschild xS^ which is the decoupled/near- 
horizon geometry of Nc D3-branes, where Nc is large and set by the (small, for relia- 
bility) curvature of the geometry. The physics of closed type IIB string theory in this 
background is dual to the physics of A/^ = 4 supersymmetric SU (Nc) gauge theory in 
four dimensions, with the supersymmetry broken by being at finite temperature [87]. 
The temperature is set by the horizon radius of the Schwarzschild black hole, as we 
will recall below. 

We introduce a D7-brane probe into the background. Four of the brane's eight 
world-volume directions are parallel with those of the D3-branes, and three of them 
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wrap an C S^. The remaining direction lies in the radial direction of the asymp- 
totically AdS5 geometry. 

Such a D3-D7 configuration controls the physics of the SU{Nc) gauge theory with 
a dynamical quark in the fundamental representation[50]. The configuration (at zero 
temperature) preserves M = 2 supersymmetry in £) = 4, and the quark is part of a 
hypermultiplet. Generically, we will be studying the physics at finite temperature, so 
supersymmetry will play no explicit role here. 

We are studying the D7-brane as a probe only, corresponding to taking the Nc 3> 
Nf limit, and therefore there is no backreaction on the background geometry. This is 
roughly analogous to the quenched approximation in lattice QCD. The quark mass and 
other flavor physics-such as the vacuum expectation value (vev) of a condensate and 
the spectrum of mesons that can be constructed from the quarks-are all physics which 
are therefore invisible in the background geometry. We will learn nothing new from the 
background; our study is of the response of the probe D7-branes to the background, 
and this is where the new physics emerges from. 

We carefully study the physics of the probe itself as it moves in the background 
geometry. The coordinates of the probe in the background are fields in an effective D7- 
brane world-volume theory, and the geometry of the background enters as couplings 
controlling the dynamics of those fields. One such coupling in the effective model 
represents the local separation, L{u), of the D7-brane probe from the D3-branes, 
where u is the radial AdSs-Schwarzschild coordinate. 

In fact, the asymptotic value of the separation between the D3-branes and D7- 
brane for large u yields the bare quark mass niq and the condensate vacuum expectation 
value (vev) {'tpip) as follows [60, 76]: 

lim L('u) = m+ + ... , (169) 

u— >oo U 
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where m = 2TTa'mq and — c = {^il)) / [Sir'^ a' N fTj) ^ in this chapter we wiU set Nj = 1. 
The fundamental string tension is defined as T = l/(27ra') here, and the D7-brane 
tension is ry = (27r)~^(a')~^. The zero temperature behavior of the D7-branes in 
the geometry is simple. The D7-brane world-volume actually vanishes at finite u, 
corresponding to the part of the brane wrapped on the shrinking to zero size. 
The location in u where this vanishing happens encodes the mass of the quark, or 
equivalently, the separation of the probe from the D3-branes. In addition, in the zero 
temperature background, the only value of c allowed is zero, meaning no condensate 
is allowed to form, as is expected from supersymmetry. 

The finite temperature physics introduces an important new feature. As is standard 
[38] , finite temperature is studied by Euclideanizing the geometry and identifying the 
temperature with the period of the time coordinate. The horizon of the background 
geometry is the place where that shrinks to zero size. The D7-brane is also wrapped 
on this S^, so it can vanish at the horizon, if it has not vanished due to the shrinking 
of the S^. For large quark mass compared to the temperature (horizon size), the 
shrinking will occur at some finite u> u^., and the physics will be similar to the zero 
temperature situation. However, for small quark mass, the world-volume will vanish 
due to the shrinking of the corresponding to the D7-branes going into the horizon. 
This is new physics of the flavor sector. 

The authors of ref. [12] explored some of the physics of this situation (the depen- 
dence of the condensate and of the meson mass on the bare quark mass) , and predicted 
that a phase transition should occur when the topology of the probe D7-brane changes. 
However, they were not able to explicitly see this transition because of poor data res- 
olution in the transition region, coming from using UV boundary conditions on the 
scalar fields on the D7-brane world-volume. The origin of this phase transition, as 
we shall see, is as follows: The generic behavior of an allowed solution for L{u) as 
in equation (169), is not enough to determine whether the behavior corresponds to 
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an S'^-vanishing D7-branc or an S'^-vanishing D7-branc. The choices of branch of 
solutions have different values of c, generically. In other words, for a given value of 
m there can be more than one value of c. There are therefore two or more candidate 
solutions potentially controlling the physics. The actual physical solution is the one 
which has the lowest value for the D7-brane's free energy. The key point is that, at a 
certain value of the mass, the lowest energy solution may suddenly come from a differ- 
ent branch, and, as the corresponding value of the condensate changes discontinuously 
in moving between branches, we find that the system therefore undergoes a first order 
phase transition. On the gauge theory side, we can imagine a similar situation occur- 
ring; two different branches of solution are competing, and the lowest energy branch 
is always picked. We are able to uncover this physics by doing a careful numerical 
analysis of the equations of motion for the probe dynamics on the gravity side of the 
AdS/CFT correspondence, by using IR boundary conditions instead of UV boundary 
conditions. 

3.2 Holographic meson melting 
3.2.1 General set up 

We begin by reviewing the physics of the D7-brane probe in the AdSs-Schwarzschild 
background solution[12]. The metric is given by: 



ds'^ = --^dt^ + ^^du^ + —,dx-dx + R^dO'^ + cos^ Od^l + R^ sit? edcj)'^ , (170) 



R^ 



,2 



where i; is a three vector. 



/(«) 




2 ' 
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and the quantity B? 



is given by: 



i?2 = ./^lujsNca' , 



where Qs is the string coupHng (which, with the inverse string tension a' sets for 
example, Newton's constant). The quantity h is related to the mass of the black-hole, 
6^ = 8G5mb.h./(37r). The temperature of the black hole can be extracted using the 
standard Euclidean continuation and requiring regularity at the horizon. Doing this 
in the metric given by equation (170), we find that f3~^ = hj-KB?. Therefore, by 
picking the value of 6, wc arc choosing at what temperature we are holding the theory. 
We choose to embed the D7-brane probe transverse to 9 and (f). In order to study 
the embeddings with the lowest value of the on-shell action (and hence the lowest 
free energy), we choose an ansatz of the form = and 6 = 0{u). The asymptotic 
separation of the D3- and D7-branes is given by L{u) = usmO. Given this particular 
choice of embedding, the world-volume of the D7-brane is given by: 



where g is the determinant of the induced metric on the D7-brane given by the pull- 
back of the space-time metric We are interested in two particular cases. First, 
there is the case where u goes to 6, which corresponds to the D7-brane probe falling into 
the event horizon. In the Euclidean section, this case corresponds to the shrinking of 
the of periodic time, we will refer to these solutions as the 'black hole' embeddings. 
Second, there is the case where goes to 7r/2, which corresponds to the shrinking of 
the S^. We will name these Minkowski embeddings. It is this change in topology {S^ 




(171) 



82 



versus shrinking) between the different solutions that will correspond to a phase 
transition. The classical equation of motion for 0{u) is: 

When u goes to infinity and the background metric becomes asymptotically AdS5 x , 
the equation of motion reduces to: 

^{v:'e'(u)) + zu^e(u) = 0, (173) 

du 

which has solution: 

0iu) = -fm+4) • (174) 

These two terms are exactly the non-nornializable and normalizable terms correspond- 
ing to a dimension 3 operator (■0V) dual field theory with source m and vacuum 
expectation value (vev) c. 

Next we solve equation (172) numerically using a shooting technique. Shooting 
from infinity towards the horizon, physical solutions are those that have a finite value 
at the horizon. This will only be accomplished for a particular m and c value from 
equation (169) which would have to be delicately chosen by hand. Therefore, if we 
instead start from the horizon with a finite solution, it will shoot towards the physical 
asymptotic solutions we desire. In order to be able to analyze the phase transition 
between the black hole and Minkowski solutions, we shoot from the horizon for the 
condensate solutions and from = tt/2 for the Minkowski solutions. This technique 
avoids having to correctly choose the boundary conditions at infinity, which is a sen- 
sitive procedure, allowing us to have many more data points to analyze the phase 



j + Zu^ cos^ e sin e^Jv? + («4 - 64) qi2 ^ q (J72) 
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transition. We impose the boundary condition that, at our starting point -the horizon, 
we have: 

^'(^)lsi^o = -^tan^(6) 

^'N 153^0 = °° (175) 

We argue that this is the physical boundary condition to take; the first is simply a 
result of taking the limit of u — > 6 in the equation of motion in equation (172), whereas 
the second is a result of requiring no conifold singularity as the shrinks to zero size 
[51]. 

We solve the equation of motion, equation (172), numerically using 6=1 and 
R= 1. These numerical choices correspond to fixing the temperature and to measuring 
lengths in units of the radius of the AdS space; the latter condition also means we are 
choosing a particular relationship between the t'Hooft coupling and the dual quantities: 

A = 5ymA^c=^, r' = l- (176) 

Several D7-branc embedding solutions are shown in Figure 20. The red (solid) lines 
correspond to Minkowski solutions, and the blue (dashed) lines correspond to black 
hole solutions. Prom each of these solutions, we can extrapolate the bare quark mass 
and quark condensate vev. 

3.2.2 The first order phase transition 

We plot the c values as a function of m in Figure 21(a). When enlarged, as shown in 
Figure 21(b), we find, as anticipated in the introductory remarks, the multi-valuedness 
in c for a given m. Physics will choose just one answer for c. There is therefore the 
possibility of a transition from one branch to another as one changes m. 



84 



L = 



u sin e 



1.75 



1.5 



1.25 



0.25 



0.75 



0.5 



1 




p = u cos e 



1 



2 



3 



4 



5 



Figure 20: Solutions for the D7-brane probe in the black hole background. 

In order to determine exactly where the transition takes place, we have to calculate 
the free energy of the D7-brane. In the semi-classical limit that we are considering, 
the free energy is given by the on-shell action times (3^^. For our case, this is simply 
given by [60]: 



where here Nf = 1. We calculate this integral numerically using our solutions for 9{u), 
and we plot the results for the energy in Figure 21(c) after we regulate the result by 
subtracting off the energy from the ^ = solution. There is again multivaluedness at 
the same value of m as before, and we zoom in on this neighborhood in Figure 21(d). 
If one follows the solutions of lowest energy for a given m, one can clearly see that 
there is a crossover from one branch to another as m changes. This was also observed 
in ref. [35]. Therefore, we find a first order phase transition — at m ~ 0.92345 — where 
the condensate's vev jumps discontinuously. This had been deduced on other grounds 
in ref. [8] . We show where the jump between curves occurs with the dashed green line 




(177) 
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Figure 21: Probe free energy (in units of 2t: t-jN jh/ E?) and condensate vev at the 
phase transition. See text for more details. 



in figures 21(b) and 21(d). We note that the shaded areas in Figure 21(b) are equal; 
this is to be expected since c and m are thermodynamically conjugate to each other, 
and the area of the graph has the interpretation as a free energy difference d£ ~ — c dm. 
This first order phase transition was studied in details in ref. [48], where the authors 
developed an appropriate formalism to study the meson spectrum of the theory before 
and after the phase transition. They were able to show that after the phase transition 
there are no bound meson states, but only quasi-normal excitations corresponding 
to a melting mesons. This is why the present understanding of the observed first 
order phase transition is that it corresponds to the confinement /deconfinement phase 
transition of the fundamental matter which in the case of consideration is the meson 
melting phase transition. 
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3.3 Concluding remarks 

We expect that, while the details of this construction will not persist in a "realistic" 
QCD string dual, the phase transition itself represents strongly coupled dynamics that 
may well persist as part of the full story of the QCD phase diagram. The transition is 
exciting in itself, of course (particularly since its dual involves a change of topology in 
the D7-branc world- volume) , but much further work is needed on several questions. 
For example, the robustness of the phase transition against corrections would be 
interesting to study. In fact it is believed that the observed first order phase transition 
is an artifact of the large Nc limit. We will come back to this issue in Chapter 7, 
when we study the critical behavior of the theory. Once again a peculiar feature of the 
theory will turn out to be an artifact of the large Nc limit. 
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Chapter 4: Phase structure of finite temperature large 
Nc flavored Yang— Mills theory in an external magnetic 
field 



4.1 The string background 

Consider the AdSs-SchwarzschildxS'^ solution that we studied in the previous section 
given by: 

where dUl = dO"^ + cos^ Od^l + siv? Odcp'^ , 
and dnl = dtjj'^ + cos^ i/jdp + sin^ edj"^ . 

The dual gauge theory will inherit the time and space coordinates t = and x = 
{x^,x'^,x^) respectively. Also, in the solution above, u G [0, oo) is a radial coordinate 
on the asymptotically AdSs geometry and we are using standard polar coordinates 
on the S^. The scale R determines the gauge theory 't Hooft coupling according to 
= oc' \J dYM^c- ^'^^ purpose of our study it will be convenient [12] to perform 
the following change of variables: 

with p = r cos 9 , L = r sin 6 . 
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The expression for the metric now takes the form: 

Following ref . [50] , we introduce fundamental matter into the gauge theory by placing 
D7-brane probes into the dual supergravity background. The probe brane is parame- 
terized by the coordinates {xq, xi,X2, X3, p, ip, f3, 7} with the ansatz from equation (23) 
that we replicate here: 

^ = const, L = L{p) . 

In order to introduce an external magnetic field, we excite a pure gauge i?-field along 
the {x'^,x^) directions as in Chapter 2: 

B = Hdx^ A dx^ (180) 

where H is a real constant. As explained in Chapter 2, while this docs not change the 
supergravity background, it has a non-trivial effect on the physics of the probe, which 
is our focus. To study the effects on the probe, let us consider the general (Abelian) 
DBI action: 

Sdbi = -NfTD7 J d^^ det'/\P[Gab + Bab] + 2WF„b) , (181) 
Ms 

where Tm = Pt/Qs = [(2vr)^a'^(7s]~^ is the D7-brane tension, P[Gab] and P[Bab] are 
the induced metric and induced S-field on the D7-branes' world-volume, F^b is the 
world-volume gauge field, and Nf = 1 here. It was shown in Chapter 2 that, for the 
AdS5 X geometry, we can consistently set the gauge field Fab to zero to leading 
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order in a', and the same argument applies to the finite temperature case considered 
here. The resulting Lagrangian is: 

(54 + 4(p2 + L(p)2)2) J 

(182) 

For large p 3> 6, the Lagrangian asymptotes to: 

C ^ -pVi + ^'(p)' , (183) 

which suggests the following asymptotic behavior for the embedding function L{p): 

L{p) = m+ ^ + ... , (184) 

where the parameters m (the asymptotic separation of the D3- and D7-branes) and 
c (the degree of transverse bending of the D7-brane in the {p, <j)) plane) are related to 
the bare quark mass ruq = m/2Tra' and the quark condensate (i'lp) oc — c respectively 
[60] (this calculation was repeated in the Appendix of Chapter 2). It was shown in 
Chapter 2 that the presence of the external magnetic field spontaneously breaks the 
chiral symmetry of the dual gauge theory (it generates a non-zero {'ipip) at zero m). 
However [12] , the effect of the finite temperature is to melt the mesons and restore the 
chiral symmetry at zero bare quark mass. Therefore, we have two competing processes 
depending on the magnitudes of the magnetic field H and the temperature T = h/irF?. 
This suggests an interesting two dimensional phase diagram for the system which we 
shall study in detail later. 
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To proceed, it is convenient to define the following dimensionless parameters: 



P 

m 



p 

b ' 

L{bp) 



m 



m = 



= m + ^ + . . . . 

p2 



(185) 



This leads to the Lagrangian: 



V 



(p2 + L(p)2) ) 



16 



16 ( p2 + 



L{pr)\' 



1 + 



4(p2 + L(p)2)' 



(186) 



For small values of rj, the analysis of the second order, non-linear differential equation 
for L(p) derived from equation (186) follows closely that performed in refs. [6, 12, 
67]. The solutions split into two classes: the first class are solutions corresponding to 
embeddings that wrap a shrinking in the part of the geometry and (when the 
vanishes) closes at some finite radial distance r above the black hole's horizon which 
is located at r = b/\/^. These embeddings are referred to as 'Minkowski' embeddings. 
The second class of solutions correspond to embeddings falling into the black hole, 
since the of the Euclidean section, on which the D7-branes are wrapped, shrinks 
away there. These embeddings are referred to as 'black hole' embeddings. There is 
also a critical embedding separating the two classes of solutions which has a conical 
singularity at the horizon, where the wrapped by the D7-brane shrinks to zero 
size, along with the S^. If one calculates the free energy of the embeddings, one can 
show [6, 12, 67] that it is a multi- valued function of the asymptotic separation m, 
which amounts to a first order phase transition of the system (giving a jump in the 
condensate) for some critical bare quark mass rricr- (For fixed mass, we may instead 
consider this to be a critical temperature.) We show in this chapter that the effect of 
the magnetic field is to decrease this critical mass, and, at some critical magnitude of 
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the parameter rjcr, the critical mass drops to zero. For r] > ?]„■ the phase transition 
disappears, and only the Minkowski embeddings are stable states in the dual gauge 
theory, possessing a discrete spectrum of states corresponding to quarks and anti- 
quarks bound into mesons. Furthermore, at zero bare quark mass, we have a non-zero 
condensate and the chiral symmetry is spontaneously broken. 

4.2 Properties of the solution 
4.2.1 Exact results at large mass 

It is instructive to first study the properties of the solution for m ^ 1. This approx- 
imation holds for finite temperature, weak magnetic field, and large bare quark mass 
m, or, equivalently, finite bare quark mass m, low temperature, and weak magnetic 
field. 

In order to analyze the case m S> 1, let us write L(p) = a + ({p) for a S> 1 and 
linearize the equation of motion derived from equation (186), while leaving only the 
first two leading terms in (p^ + m^)~^. The result is: 

9p(p'0 - 2^^-2.3 ^ + ^oT^rrt^ ^ + o{0 = ■ (i87) 

(m^ + p^)"^ 2{m^ + p^)^ 
Ignoring the 0{Q terms in equation (187), the general solution takes the form: 

^2 2(772 -I- 1)2 _ 2 

4/9"=(m^ + p'^) 96/9^(771"^ + p'^)-^ 

where we have taken ^'(0) = C(0) = 0. By studying the asymptotic behavior of this 
solution, we can extract the following: 

rf l + Arf + 2r]'^ ^ / 1 
= "-453 + ^2^^ + 

c = l^-i^^^ + ofl) . (189) 
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By inverting the expression for m, we can express c in terms of fh: 



Finally, after going back to dimensionful parameters, we can see that the theory has 
developed a fermionic condensate: 

m)^-c=--H''^ 96iS^ ■ (1^^ 

The results of the above analysis can be trusted only for finite bare quark mass and suf- 
ficiently low temperature and weak magnetic field. As can be expected, the physically 
interesting properties of the system should be described by the full non-linear equation 
of motion of the D7-brane. To explore these we need to use numerical techniques. 

4.2.2 Numerical analysis 

We solve the differential equation derived from equation (186) numerically using Math- 
ematica. It is convenient to use infrared initial conditions [3, 6]. For the Minkowski 
embeddings, based on symmetry arguments, the appropriate initial conditions are: 

L(p)|p=o = -^>in, L'(p)|p=o = 0. (192) 

For the black hole embeddings, the following initial conditions: 

L 



-^(P)le.h. = An, -L'(p)|e.h. ^ 



(193) 
e.h. 



ensure regularity of the solution at the event horizon. After solving numerically for 
L(p) for fixed value of the parameter 77, we expand the solution at some numerically 
large Pmax, and, using equation (184), we generate the plot of — c vs m. It is instructive 
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to begin our analysis by revisiting the case with no magnetic field (r/ = 0), familiar 
from refs. [6, 12, 67]. The corresponding plot for this case is presented in Figure 22. 
Also in the figure is a plot of the large mass analytic result of equation (190), shown 




Figure 22: The solid curve starting far left (red) represents solutions falling into the 
black hole, the dotted (blue) curve represents solutions with shrinking S"^. The vertical 
dashed line corresponds to the critical value of rh at which the first order phase tran- 
sition takes place. The solid black curve dropping sharply from above is the function 
derived in equation (190), corresponding to the large mass limit. 

as the thin black curve in the figure, descending sharply downwards from above; it can 
be seen that it is indeed a good approximation for m > rhcr- Before we proceed with 
the more general case of non-zero magnetic field, we review the techniques employed 
in ref. [6] to determine the critical value of rh. In Figure 23, we have presented the 
region of the phase transition considerably magnified. Near the critical value rhcr, the 
condensate c is a multi- valued function of rh, and we have three competing phases. 
The parameter c is known [60] to be proportional to the first derivative of the free 
energy of the D7~brane, and therefore the area below the curve of the — c vs rh plot is 
proportional to the free energy of the brane. Thus, the phase transition happens where 
the two shaded regions in Figure 23 have equal areas; furthermore, for rh < rhcr, the 
upper-most branch of the curve corresponds to the stable phase, and the lower-most 
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Figure 23: The area below the (— c, m) curve has the interpretation of the free energy 
of the D7-brane; thus the phase transition pattern obeys the "equal area law" — the 
area of the shaded regions is equal. 

branch of the curve corresponds to a meta-stable phase. For rh > rhcr, the lower-most 
branch of the curve corresponds to the stable phase, and the upper-most branch of 
the curve corresponds to a metastable phase. At m = rricr we have a first order phase 
transition. It should be noted that the intermediate branch of the curve corresponds 
to an unstable phase. 

Now, let us turn on a weak magnetic field. As one can see from Figure 24, the effect 
of the magnetic field is to decrease the magnitude of rhcr- In addition, the condensate 
now becomes negative for sufficiently large m and approaches zero from below as m — > 
00. It is also interesting that equation (190) is still a good approximation for m > rhcr- 
For sufficiently strong magnetic field, the condensate has only negative values and the 
critical value of rh continues to decrease, as is presented in Figure 25. If we further 
increase the magnitude of the magnetic field, some states start having negative values of 
m, as shown in Figure 26. The negative values of m do not mean that we have negative 
bare quark masses; rather, it implies that the D7-brane embeddings have crossed L = 
at least once. It was argued in ref. [12] that such embeddings are not consistent with a 
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Figure 24: The effect of the weak magnetic field is to decrease the values of rhcr and 
the condensate. Equation (190) is still a good approximation for rh > rhcr. 

holographic gauge theory interpretation and are therefore to be considered unphysical. 
We will adopt this interpretation here, therefore taking as physical only the m > 
branch of the —cvsm plots. However, the prescription for determining the value of 
rricr continues to be valid, as long as the obtained value of rhcr is positive. Therefore, 
we will continue to use it in order to determine the value of = r/cr for which rhcr = 0. 
As one can see in Figure 27, the value of r/cr that we obtain is rycr ~ 7.89. Note also 
that, for this value of rj, the Minkowski m = embedding has a non-zero fermionic 
condensate Ccr, and hence the chiral symmetry is spontaneously broken. For r] > 
rjcr, the stable solutions are purely Minkowski embeddings, and the first order phase 
transition disappears; therefore, we have only one class of solutions (the blue curve) 
that exhibit spontaneous chiral symmetry breaking at zero bare quark mass. Some 
black hole embeddings remain meta-stable, but eventually all black hole embeddings 
become unstable for large enough ij. This is confirmed by our study of the meson 
spectrum which we present in later sections of the chapter. The above results can be 
summarized in a single two dimensional phase diagram which we present in Figure 28. 
The curve separates the two phases corresponding to a discrete meson spectrum (light 
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Figure 25: For strong magnetic field the condensate is negative. The value of rhcr 
continues to drop as we increase r]. 

mesons) and a continuous meson spectrum (melted mesons) respectively. The crossing 
of the curve is associated with the first order phase transition corresponding to the 
melting of the mesons. If we cross the curve along the vertical axis, we have the phase 
transition described in refs. [6, 12, 67]. Crossing the curve along the horizontal axis 
corresponds to a transition from unbroken to spontaneously broken chiral symmetry, 
meaning the parameter c jumps from zero to Ccr ~ 4.60, resulting in non-zero quark 
condensate of the ground state. It is interesting to explore the dependence of the 
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Figure 26: For sufficiently liigli values of rj there are states with negative m which are 
considered non-physical. However, the equal area law is still valid as long as rricr > 0. 

fermionic condensate at zero bare quark mass on the magnetic field. From dimensional 
analysis it follows that: 

ccr = b'cU^) = ^-^R'H^" ■ (194) 

In the T ^ limit, we should recover the result from Chapter 2: Ccr ^ 0.226R^H^/'^, 
which implies that Ccr{rj) ~ 0.226r/^/^ for ?] S> 1. The plot of the numerically extracted 
dependence Ccr(f?) is presented in Figure 29; for rj > r]cr, Ccr(f?) very fast approaches 
the curve 0.2267/'^/^. This suggests that the value of the chiral symmetry breaking 
parameter Cqt depends mainly on the magnitude of the magnetic field H, and only 
weakly on the temperature T. 

4.3 Thermodynamics 

Having understood the phase structure of the system, we now turn to the extraction 
of various of its important thermodynamic quantities. 
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Figure 27: For rj = r/cr the critical parameter mcr vanishes. There are two m = 
states with equal energies, one of them has non-vanishing condensate — Ccr ~ —4.60 
and therefore spontaneously breaks the chiral symmetry. 

4.3.1 The free energy 

Looking at our system from a thermodynamic point of view, we must specify the 
potential characterizing our ensemble. We are fixing the temperature and the magnetic 
field, and hence the appropriate thermodynamic potential density is: 



where // is the magnetization density and S is the entropy density of the system. 
Following ref. [68], we relate the on-shell D7-brane action to the potential density F 
via: 



dF 



SdT - fidH 



(195) 



F = 2TT^NfTD7lD7 , 



(196) 
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Figure 28: The curve separates the two phases corresponding to discrete meson spec- 
trum (light mesons) and continuous meson spectrum (melted mesons). 

where (here, A'^^ = 1): 

pmax 2 4 ~ 

Id, = h' j d~pf (l - (l + (i?^) ' V'^n^ + 4ound; (197) 

Pmin 

^2 

7? = -^F; f = r/6; ~p = /j/ft; L = L/6; = + ^2_ 

In principle, on the right hand side of equation (196), there should be terms propor- 
tional to which subtract the energy of the magnetic field alone; however, as we 
comment below, the regularization of //jy is determined up to a boundary term of the 
form const x . Therefore, we can omit this term in the definition of F. The boundary 
action /bound contains counterterms designed [13] to cancel the divergent terms coming 
from the integral in equation (197) in the limit of Pmax cxd. A crucial observation is 
that the finite temperature does not introduce new divergences, and we have the usual 
quartic divergence from the spatial volume of the asymptotically AdSs spacetime [53] . 
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Figure 29: The solid curves are the numerically extracted dependence Ccr(??), while 
the dashed curve represents the expected large r] behavior Ccr(??) ~ 0.226ry^/^. The 
solid curve segments at the bottom left and to the upper right (blue) are the stable 
states. The straight segment and the arc that joins it (lower right, red) and red are 
the unstable states. The rest (cyan) are meta-stable states. 

The presence of the non-zero external magnetic field introduces a new logarithmic 
divergence which can be cancelled by introducing the following counterterm: 



-^log (^) / d'xV^^B.^B^'^ , (198) 
where 7 is the metric of the 4-dimensional surface at p = pmax- Note that in our case: 



^V^B^^B^'' = , (199) 



which gives us the freedom to add finite terms of the form const x at no cost to the 
regularized action. This makes the computation of some physical quantities scheme 
dependent. We will discuss this further in subsequent sections. The final form of /bound 
in equation (197) is: 



/bound = -\pL. - \r'h' log ^ . (200) 
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It is instructive to evaluate the integral in equation (197) for the L = embedding at 
zero temperature. Going back to dimensionful coordinates we obtain: 



Pmax I 



^ - |/5max+2-'^ ^ 1 §" 



rPmax+^i^'^'log ^ + ^(l+log4-log i72)+0(p-L) . 

(201) 

The first two terms are removed by the counter terms from /bound) and we are left 
with: 

F{b = 0,m = 0,H) = 27r^NfTD7^^{l + log4:-logH^) . (202) 

8 

This result can be used to evaluate the magnetization density of the Yang-Mills plasma 
at zero temperature and zero bare quark mass. Let us proceed by writing down a more 
general expression for the free energy of the system. After adding the regulating terms 
from /bound) we obtain that our free energy is a function of m, b, H: 

F{b, m, H) = 27r^ NfTmb'^lDrirh, rj^) + F(0, 0, H) , (203) 

where /£)7(m, rj) is defined via: 



Pmax / ^ 1 \ 

in. = / .p(^=(l-J^)(l + ^i55^)^ym^-^j-A,./4 

Pmin 

-^r?2logp^ax - ^?7'(1 + log 4 - log 7/2); f^ = p^ + L{pf . (204) 



In order to verify the consistency of our analysis with our numerical results, we derive 
an analytic expression for the free energy that is valid for m » y/fj. To do this we use 

that for large m, the condensate c is given by equation (190) which we repeat here: 
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as well as the relation dlDi/dm = —2c. We then have: 

m 

/I 1 _|_ 4^2 _|_ Q„4 

c(m, r7)dm + e(r?) = e(^)--r?^ log m- ^^^^^^^ ^ +0(l/m^) , (206) 

where the function ^(ry) can be obtained by evaluating the expression for ID^ from 
equation (204) in the approximation L ^ rh. Note that this suggests ignoring the term 
L'^ which is of order c^. Since the leading behavior of at large m is this means 

that the results obtained by setting L'^ = can be trusted to the order of 1/m, and 
therefore we can deduce the function £,{f]), corresponding to the zeroth order term. 
Another observation from earlier in this chapter is that the leading behavior of the 
condensate is dominated by the magnetic field and therefore we can further simplify 
equation (204): 



Pmax j 

Id7 ^ . lim / dpp^{\ l+ -l)--??^logpniax- -77^(1 -log—) 

' 

2 2 2 

= -|-iog^-|-(3-log^) + 0(l/m=^) . (207) 



Comparing to equation (206), we obtain: 



e(^) = -^(3-log^), (208) 



and our final expression for //jy, valid for rh S> ^Jfi^. 

in. = 4 (3 - log f ) - i,^lc^™ - iijl^ + OiUnf ) . (209) 
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4.3.2 The entropy 

Our next goal is to calculate the entropy density of the system. Using our expressions 
for the free energy we can write: 



dT J jj db \ / 

= -2Tr^ R^NfTorb^ I 4/d7 + 2cm - 4^^??^ j = 2iT^B?NfTD7b^S{m, rj^) . (210) 



It is useful to calculate the entropy density at zero bare quark mass and zero quark 
condensate. To do this, we need to calculate the free energy density by evaluating the 
integral in equation (204) for L = 0. The expression that we get for /d7(0, rf') is: 

^7(0, n') = \{^- - log ^ "^'^' ^ ■ (211) 

The corresponding expression for the entropy density is: 

S\m=o = 2n^R'NfTD7T^ ("2 + + Sl^) " ^^^^^ 

One can see that the entropy density is positive and goes to zero as T — 0. Our next 
goal is to solve for the entropy density at finite m for fixed 77. To do so, we have to 
integrate numerically equation (210) and generate a plot of S versus m. However, for 
rh ^ y/rj we can derive an analytic expression for the entropy. After substituting the 
expression from equation (209) for Id7 into equation (210) we obtain: 



or if we go back to dimensionful parameters: 

Sib, m, H) = 2n'R^NfTnrb' (^^^^^^) • (214) 
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One can see that if we send T ^ 0, while keeping r] fixed we get the T'' behavior 
described in ref. [68], and therefore the (approximate; Nf/N <^ 1) conformal behavior 
is restored in this limit. In Figm'e 30, we present a plot of S versus rh for rj = 0.4 . 
The solid smooth black cm've corresponds to equation (213). For this S is positive and 
always a decreasing function of m. Hence, the entropy density at fixed bare quark mass 
m = mb, given hy S = R"^ N fT^jm^ S / , is also a decreasing function of m and 
therefore an increasing function of the temperature, except near the phase transition 
(the previously described crossover from black hole to Minkowski embeddings) where 
an unstable phase appears that is characterized by a negative heat capacity. 




Figure 30: A plot of S versus rh for rj = 0.4. The thin (sharply descending and extend- 
ing to the right) black curve corresponds to that large mass result of equation (213). 



4.3.3 The magnetization 

Let us consider equation (202) for the free energy density at zero temperature and zero 
bare quark mass. The corresponding magnetization density is given by: 

^) = R^NfTm -log-. (215) 




105 



Note that this result is scheme dependent because of the freedom to add terms of the 
form const x H"^ to the boundary action that we discussed earlier. However, the value 
of the relative magnetization is given by: 

M-Mo = -(^|^) -Ho = -'^T^''R^NfTD7l?{^\ = 2Tr^ NfTorb^ , (216) 

is scheme independent and is the quantity of interest in the section. In equation (216), 
we have defined fi = —dlD7/dr)\rh as a dimensionless parameter characterizing the 
relative magnetization. The expression for jl follows directly from equation (204): 

Pmax 

/p (Ar — 1) v V 

^^^r77f^^W=T^ + ^ P^-- - ^ log ^. (217) 
r^y^(4r4 + 1)^ + 16r7r^ 2 2 

Pmin 

For the large m region we use the asymptotic expression for I^^ from equation (209) 
and obtain the following analytic result for jl: 

A = I - I log I + r?logm + ^^^^^ + 0{l/m') . (218) 

We evaluate the above integral numerically and generate a plot of jl versus rh. A plot 
of the dimensionless relative magnetization fl versus m for 77 = 0.5 is presented in 
Figure 31. The black curve corresponding to equation (218) shows good agreement 
with the asymptotic behavior at large m. It is interesting to verify the equilibrium 
condition djl/dT > 0. Note that since /iq does not depend on the temperature, the 
value of this derivative is a scheme independent quantity. From equations (216) and 
equation (218), one can obtain: 

I = ^r'R'N^To., (2, - 2|,) = 2.VN,To.^, > , (210) 
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which is vahd for large m and weak magnetic field H. Note that the magnetization 
seems to increase with the temperature. Presumably this means that the temperature 
increases the "ionization" of the Yang-Mills plasma of mesons even before the phase 
transition occurs. 




Figure 31: A plot of the dimensionless relative magnetization fi versus m for r] = 0.5. 
The thin black curve (starting with a steep descent) corresponds to the large mass 
result of equation (218). 



4.4 Meson spectrum 

In this section, we calculate the meson spectrum of the gauge theory. The mesons 
we are considering are formed from quark-antiquark pairs, so the relevant objects to 
consider are 7-7 strings. In our supergravity description, these strings are described 
by fluctuations (to second order in a') of the probe branes' action about the classical 
embeddings we found in the previous sections [59]. Studying the meson spectrum 
serves two purposes. First, tachyons in the meson spectrum from fluctuations of the 
classical embeddings indicate the instability of the embedding. Second, a massless 
meson satisfying a Gell-Mann-Oakes-Renner (GMOR) relation [34] will confirm that 
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spontaneous chiral symmetry breaking has occurred. As a reminder, in ref. [59], the 
exact meson spectrum for the AdSs x background was found to be given by: 



M{n,i) = '^^in + e+l){n + e + 2) , (220) 



where £ labels the order of the spherical harmonic expansion, and n is a positive integer 
that represents the order of the mode. The relevant pieces of the action to second order 
in a' are: 



Fro^ A B(o, A P 



"(2) A -D(2) 
iMs 



C(4) 



S/Nf = -Tor J (f^gab + Bab + 2Tra'Fab + (27ra') fij 

+ (27ra') / i^{2) A F(2) A P [C(4)] , (221) 



'Ms 



1 



C(4) = ——dtAdx^Adx^Adx^, (222) 



^4 

(7(4) = ~ (l ~ cos^ 0) sin cos dtp A d(f)2 A d(j)z A d^i , (223) 

9s 

where P [C(4)] is the pull-back of the 4-form potential sourced by the stack of Nc 



D3-branes, P 



C(4) 



is the pull-back of the 4-form magnetic dual to C(4), and F^2) 
is the Maxwell 2-form on the D7-brane world-volume. At this point, we resort to a 
different set of coordinates than we have been using. Instead of using the coordinates 
(p, L) introduced in equation (179), we return to the coordinates {z = 1/n^, 9) because 
the analysis is simpler. We consider fluctuations of the form: 

e = eoiz) + 27ra'xie) , (224) 
ct>i = 27ra'$(r) , (225) 

where the indices a,b = ... 7 run along the world-volume of the D7-brane. 6o{z) 
corresponds to the classical embedding from the classical equations of motion. Plugging 
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the ansatz in equations (224) and (225) into the action and expanding to second order 
in (27ra'), we get as second order terms in the lagrangian: 



—Cp2 

-Cf-x 

L.p2 



t-'P-^ 



'^S'^'R'daxdbX - ^^i?" [d'oT E^'S'^daxdbX 



+ 2^ 



ES^^R^siu^Ooda^db^ , 

a, (^V^rXe^'j^^) - J^^deV^] = xF23f , 

fP rr mnop 



2 

1 

4 

XF23 

1 1 



(226) 



— ^FoiB23R'^ sin tp cos tpdz (l — cos^6'o) = —^FoiB23R'^ sintp cos tpdzK . 



We have taken Eab = g^^b + ^"-^ zeroth order contribution from the DBI 

action. In addition, we use that E"-^ = S"-^ + J"^ where 5"^ = 5''" and J"^ = - J^". 
We use this notation for brevity. The indices m,n, o,p = A . . .9 run in the transverse 
directions to the D3-branes. Prom these lagrangian terms, we derive the equation of 
motion for x to be: 



= da\ V-ES'^'R 



,,^Jl + 46V(0(,)2~ 



dbX 



-F23 



y l + 4z2(0(,)2 ^ 

ER^OqE'^^J^^^ - J^^de 



X 



-E 



dz [E'^R^'^de 



-E 



-E 
(227) 



The equation of motion for $ is given by: 



dai^^f^S^-^R^ sm^ dod^^^ - FQiB23R'^shi^cos'4;dzK = 0. (228) 
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The equation of motion for Af, is given by: 



da (-V^S'^'^' S'^' Fa>b' - Xf {SISI - 6^6^) + B23^d,K (^6^S\ - S^S^) (229) 

We are allowed to set = with the constraint (using that 5"^^ = S^^): 

S''^dmdoAo + S'^dmdiAi + S^^dm{d2A2 + d3As) = (230) 

Therefore, we can consistently take = ^i^i = 0, ^2^2 = —dsA^. With this 
particular choice, we have as equations of motion for the gauge field: 

-^0 (y^S^'^S^^doAi^ + d,KB23do^ - d, [^^^S''S^^^,A^) (231) 

-dfa(V^S^^S^^dnAr) = 0, 
-9o (7=^^50052290^2) + fd^x - Qz {sl^S''S''''d,A^ 

-drn{V^S^^S^^dnA2) = 0, 

-do (V^S'^'^S^^doAs) - fd2X - (^V^S^^S^'d^A^) 

-drnly^S^^S^^dfiAs) = 0, 

where the indices rh, h run over the that the D7-brane wraps. If we assume that 
diX = 0, we find that the equations for A2 and A3 decouple from x- Therefore, we 
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can consistently take i<23 = 0, or, in other words, A2 = = 0. This simphfies the 
equations of motion that we need to consider to: 



= da 



'-ES'^'R 



1 + 4^2 ^Ql^ 



2 /'flM2 



dhX 



X 



(232) 



= 



-drniy^S^'^S^^dnAi) , (233) 

da (y^S'^B? sin^ ^o5&^) - -^01^23^?^ sin ^ cos i^d^K . (234) 



In the proceeding sections, we will work out the solutions to these equations numer- 
ically using a shooting method. With an appropriate choice of initial conditions at 
the event horizon, which we explain below, we numerically solve these equations as an 
initial condition problem in Mathematica. Therefore, the D.E. solver routine "shoots" 
towards the boundary of the problem, and we extract the necessary data at the bound- 
ary. 

4.4.1 The X meson spectrum 

In order to solve for the meson spectrum given by equation (232), we consider an 
ansatz for the field x of the form: 



X = h{z) exp {-iut) , 



(235) 



where we are using the same dimensionless coordinates as before, with the addition 
that: 

z = b~'^z , u = R~% Cj . 

In these coordinates, the event horizon is located at 5 = 1. Since there are two different 
types of embeddings, we analyze each case separately. We begin by considering black 
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hole cmbeddings. In order to find the appropriate infrared initial conditions for the 
shooting method we use, we would like to understand the behavior of h{z) near the 
horizon. The equation of motion in the limit of 5 — ^ 1 reduces to: 

h"{z) + j^h'Cz)+ f h{z) = 0. (236) 
z — L 16 (z — 1} 

The equation has solutions of the form (1 — z)'^'^'^/^, exactly of the form of quasinormal 
modes [82]. Since the appropriate fluctuation modes are in-falling modes [48], we 
require only the solution of the form (1 — 5)""*'^/^. This is our initial condition at the 
event horizon for our shooting method. In order to achieve this, we redefine our fields 
as follows: 

hiz) = yiz)il-zr^/\ 

which then provides us with the following initial condition: 

y{z^l) = e, (237) 

where e is chosen to be vanishingly small in our numerical analysis. The boundary 
condition on y'{z — > 1) is determined from requiring the equation of motion to be 
regular at the event horizon. The solution for the fluctuation field y{z) must be com- 
prised of only a normalizable mode, which in turn determines the correct value for lo. 
Since we are dealing with quasinormal modes for the black hole embeddings, uj will 
be complex; the real part of oj corresponds to the mass of the meson before it melts, 
and the imaginary part of a; is the inverse lifetime (to a factor of 2) [48] . We begin by 
considering the trivial embedding 9o{z) = (a black hole embedding). This embedding 
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corresponds to having a zero bare quark mass. The equation of motion (232) simphfies 
tremendously in this case: 



h"(2) + 



2z 



1 



— 1 z{\ + z'^rj^) 



h'(z) + 



3 + 5(-3z + lD2) 

(52 _ 1)2 



hCz) 



. (238) 



We show solutions for lv in Figure 32 as a function of the magnetic field r/. In par- 
ticular, we find the same additional mode discussed in ref. [28]. This mode becomes 
massless and eventually tachyonic at approximately ij ~ 9.24. This point was origi- 
nally presented in Figure 29, where the intermediate unstable phase joins the trivial 
embedding. This is exactly when the — c vs m plot has negative slope for all black hole 
embeddings. 
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Figure 32: The x meson mass as a function of magnetic field for the trivial embedding. 
The upper (red) curve is the generalization of the mode described in ref. [48]. The 
lower (blue) curve is the generalization the mode discussed in ref. [28]. We do not 
extend this second curve to small r] because the numerics become unreliable. 



We now consider embeddings with non-zero bare quark mass. This means solving the 
full equation (232). We have both embeddings to consider; for the black hole embed- 
dings, we will follow the same procedure presented above to solve for the complex ui. 
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We can still use the same procedure because in the limit of 5 1, the equation of 
motion still reduces to equation (236). For the Minkowski embeddings, we do not have 
quasinormal modes, and a) is purely real. Therefore, we use as initial conditions: 

X{z Zmax) = e , (239) 
X'{z Zmax) = OO . (240) 

In Figure 33 and Figure 34, we show solutions for very different magnetic field values. 
In the former case, rj is small, and we do not have chiral symmetry breaking; in the 
latter, r] is large, and we have chiral symmetry breaking. It is important to note that 
in neither of the graphs do we find a massless mode at zero bare quark mass. In Fig- 
ure 33, we find that fluctuations about both the black hole and Minkowski embeddings 
become massless and tachyonic (we do not show this in the graph). The tachyonic 
phase corresponds exactly to the regions in the — c vs rh plot with negative slope. 

Re[w] 




Figure 33: The x meson mass as a function of bare quark mass for tj = 1. The dashed 
(blue) curve corresponds to fluctuations about black hole embeddings. The solid (red) 
line corresponds to fluctuations about Minkowski embeddings. These modes have a 
purely real lo. The straight dashed (black) line corresponds to the pure AdSs x 
solution. 
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Figure 34: The x meson mass as a function of bare quark mass for rj 
(black) line corresponds to the pure AdSs x solution. 



10. The dashed 



4.4.2 The $ and A meson spectra 

Let us now consider the coupled fluctuations of $ and A in equations (233) and (234). 
We consider an ansatz (as before) of the form: 

$ = (p{z) exp {—iujt) , 
Ai = A{z) exp {—iujt) . 

It is interesting to note that, for the trivial embedding 9o{z) = 0, one of the coupled 
equations is equal to zero, and we simply have: 



z 



(2 + ^M3^^-l)),,,,, , 



2 



A"{z) + — i-i- 37r#^'(5) + oA{z) = . 

^ ' (z2 - 1) (1 + z2j^2) V J 4^ (^2 _ 1)2 V J 
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Again, we note that in the hmit of 5 ^ 1, we have: 

A"{z) + ^^A'{z) + — ^A{z) = . 

This is exactly the form of equation (236), so A(z) has the same solutions of in-falling 
and outgoing solutions, which provides us with the necessary initial conditions for our 
shooting method. We show solutions for cj in Figure 35. Unfortunately, we do not 
Re[(b] 




Figure 35: The A meson mass as a function of magnetic field for the trivial embedding. 

know how to solve for the quasinormal modes for other black hole embeddings. Since 
the equations of motion are coupled, we are unable to find an analytic solution for the 
fluctuations near the event horizon. This prevents us from using infrared initial condi- 
tions for our shooting method. However, we are actually more interested in searching 
for the "pion" of our system, which will occur when we have chiral symmetry breaking. 
In those cases, we are only dealing with Minkowski embeddings with a pure real to, so 
we may ignore the black hole embeddings. 
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Since the equations of motion are coupled, it turns out that only for specific initial con- 
ditions will both fluctuations only be comprised of normalizable modes. We represent 
this by a parameter a (which we must tune) as follows: 

^(z^Zmax) = icosa, (241) 
(/)(z ^ Zmax) = sina . (242) 



The initial conditions on A! and 0' are determined from the equations of motion. We 
show several solutions in Figure 36. There are several important points to notice. 
First, we find that the lowest mode satisfies an GMOR relationship given by: 



l.lm^/2 _ (243) 



Therefore, we have found the Goldstone boson of our system related to the breaking 
of chiral symmetry. Second, the modes exhibit the Zeeman splitting behavior that was 
discussed in ref. [29]. It is interesting that the various modes always cross each other 
at approximately m ~ 2. We have no intuitive explanation for this behavior. 




Figure 36: Coupled A — cj) fluctuations for rj = 10. The dashed (black) line corresponds 
to the pure AdSs x solution. 
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Figure 37: Zoom near zero bare quark mass for the lowest mode. 

4.5 Concluding remarks 

We have extended the holographic study of large gauge theory in an external 
magnetic field described in Chapter 2, to the case of finite temperature, allowing us 
to study the properties of the quark dynamics when the theory is in the deconfined 
plasma phase. 

The meson melting phase transition exists only below a critical value of the applied 
field. This is the critical value above which spontaneous chiral symmetry breaking is 
triggered (in the case of zero mass). Above this value, regardless of the quark mass 
(or for fixed quark mass, regardless of the temperature) the system remains in a phase 
with a discrete spectrum of stable masses. Evidently, for these values of the field, it is 
magnetically favorable for the quarks and anti-quarks to bind together, reducing the 
degrees of freedom of the system , as can be seen from our computation of the entropy. 
Meanwhile, the magnetization is greater in this un-melted phase. 

There have been non-perturbative studies of fermionic models in background mag- 
netic field before, and there is a large literature (see e.^., the reviews of refs. [70, 85], 
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and the discussion of ref. [81] and references therein). Generally, those works use quite 
different methods to examine aspects of the physics — some primary non-perturbative 
tools are the Dyson-Schwinger equations in various truncations). Our results (and the 
zero temperature result obtained with these methods in the zero temperature case [29] ) 
are consistent with the general expectations from those works, which is that strong 
magnetic fields are generically expected to be a catalyst for spontaneous chiral sym- 
metry breaking in a wide class of models (see e.g., refs. [43, 81, 70] for a discussion of 
the conjectured universality of this result). 
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Chapter 5: Phase structure of finite temperature large 
Nc flavored Yang— Mills theory in an external electric fleld 



5.1 Introductory remarks 

In this chapter, we study large Nc SU (N) gauge theory with non-backreacting hyper- 
multiplet quark flavors in the presence of a background electric field, at both zero and 
finite temperature. The holographic description of this theory at finite baryon chemical 
potential was first studied in ref. [52] , where the existence of a global electric current 
induced was demonstrated. Further study of this set up was considered in ref. [72]. 
The phase structure of the theory at zero baryon chemical potential was studied in 
ref. [5]. 

We may expect the behavior of the theory to be quite different from the magnetic 
case considered in the previous chapter. In a sense the electric field is pulHng the 
quarks apart, decreasing the bound energy and thus we may expect that the onset 
of the deconfinement phase transition will happen at lower temperature, than in the 
absence of an electric field. Furthermore, at zero temperature we may expect that for 
sufficiently strong electric field the binding energy of the quarks will be completely 
overcome and the mesons will dissociate into their constituent quarks. 

It is worth noting that the phase transition will be driven by the quantum fluctua- 
tions, which will result in a quantum phase transition. On the other side, the mesons 
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are electrically neutral, while the dissociated quarks aren't. Hence the dissociation can 
be seen as a conductor/insulator phase transition [5]. 



5.2 General set up 

In this section we describe the gravitational background and the technique employed 
to introduce an external electric field. We consider the following form for the metric 
of the AdSs-Schwarzschild X 5^ background: 

ds^a' = -'^L_^dt^ + ^dx^ + ^^du'' + R'cos''ednl (244) 

+R^de'^ + sin^ Odcf? , 
e* = Qs', C(4) = ~^dt A dxi A dx2 A dx^ . 

In the above, there are time and space coordinates t = x^ and x = (x^, x^, x^) respec- 
tively (which the dual gauge theory will also have) and also u G [0, oo) which is a radial 
coordinate. Those are the coordinates on the asymptotically AdSs geometry. We are 
using standard polar coordinates on the S^, with: 

dnl = d9^ + cos^ 9dnl + sin^ edcj)"^ , 
and dnl = + cos^ ^d/? + sin^ Odj^ . (245) 

The scale R determines the gauge theory 't Hooft coupling according to R^ = a' ^ Qy^^N . 
The parameter h sets the radius of the horizon of the black hole, which in turn sets the 
temperature, according to T = b/-KR^. Now following the usual technique to introduce 
flavors to the theory we add Nf D7-branes to the theory. The D7-branes are extended 
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along the xq, xi,X2,X3, u, ip, (3, 7 directions of the geometry and have a non-trivial pro- 
file along 6, namely 9 = 9{u). In this chapter we do not study the meson spectrum of 
the theory, thus the relevant part of the lagrangian is: 



-Tb7 J deti/2 (p + P [Bab] + 27ra'Fab) (246) 



In order to introduce electric field we consider a constant S-field along the {xo,xi) 
plane: 

B = Edxo A dxi (247) 

It is convenient to work with the coordinate u. The corresponding lagrangian is given 
by: 

/,,4 _ „,4 j ,,4 _ k4 

jCo,J V^ti' cos3 9\ 1 + 

COS sin (248) 

V — 0* V 

where, 

nt = 6^ + . (249) 

As one can see at u = the action (248) vanish and for tx < u* it becomes imaginary. 
Often we will refer to it* (or equivalently r* ) as the "pseudo- horizon" or the "vanishing 
locus". This implies instability of the D7-brane embedding and is a signal that we 
should generalize the ansatz to include non-trivial U{1) gauge field on the D7-brane. 
As it was shown in ref. [52], one can heal the on-shell action (248) by generalizing the 
ansatz to: 

Ai = Ai{u) . (250) 

using Furthermore, it can be shown that the asymptotic behavior of the gauge field 
at infinity u — > 00 can be related to the global electric current of the theory along xi , 
namely = {'i/j^^i/j). 
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After substituting in the lagrangian (246) we obtain the following lagrangian: 

/,/4 _ ,.4 ,,4 _ U4 „4 _ >,4 

Cocu^cos^eJ . ,1 (1+ . f^si^^cos^ ^ (251) 

V — 0* 

where: 

/ = 2Tra'Ai{u) . (252) 

Now using the fact that f{u) is a cyclic field, we arrive at following the equation of 
motion for f{u): 

which results to: 



tt^COS^ I 



= K = const . (254) 



We will show now that the constant K is related to the global electric current = 
{'tj^^^il)). Indeed, the four dimensional action of the dual gauge theory is obtained after 
integrating along u, VL-^ and regularizing the corresponding integrals. Therefore for the 
variation of the four dimensional action we obtain: 

^54D = j d}xJ^5A^{oo) , (255) 

where we have identified >l^(oo) with the U{1) gauge field in the 4D dual gauge theory. 
On the other side, we have that: 

SS4B = I d^xdnsduSC = f d^xdn^du (^5Ai' + . dJA^ = 

J J \^A^' d{dyArau) J 

= /<i*-*3«.(5(|^M,)=/<iSd!)3^M,(oo), (256) 
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where we have considered variation such that 5^i|horizon = 0. Now after comparing 
equations (255) and (256) we obtain that: 

/dC r BC 

dn. = 27ra' d^^— = iMNfTorK , (257) 

d{duAi) J df 

therefore K is indeed proportional to J^. Now let us proceed by solving for /' from 
equation (254) and substituting in equation (251). The lagrangian of the resulting 
on-shell action is: 

As one can see, the on-shell lagrangian from equation (258) is regular at u = 
provided that we have: 

= ulR'^E'^cos^eo . (259) 

Here is the angle at which the D7-brane enter the vanishing locus (at u = u*). One 
can also see that if the D7-brane close at some Umin > u>t above the vanishing locus 
(Minkowski type of embeddings) , we have to set = to avoid imaginary action. This 
implies that the current is equal to zero. And it better vanish because Minkowski 
type of embeddings describe the confined (meson gas) phase of the fundamental matter 
and hence, there are no free charges in the system. On the other hand, as we arc going 
to demonstrate in the next sections, embeddings which reach the vanishing locus will 
reach the horizon and are thus classified as a black hole embeddings, the latter describe 
deconfined quarks and therefore the theory is in a conductive phase. 

From equation (251), we can derive the equation of motion for 6{u) and substitute 
the expression for /' which is obtained from solving equation (254). Its exact form is 
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not particularly illuminating, this is why wc will not display it fully here 
of large u, the equation of motion asymptotes to: 

4- (u^O'(u)) + 3u^9(u) = , 
du ^ ' 

which has solution: 

e{u) = - + 4 . (261) 

The constants m and c are related to the bare quark mass and the condensate respec- 
tively, in a manner that is by now very standard [6, 12, 50, 59]. In our notation, the 
exact relationship is given by: 

"^1 = TT^, ' (V'V') = -STT^a'NfTmc . (262) 



. In the limit 



(260) 



5.3 Properties of the solutions 
5.3.1 Exact results at large mass 

It is instructive to study the properties of the quark condensate as a function of bare 
quark mass for large mass. This corresponds to m S> i^^-E. Note that this corresponds 
to the so-called Minkowski embeddings for the probe D7-brane for which the constant 
K proportional to the electric current vanishes by virtue of equation (259) since 
Minkowski embeddings correspond to = 7r/2. To extract the behavior of quark mass 
and condensate we linearize the equation of motion obtained from equation (246) in 
the same way as described in the previous chapter on the magnetic case. We find the 
following analytic behavior: 

^^^^ = ^ + 96^;^^ + ^ b • ^'''^ 
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This suggests that for high enough quark mass the condensate vanishes. However, 
unhke the magnetic field case from the previous chapter, the vanishing occurs from 
the positive side of the condensate axis. We will explore these features and more in 
the next section using numerical techniques. 

5.3.2 The case of vanishing temperature 

It is instructive to first study the case of vanishing temperature, since a number of 
important features will already appear in this case. For this simple case, the energy 
scale is set by RVe, and therefore it is convenient to introduce the dimensionless 
quantities u and rh via: 

u = RVEu; m = RVEm. (264) 

We solve the equation of motion for 6{u) using shooting technique and infrared 
boundary condition described in refs. [3, 6, 51]. For Minkowski embeddings we impose 
the following boundary condition (motivated by the avoidance of conical singulari- 
ties [51]): 

e'{u)\0=^/2 = <x . (265) 

For embeddings reaching the pseudo-horizon we find it numerically convenient to shoot 
forward (towards infinity) and backward (towards origin) from the vanishing locus 
= 1. The boundary condition ensuring the smoothness of solutions across the 
vanishing locus is found from the equation of motion itself to be: 



^0 = ^(1); 



89 
du 



■tan 



(266) 
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From the above boundary condition we can see that the embeddings reaching the 
vanishing locus at = t: /2 has a conical singularity. Next we proceed to discuss our 
numerical findings. 

Numerical analysis shows that there are three different classes of embeddings for 
finite electric field and zero background temperature, which can be classified by the 
topology of the probe D7-brane. There are the smooth Minkowski embeddings that 
close before reaching the vanishing locus, due to shrinking of wrapped by the probe 
brane. There are also the embeddings that reach the vanishing locus before the 
shrinks. These can in turn be classified in two different categories. The first are singular 
solutions that reach the vanishing locus but close before reaching the origin. These 
therefore have a conical singularity at the closing point. We will discuss these further 
below. The second are smooth solutions that pass through the vanishing locus and 
reach all the way to the origin with no singular behavior. These different embeddings 
are summarized in figures 38 and 39. 



rsmO 




Figure 38: The solid curves starting far left (red) represent solutions smoothly closing 
before reaching the vanishing locus. The solid remaining curves (blue) correspond 
to embeddings passing the vanishing locus (denoted by a semi-circular dashed black 
curve). 

We can explicitly see the conical singularity appearing for the brane (second curve 
in Figure 39) since it intersects the vertical axis at some non-zero value of 9. De- 
pending on the magnitude of the electric field this class of singular solutions starts 
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Figure 39: Two kinds of embeddings that pass the vanishing locus corresponding to D7- 
brane closing at the origin (left), and forming a conical singularity (right) respectively. 
The semi-circular (red) curve corresponds to the vanishing locus. 

appearing after a certain value of (^o)min' close to 7r/2 (which corresponds to the 
critical embedding), is reached and persists until = ■7r/2. 

The existence of these embeddings suggests the possibility of a transition in topol- 
ogy of the probe brane as a function of the parameters, as happened for finite temper- 
ature. Here, it can happen as a result of the external electric field that we've applied. 
The Minkowski embeddings simply have a shrinking S^, while the embeddings reach- 
ing the origin are distinguished by having an shrinking as well as touching the AdS 
horizon. However, the presence of singular solutions makes the nature of the transition 
subtle, as we discuss later. We will later see that finite temperature case also reveals 
similar classification of embeddings. 

From the asymptotic behavior of these embeddings we can extract condensate as a 
function of the bare quark mass. In Figure 40 we show this dependence. There are two 
important mass scales. One, rhcr, is where the phase transition between the two types 
of embedding occurs, and the other rh* , is where the singular solutions first appear. 
These values are in turn determined by the parameters (^o)cr and (0o)min) respectively. 
If ^Omin > (^o)cr) then the singular solutions are never thermodynamically favored. 
We see the presence of a first order phase transition speculated based on the general 
arguments. Furthermore, we observe that there is no chiral symmetry breaking, since 
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Figure 40: The condensate as a function of bare quark mass rh. The curve segment 
coming in from the left (blue) corresponds to embeddings reaching the vanishing locus 
and the curve segments going out to the right (red) correspond to Minkowski embed- 
dings. There is a family of solutions that contain a conical singularity between the 
vanishing locus and the origin. They exist for (0o)mm < < (^o)cr- On the right is 
a magnification of the turn-around region where these segments join, showing multi- 
valuedness. An analysis of the free energy reveals a first order phase transition at rhcr 
where there is a jump from one type of embedding to another. 

in all cases the condensate vanishes for vanishing mass. This is to be contrasted with 
the case of external magnetic field studied with these methods in the previous chapters. 

This latter observation fits our intuition that the mesons in this theory have a 
binding energy that grows with (it is proportional to) the constituent quark mass [59] . 
For a given quark mass there should exist sufficiently high electric field that can reduce 
the binding energy the quarks. This allows two things to happen: First, this inhibits 
the formation of chiral condensate, and second, this ultimately will dissociate the 
mesons into its constituent quarks. This dissociation is in fact a transition from an 
insulating to a conducting phase, mediated by the external electric field. On the dual 
gravity side this corresponds to a transition from Minkowski embeddings that do not 
reach the origin to embeddings that do. 

We must note that the appearance of the singular solutions (those that have a 
conical singularity) are not well understood at the moment. There is therefore the 
possibility that there is an as yet to be identified intermediate phase right after the 
meson dissociates. We can make our dissociation transition explicit. From equa- 
tion (264) we can deduce the exact dependence of m* and rricr on the electric field E 
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as m* = RsfEm* and m-cr = RsfE rhcr- This latter is the value of the critical mass for 
a given electric field, Using the free energy considerations that we have used in previous 
work (it amounts an equal-area law, see ref. [6]), we can also independently (as a test 
of our methods for later) determine this dependence numerically for our solutions, and 
we display this in Figure 41. We also include the critical mass for which the singular 
solutions appear as a function of the field. We see that the analytic behavior deduced 
above is nicely confirmed by the numerics. Indeed, the singular solutions seem to be 




Figure 41: The phase diagram showing the electric field induced first order phase 
transition separating the dissociated phase from the phase with stable mesons. The 
blue (dashed) curve separates the phases. The region within the blue (solid) curve and 
red (dashed) curve represents the embeddings with conical singularity. 

unavoidable in the phase diagram (at least at the level of our analysis). However, 
they seem to lie very close to the dissociation transition and in the limits — > and 
E ^ oo the two curves tend to merge. 

Presently we do not completely understand the role played by the singular solutions 
in the physics. One simple possibility (but not the only one) is that stringy corrections 
smooth out the singularity in the interior, while preserving the asymptotic behavior. 
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This would mean that we have a simpler phase diagram that that in Figure 41 (i.e., 
with no solid line) . 

We now turn to the finite temperature case. It is worth noting that a simple 
check on our results there will be that the limit of large electric field (compared to 
temperature) should recover the physics that we have seen here. 

5.3.3 The case of finite temperature 

In order to proceed, it is convenient to introduce dimensionless parameters: 

f = r/h; u = u/b; E = R^E/b'^; rh = m/b; c = c/b^; (267) 

we again solve the equation of motion for 6{u) numerically using Mathematica. 
We use similar shooting technique and boundary conditions outlined in the previous 
section. For Minkowski embeddings we impose equation (265); and for for black hole 
embeddings the appropriate initial condition to ensure smoothness of the embedding 
when we stitch the solutions across u* is: 

^0 = 0{u,y, (268) 

de ^ {3ut - 1) - ^/{3ui - ly + mijuj - l) tan^ 00 
du^ 3Uit{uj — 1) tan 6*0 

We show several solutions for £^ = 1 in Figure 42. We can extract from these em- 
beddings the condensate c as a function of the bare quark mass m by analyzing the 
asymptotic behavior of the solutions. We show several of these solutions for various 
dimensionless electric field values in Figure 43. The first order phase transition de- 
scribed in refs. [6, 12] persists in the presence of an electric field. We determine the 
critical mass at which this transition occurs using the free energy analysis (equivalent 
to an equal-area law) described in ref. [6], and we present the phase diagram in Fig- 
ure 44. The numerical results confirm our intuition: at zero temperature, the mass 
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Figure 42: Several solutions for classical D7-brane embeddings with E = 1. 

of the meson is proportional to the bare quark mass, and hence the binding energy 
is proportional to the bare quark mass [59]. Turning on an electric field will decrease 
the binding energy (the electric field will pull the quark-antiquark pairs apart), and 
only mesons built out of sufficiently heavy quarks-the origin of the critical mass-will 
survive. At finite temperature, we can think of the increasing of the critical mass as 
equivalent to the decreasing of the critical temperature. Therefore, lowering the bind- 
ing energy causes mesons to melt at lower critical temperatures. In addition, since 
c = c/b^, the magnitude of the dimensionless condensate increases with decreasing 
critical temperature (as we see in Figure 43). Furthermore, the stronger the electric 
field, the more initial binding energy is required and therefore a higher critical mass is 
needed. 

We can deduce the asymptotic behavior of the phase diagram by considering large 
electric fields; the dominant energy scale is set by the electric field, and therefore, we 
are effectively at zero temperature. To see this, one should note that the energy scale 
in the geometry is set by n=K, and for strong electric field, n* is much larger than the 
position of the event horizon horizon. Therefore, we should expect to reproduce the 
results for pure AdSs x in an external electric field. In the case of pure AdSs x 
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Figure 43: The condensate as a function of bare quark mass for increasing electric 
field. From bottom to top, the electric field values are E = 0.01, 0.1, 0.3, 0.5, 0.7, 0.9. 



in an external electric field, the only energy scale is given by R^/E, and therefore we 
can easily show that the bare quark mass and condensate scale as: 

m oc R^/E , 
c oc R^E^/^ . (269) 

This result suggests that: 

1/mcrit oc 1/71. (270) 

This dependence is checked numerically (and shown in Figure 44) and appears to be 
vahd for E > 2.f 

It is instructive to compare the phase transition to the one for zero electric field 
and finite temperature. In the latter case the Minkowski and black hole embeddings 
are separated by a "critical" embedding with a conical singularity at the event horizon 
as shown in Figure 56. The straightforward generalization of this picture to the case 
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Melted Mesons 




Figure 44: The phase diagram depicting the first order phase transition, denoted by a 
heavy hne (blue) that separates the melted and stable meson phase. The thin black 
curve is the best fit curve given by const, x for large electric field. 



of finite electric field suggests the existence of a "critical" embedding with a conical 
singularity at the vanishing locus {u = u^,). However, similarly to the zero temperature 
case a more detailed study of the transition from Minkowski to black hole embeddings 
reveals an interesting third class of embeddings. These are embeddings which enter 
the vanishing locus at n*, and have a conical singularity above the event horizon. See 
Figure 57. 

These solutions correspond to embeddings which enter the vanishing locus at 
close to 7r/2. Such solutions persist until the angle reaches some minimal value 
(^o)min below which the embeddings simply fall into the black hole. In Figure 48 one 
can see examples of embedding for different values of ^o- Finite temperature adds 
richness to the picture by adding the possibility that these singular solutions may or 
may not be bypassed by the otherwise melting transition induced by temperature. 
This depends entirely on the value of the electric field for a given temperature, and 
for sufficiently high electric field the singular solutions appears unavoidable (at least 
at this level of analysis — as mentioned earlier, there may well be stringy corrections 
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Figure 45: At zero electric field and finite temperature the Minkowski and black hole 
solutions are separated by a "critical" embedding which has conical singularity at the 
event horizon. 




Figure 46: At finite electric field there are two "critical" solutions corresponding to 
embeddings with a conical singularity at the vanishing locus = it*) and at the event 
horizon respectively. The intermediate solutions correspond to black hole embeddings 
which have conical singularity between the event horizon and the vanishing locus. 

to the conical geometry in the interior). The two possible scenarios mentioned above, 
are presented in Figure 47. Let us summarize them again. 

For sufficiently weak electric field or large temperature {E < Ecr) the (0o)min 
embedding is in the vicinity of the "critical" (^o = 7r/2) embedding and is bypassed 
by the phase transition. The thermodynamically stable phases correspond to smooth 
Minkowski and black hole embeddings on gravity side or meson gas and quark gluon 
plasma states in the dual gauge theory. 

For E > Ecr the (^o)min embedding is thermodynamically stable and corresponds 
to some bare quark mass fh* < rhcr- This means that some of the classically stable 
black hole solutions, namely the one corresponding to bare quark mass in the range 
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Figure 47: The physical scenarios described in the text for (top) E < and (right) 
E > Ecr- (These are magnifications, and the small scale wiggles are due to the onset 
of numerical fluctuations at the scale we are observing, and can be ignored.) At the 
top, we have Oq < (^o)min and so the phase transition bypasses the conically singular 
branch of solutions. On the bottom, we have > (^o)mini and so the conically singular 
solutions appear on the resulting phase diagram. See text for further discussion. 

rh* < rh < fhcr, will have a conical singularity. This suggests a richer structure for 
the phase diagram than the one presented in Figure 44 and in particular the existence 
of a special point at E = E^ and ?7icr(-E'cr)- The corresponding phase diagram is 
presented in Figure 49. The new choice of dimensionless variables on the axes is made 
in order to emphasize the existence of the non-smooth {i.e., conically singular) black 
hole embeddings. The solid curve corresponds to the phase curve from fig 44 (but note 
that the axes are different here, and hence the shape) and the dashed curve separates 
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Figure 48: Examples of D7-brane embeddings for > (^o)mirn = (^o)min and < 
(^o)min respectively. The inner semi-circular segments represent the event horizon, 
while the outer semi-circular arcs correspond to the vanishing locus at u = u^,. The 
numerical value of the minimal angle is Oomin ~ 1-345 and the parameter E = 3. 

embeddings with conical singularity from smooth black hole embeddings (that lie below 
it). The area above the solid curve is in the stable meson phase and the area below 
the dashed curve is in the melted/dissociated phase. The area between the curves 
corresponds to the embeddings that have a conical singularity before ending on the 
horizon. The vertical dashed line corresponds to the critical value Ecr ~ 1-26 below 
which the conical solutions are not part of the energetically favorable solutions and so 
do not complicate the story. Note that in the T ^ limit the only energy scale is 
R\fE and the appropriate dimensionless variable is: 



while in the same limit E oo. Therefore the phase curves from Figure 49 qualita- 
tively approach the zero temperature (large electric field) case of the previous section, 
with the critical parameters m* ~ 1.22 and rhcr ~ 1.32. 



m 



= lim 



m 
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Figure 49: The solid curve correspond to the phase curve from Figure 44 and the 
dashed curve separates embeddings with conical singularity from smooth black-hole 
embeddings. The area above the solid curve is in meson gas phase and the area below 
the dashed curve is in quark gluon plasma phase. The area between the solid and 
dashed curves corresponds to embeddings which reach the black hole horizon but have 
a conical singularity along the way. Their physical status remains to be determined. 
(See text.) The vertical dashed line corresponds to the critical value Ecr ~ 1.26. 

5.4 Geometric aspects of the instability 

As we saw in the previous sections, the addition of external electric field results in the 
formation of a new vanishing locus defined by = 6^ + R^E'^. Minkowski embeddings 
which reach that vanishing locus become unstable. To rectify this, the Ai component 
of the U{1) gauge field on the D~brane should be turned on with a non-trivial profile 
along the holographic direction u. The interpretation is that for sufficiently strong 
electric field (equivalently sufficiently small bare quark mass) the binding energy of the 
quarks is completely overcome by the external electric field and mesons get dissociated 
to their constituent quarks, which provides free electric charges. The additional gauge 
field is reflecting this by supporting global electric current along the direction of the 
external electric field. In this way the vanishing locus at u = u^, naturally refiects the 
existence of an insulator/conductor phase transition. 
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However, it would be nice if we could shed more light on the geometrical aspects 
of that transition and in particular on the physical reason for the instability of the 
Minkowski embeddings at the vanishing locus. In this section we will show that in 
an appropriate T-dual description this corresponds to an over spinning of the D7- 
branes, causing them to become superluminal. To illustrate this let us T-dualise the 
gravitational background along the xi direction. The relevant part of the Buscher 
transformation rules [49] is: 



Gii = 


1 




— Gf^i, + 


Bi^Biv — Gi^Giu 
Gn 


Gfxi = 


B^i 
Gn' 




G^i 
Gn' 


" Gn • 



(272) 



Applying the transformations (272) to the background (244) with the constant Kalb- 
Ramond B-field from equation (247) results to the following geometry: 

ds^a' = -''^-p^df+'^idxl + dx'l) + ^^^du^ + R'de^ (273) 
li u It u — b 

+R^ cos^ edn\ + R^ sin^ OddP' + -^dxl + 2^^dtdxi; 

e = gs — ; 5^1/ = . 
u 



In this way we trade the non-zero -B-field for an off diagonal component of the metric. 
Furthermore, the metric of the T-dual background develop an ergosphere at u = -u*, as 
can be seen from equation (273). This suggests that the D6-brane probes (the probes 
loose one dimension because of the T-dualization), which reach the ergosphere, will 
develop speed (along the xi direction) equal to the local speed of light and hence will 
become superluminal beyond the ergosphere, which on the other side, will result in 
an imaginary action. Naively it seems that the D6-branes embeddings should become 
non-static beyond the ergosphere and develop angular velocity along the xi direction. 
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However, there is another possibihty stemming from the fact that the D6-branes are 
extended objects. Indeed, for extended objects the only condition for stability that we 
need to satisfy is that they should not be spacelike, hence if they extend along the xi 
direction they can still describe static configurations. 

Let us make the discussion more explicit by writing down the DBI action for a D6- 
brane which is extended along the t,X2,X3,u directions of the background and wraps 
an described by tp, P, 7. In addition we will consider an ansatz of the probe having 
non-trivial profile along the 6 and xi directions: 



After substituting into the action (246) and using equation (273) we get the lagrangian: 



if one compares equation (275) to equation (251) one can see that they are indeed 
identical (as they should) provided that we define f{u) = xi{u). The rest of the analysis 
is completely analogous to the analysis of the lagrangian (251). This geometric point 
of view can be particularly useful if one is trying to address the universal properties of 
the insulator /conductor phase transition subject of our study. We shall come back to 
that in a later chapter when we will, show that this electrically driven quantum phase 
transition falls into the same universality class as the quantum phase transition driven 
by an R-charge chemical potential. 

Another point that is worth noting is that since the S— field is a pure gauge, 
it is equivalent to a fixing of the U{1) gauge field of the flavor brane, actually this 
is the way an external electric field was introduced in ref. [5, 51]. In the T-dual 
description this equivalence corresponds to the freedom to change coordinates to a 
rotating frame. Indeed, one can easily check that the background (273) is actually 



9 = 6{u); xi = xi{u); 



(274) 




(275) 
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just the background T-dual to the AdSs-BHx S'"'' background (244), but in a rotating 
frame. Note that being able to work consistently in a rotating frame is due to the 
fact that our backgrounds are asymptotically AdSs x and the rotation is along the 
compact directions of the background. In this way a D6-brane probe, corotating with 
the frame, becomes superluminal as it goes deeper into the bulk of the geometry and 
is regular at infinity, which is in accord with the holographic requirements. 

5.5 Concluding remarks 

It is very encouraging that non-trivial non-perturbative phenomena resulting from 
external fields such as those we have seen here (a dissociation phase transition, metal- 
insulator transition and the associated response current) can be so readily extracted 
in this kind of holographic study. We found that since the electric field works together 
with the presence of finite temperature, the resulting phase diagram which accounts 
for the effects of both is rather simple (ignoring the complication of the special conical 
solutions — see below) . 

Clearly, the story is not quite complete in this electric case, since above a certain 
value of the electric field, some of our solutions develop conical singularities in the 
interior. We are not entirely sure about the nature of these solutions. As already 
stated, one possibility is that the solutions are locally (near the singularity) corrected 
by stringy physics, perhaps smoothing the conical points into throats that connect to 
the horizon. We would expect in that case that our phase diagram would largely remain 
intact, since the values of the masses and condensates for each solution are read off at 
infinity, far from where the conical singularity develops. The results for the relative 
free energies of those solution (compared to the other solutions at the same values of 
the mass) would be the same, and so the complete story would be unaffected. A more 
drastic possibility is simply that the physics is considerably modified by instabilities 
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associated with those conical solutions that our current study has not revealed — the 
appropriate part of the meson spectrum is particularly difficult to study numerically 
here, and so the tell-tale signs of tachyonic modes that represent an instability are 
hard to check for in this case. Such a modification could remove quite a significant 
part of the phase diagram for large enough values of the electric field, allowing the 
possibility of a richer structure than we have seen so far. 
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Chapter 6: Phase structure of finite temperature large 
Nc flavored Yang- Mills theory with an R— charge 
chemical potential 

6.1 Introductory remarks 

Not long after the birth of AdS/CFT [41, 66, 86], where the finite temperature phase 
structure of A/" = 4 SU (Nc) (at large Nc) was understood in terms of the thermody- 
namics of Schwarzschild black holes in AdS [86, 87], it was recognized [18] that progress 
could be made in understanding the properties of gauge theories in the presence of a 
global current, akin to a chemical potential for baryon number or isospin, by study- 
ing the physics of charged black holes, such as Riessner-Nordstrom [18, 19], or more 
general charged black holes with non-trivial scalars [22], in AdS. 

The resulting (/xr, T) or {q, T) phase diagrams (where q is conjugate to /XR-the 
analogue of Baryon number) of the A/" = 4 gauge theory were found to be rather rich, 
with e.g., a first order phase transition line, and even a second order critical point in 
the {q,T) plane [18, 19]. This led to the hope that such studies might lead to insights 
into the physics of the QCD phase diagram. The idea is that studies of the effect of 
the global U{1) symmetry-even though it is not exactly baryon number-might lead 
to physics in the same universality class as the more realistic gauge theories, giving 
insight into QCD at finite temperature and density. 
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A more firm footing for tliis idea should be obtained by the study of those dynamics 
in the presence of fundamental flavors of quark, a key feature of QCD that distinguishes 
it from the Af = A gauge theory (which of course only has adjoint matter of very 
specific types and quantities allowed by supersymmetry) . To this end in this chapter 
we consider D7-brane probing of the spinning D3-brane geometry. 



6.2 Spinning D3— branes 
6.2.1 The string background 

To source subgroups of the U{1)^ C SO (6) global symmetries corresponding to ro- 
tations in the transverse M^, we consider D3-branes with angular momentum. The 
supergravity description/metric for spinning D3-branes is [21, 58]: 



ds^ = H~^'^ \ - (l - ^ \ dt^ + dx ■ dx \ + hI'^ 
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In order to explore the decoupling limit {a' 0), which we are interested in, we define: 

r = a'u , rn = a'uu, 4 = , 



such that: 



^ - sinh(/?3)-4V> 



With these definitions, in the decoupling limit, the metric is: 



2 



+A-U^nifj^{Rd(()i-Aidtf + R'^nidiiH , (277) 



, i=l 



where we have defined: 

2 4 2 2 

f=^mn2ns-^, n. = i + ^, ai = '^-^, (278) 

A = 7^i7^2'H3 ^ ;r^, /Xl = siu 0, 1^2 = COS 6 Slli 1p , fls = COS 9 COS tp . 

i=l 

Of immediate interest for us is to compute the temperature of the background and to 
examine its dependence on the different choices of angular momenta qi in that we 
may have. As we are going to show, the case of three equal charges qi = q is the simplest 
scenario containing extremal horizon and hence zero temperature at finite chemical 
potential. This is important because this will enable us to study the first order phase 
transition of the fundamental matter at zero temperature, when the fluctuations of the 
system are purely quantum. Another interesting property of the at three equal charges 
case is the extension of the phase diagram of the dual flavored gauge theory to the 
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case of non-zero R-chargc chemical potential. As we are going to demonstrate later, 
it has the same structure as the one for the adjoint sector of the theory, corresponding 
to the Hawking-Page transition of the gravitational background [18]. Let us proceed 
with the calculation of the temperature. 

6.2.2 Calculation of the temperature 

In order to calculate the temperature of the background we could evaluate the surface 
gravity at the event horizon. This calculation was performed in ref. [78]. Here we are 
going to undertake a more explicit approach. First we will regularize the background 
(277) by going to a rotating frame in which the geometry does not have an ergosphere. 
After that we will re-scale the coordinates to focus on geometry near the event horizon 
of the background, thus going to an appropriate Rindler frame. 

The radius ue of the actual event horizon of the background (277) is determined 
by the largest root of f{u) = 0. In addition it has an ergosphere determined by the 
expression: 

3 
i=l 

Since the background (277) is asymptotically AdSs x S^, we can "remove" the ergo- 
sphere (279), by going to a rotating frame. This is equivalent to gauge shifting A] 
from (278) such that: 

A'; = -Ri,\, + Al (280) 
The parameters /Lt|^ are set by the condition ^t*|uB = hence: 

Prom the asymptotic behavior at infinity {u — ^ oo), it is clear that /Xj^ correspond to 
the angular velocities of the frame along 0^. In the dual gauge theory these correspond 



146 



to having time dependent phases of the adjoint complex scalars or equivalently to an 
R-charge chemical potentials for the corresponding scalars [18]. 
Next we consider the following change of coordinates: 

u = UE + az'^; 6 = 9E-yi; i/^ = ipE - y2 ■ (282) 

Note that the choice of ^Ej V'e is not important for the computation of the temperature 
and the constant a is to be determined. In order to compute the temperature we should 
focus on the following two terms in the metric (277) : 

ds'^/a' = Gttdf + Guudu^ + ... (283) 

and make sure that near the horizon it is of the form: 

ds'^/a'^ = -{27rTfz'^dt'^ + dz'^ + ... , (284) 

which corresponds to a going to a Rindler frame. From equation (277) one can see 
that we have: 

G„„ = AV2/(^^)-i, (285) 

where: 

f{u) = j^,{{u^ + ql){u^ + <i){u^ + ql)-u^ui,) (286) 
= ~ '"e) ~ ""i) - '"2) • 
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Here uf and U2 refer to the other two roots of the polynomial in u^. Note that they do 
not need to be real. After substituting equations (282) and (286) into equation (285) 
and leaving only the term leading in z, yi and y2 we obtain : 

where Ae = A|(„j^_gjj ,0^^). The corresponding value of a is: 

a = . (288) 

Now we focus on the Gu term in (277) but with the new gauge choice for Al given in 
equation (280). The exact expression is: 

3 

Gu = -A"\H{H2H^r^f + A-V2 ^ Hi^4{A!\f (289) 

1=1 

Note that close to the horizon to the leading order in z we have: 

R W + q^ ul + qfj R{ul + qf) 

and hence the second term in equation (289) is of order z^. Thus at leading order in 
z, yi and y2 we have: 

Gudt^ = - ^i)% 2^^2 ^ _^27rTfdt^ , (291) 

u^u^R^ 

where we have used equation (288) for the value of a. Therefore our final expression 
for the temperature of the background is: 
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which is the same as the expression obtained in ref. [78] by calculating the surface 
gravity at the event horizon. 

6.2.3 The extremal ceise and the case of three equal charges qi — q 

It is worth exploring the case when we have extremal horizon of the background (277) 
in some more details. In particular we will consider the simplest case, when have three 
equal charges. The corresponding expression for the temperature is: 

^-2J4('"^^ ''''' 

In order to explore the extremal case let us study the function f{u) determining the 
position of the event horizon of the geometry. It is convenient to introduce the following 
dimensionless parameters: 

u = u/ur; q = q/un. (294) 
The expression for f(u) from equation (286) is then: 

^^""^ = + ^^^^^ + ~ + ■ ^^^^^ 

Therefore we need to study the behavior of the zeros of the polynomial in equation 
(295), as we change the parameter (f. In particular we want to examine for which 
values of the polynomial has a double zero. Therefore we need /{ue) = and 
/'{ue) = 0. The second condition amounts to: 

3ui + 6ulf + 3g'' - 1 = , (296) 

which has the solution: 




(297) 
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Substituting from equation (297) into equation (295) results to the following simple 
expression: 

and therefore we learn that for charges q satisfying: 



^4, (299) 



the background has an extremal horizon at: 



4 = ^4 • (300) 



It is easy to verify that for these values of q and the temperature of the background 
given by equation (293) is zero (as it should be). 

An interesting question is, what is the physical meaning of the parameter uu in 
the dual gauge theory. Recall that in Chapters 3,4 and 5 the parameter ur was 
proportional to the temperature of the dual theory. However, here we have an extremal 
horizon and therefore the dual gauge theory is at zero temperature. The answer to that 
question becomes obvious if we calculate the R-charge chemical potential of the theory 
given in equation (281). Indeed, for the extremal values of q and ue from equations 
(299) and (300) we have: 

= ^6V^g ■ (301) 

and hence we can think of the parameter un as determining the R-charge chemical 
potential of the dual gauge theory. This provides a convenient framework to study 
of the confinement/deconfinement phase transition of the fundamental matter at zero 
temperature. In the next section we will focus on the introduction of fundamental 
matter to the spinning D3-branes background. 
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6.3 Introducing fundamental matter 

In this section we focus on the introduction of fundamental matter to the dual gauge 
theory. The technique that we employ is the same used in the previous chapters, 
namely introducing a stack of Nf D7-branes in the probe limit {Nf <C Nc). Before 
we do that, we should notice that, unlike all other cases that we have considered so 
far, the general spinning D3-brane background has a deformed compact part of the 
geometry. In particular the internal sphere in (277) is deformed by the rotations 
along ^1, (^2 and 03 . Therefore, we cannot directly apply the set up that we have used 
so far, namely the D7-brane extendeds along the non-compact AdS-like part of the 
geometry and wraps an internal inside the part of the geometry. However, there 
is a way out of it if we consider the slightly less general case of having q2 = qs- Indeed, 
in order to restore some of the symmetry of the metric (277), we will consider the case, 
when q2 = qs- This corresponds to having an (parameterized by 02,03) inside 
the deformed S^. 

Now if we introduce D7-branes filling the AdS-like part of the geometry and wrap- 
ping the S^, we will add fundamental matter to the gauge theory. Furthermore, we 
are free to rotate the D7-branes along 0i and the corresponding angular velocity is 
interpreted as a time dependent phase of the bare quark mass. (Recall that in intro- 
ducing D7-branes to the D3-brane system we actually add flavors as chiral superflelds 
into the N = 2 gauge theory) . If that phase is the same as the phase of the complex 
adjoint scalar, /x^i, it is equivalent to an R-charge chemical potential for both the 
adjoint scalar and the chiral field. 

On the gravity side of the description this is equivalent to letting the D7-branes 
have the same angular velocity as the rotating frame of the background. Moving to 
the frame corotating with the D7-brane corresponds to going back to the gauge choice 
for A\ from equation (278). The price that we pay is that we again have an ergosphere 
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in the bulk of the background. It will be also convenient to gauge shift A'^ and A'^ to 
A"f and A"^, correspondingly, so that the ergosphere is at: 

A(7^iHi)~ V - sin^ = 0. (302) 

The shifted forms, A"'^ and A"^, vanish at the locus given by equation (302). In 
analogy to the T-dual description of the previous chapter for the external electric field 
case, the ergosphere embeddings will have to be extended along so that they can 
stay non-space-like beyond the ergosphere. The possible D7-brane embeddings then 
naturally split into two classes: Minkowski embeddings that have a shrinking S'^ above 
the ergosphere and ergosphere embeddings which reach the ergosphere. These classes 
are again separated by a critical embedding which has a conical singularity at the 
ergosphere. The set up is very close to the one employed in the T-dual description of 
the electric case, considered in Chapter 5. In Chapter 7, we shall shed more light on 
this similarity. 

Now let us focus on the properties of D7-branes embedding. There are two different 

cases that we consider: 

(1) . A single charge in the 4>i direction: qi = q and q2 = qs = 0. This case is a useful 

test case since it is the simplest problem to explore and provides valuable insight 
into the next, more complicated case. 

(2) . Three equal charges: qi = q2 = qs = q- This case is of greater interest since it 

is the simples case involving extremal horizon and thus enables us to study the 
gauge theory at zero temperature and finite R-charge chemical potential. 
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6.3.1 Single charge: The general set up 

For the single charge case, we choose qi = q and ^2 = Q'3 = in equation (277), and 
the metric becomes: 

+A-1/2 j-^ gi^2 ^^^^^ _ ^1(^^)2 + ^2 pQg2 ^^j^2| ^ (3Q3^ 



where 



Although there is a change of variables that allows us to rewrite this metric in a more 
conventional form, it will be computationally convenient to work with the metric in 
equation (303). In order to introduce dynamical quarks into our gauge theory, we 
probe the geometry with D7-branes [50]. We choose the coordinates {0,(f)i) to be the 
transverse coordinates to the D7-brane. However, we will define radial coordinates in 
the M^-subspace transverse to the D-3 branes via: 

p = ucos9, L = usin9. (305) 

Although these radial coordinates are not transverse to each other, the non-zero Glp 
component of the metric vanishes at the asymptotic boundary. Therefore, by expand- 
ing: 

L = m+^^+o(^), (306) 



we can still make sense of the coefficient m as being the asymptotic separation of 
the D-3 and D-7 branes. This will allow us to extract the bare quark mass and the 
condensate value. However, it should be noted that unlike the cases studied so far. 
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we cannot directly identify the coefficient c with the quark condensate. The correct 
identification can be made after an appropriate change of variables that restores the 
symmetry of the geometry at the asymptotic boundary. We will explicitly show this 
later. We choose an ansatz 9 = 6{u) and d/^cpi = 0. With these choices, the lagrangian 
of the D7-brane is given by: 



Note that the ergosphere of the background (303) is given by equation (302) which for 
the case in consideration simplifies to: 

u^{u^ + cos^ e)-u^ = Q . (308) 

The event horizon given by the zeros of f{u) from equation (286) is determined by the 
equation: 

4(4 + g^) - 4 = . (309) 

Now it is straightforward to see that, just like in T-dual description of the external 
electric field in Section 5.4, the D7-brane probes, which reach the ergosphere, develop 
speed (along the 4>i direction) equal to the local speed of light and hence become 
superluminal beyond the ergosphere, which results in an imaginary on-shell action and 
hence free energy. Naively it seems that the D7-branes embeddings should become 
non-static beyond the ergosphere and develop angular velocity along the direction. 
However, there is another possibility stemming from the fact that the D7-branes are 
extended objects. Indeed, for extended objects the only condition for stability that we 
need to satisfy is that they should not be sace-like, hence if they extend along the 0i 
direction, they can still describe static configurations. 
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Let us make this discussion explicit and consider the ansatz cpi = 4>i{u) for the 
D7-brane action. The corresponding lagrangian is given by: 



The equation of motion for can be integrated to give us the constraint: 



(1 - sin2 e){l + i?2/^'2) + R^f sin^ e<t)'^ 



Ki . (311) 



We shall discuss the physical meaning of the constant of integration Ki in the dual 
gauge theory in Subsection 6.3.7. Let us use equation (311) to solve for (p[ and obtain 
the on-shell lagrangian: 



^"^X u2,f^ ef^'ff + ^V^^^^'cos^g sing sin V cos . (312) 

If we require that the lagrangian (312) be regular at the ergosphere, we can fix the 
constant of integration: 

Kf = q^u^ cos^ ^0 sin'' Oo . (313) 

The parameters uo and in equation (313) correspond to the coordinates in the {u, 9) 
plane for which the D7-brane enters the ergosphere. Note that uq and 9o are related 
via equation (308). One can sec that for the critical embedding at = ^/2, as well as 
for Minkowski type of embeddings, we have Ki = and the lagrangian (312) reduces 
to the one given in equation (310). 

Now let us address the question about the quark condensate of the theory and how 
is it encoded in the profile of the D7-brane embedding in the (p, L) plane. As we 
commented above we cannot directly relate the coefficient c in the expansion (306) to 
the condensate of the theory. The reason is that because of the deformed it is not 
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so simple to calculate the variation of the action with respect to L. Furthermore, one 
can show that for the coordinates p and L defined in equation (305) the on-shell has 
a logarithmic divergences which require counter term depending on the field L. As we 
saw in Chapter 4, the main advantage of using L as a transverse field is that the counter 
terms are independent on L and we can regularize the action be simply subtracting 
the action for the L = embedding. This is why it is convenient to perform change 
of coordinates that will remove the logarithmic divergences of the action. To this end 
we will restore the symmetry of the as much as possible as we move away from the 
horizon. 

By looking at the distorted part of the background in equation (277), one can 
see that, away from the horizon (i.e. when we see the effect of -uh less and less), it has 
the conformal form: 

3 

dss,'' = Y.^i i^ixj + d^jixj) , (314) 

i=l 

where 

jjLi = sin0 , /X2 = cos sin '0 , = cos cos -0 . (315) 

Following ref. [77], the change of variables that would restore the symmetry is given 
by: 

yi = \/r^ + Qi IJ-i cos 01, 1/2 = \Jr'^ + ql Mi sin 01, ys = ^^r^ + g| /X2 cos 02, 

yA = ^Jr"^ + /^2 sin 02, 2/5 = \/r^ + qI /X3Cos03, ye = \l r"^ + qI At3sin03 . (316) 

In view of this, for the single charge case, one defines the following radial coordinate: 



r = ^ . (317) 

Jq2 + + 2^v?{v? + g2) - u% 
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This allows us to rewrite: 



f-^du^ + R^de^ = ^ + de^] . (318) 



dr'^ 
J.2 _j_ q2 



We can now define new coordinates: 



L = -^7-2 _|_ g2 gjj^ Q ^ p — j-QQgQ ^ (319) 



so that we have: 



dr^ 



R\-^ 2 + '^^ > = ; {dL^ + dp'') . (320) 



Then from the asymptotic behavior of r for large u: 



./:p^^ = u+p + ... , (321) 

^ Lb 



we can argue that L has the expansion: 



~ c + m%- 

L = m+ ^ + . . . . 322 

p2 



Using the standard argument from ref. [60] we can show that: 



- 

{ijj^) ~ -c - my = c' . (323) 



6.3.2 Single charge: Properties of the solution 

Let us now study the properties of the D7-brane embeddings. As we already men- 
tioned, the solutions fall into two different classes, namely Minkowski embeddings 
which have a shrinking and correspond to meson states and ergosphere embeddings 
which reach the ergosphere and eventually fall into the black hole. 
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We will consider first the Minkowski type of embeddings. The relevant lagrangian 
is given in equation (307). The equation of motion for 9{u) derived from (307) is a 
second order non-linear differential equation, thus similar to the cases studied in the 
previous chapter, we will rely on numerical methods to solve it. However, It turns 
out that for values of the parameter g 3> uh one can derive a simple analytic solution 
which proves to be a very good approximation even at q of the order of uh- Indeed, 
for Minkowski embeddings the D7-brane close above the ergosphere and the minimum 
value of the coordinate u is u„iin > uu- This suggests that for large q one can ignore 
the terms in (307) and the lagrangian simplifies to: 

jC(xu^ cos^ e^/v? + q'^cos^eJ ^ ^ ^ + 9''^ . (324) 

V + 

If we change coordinates as in equation (319), the action becomes simply the action 
for a regular AdSs x space, namely: 



We can determine the constant in (326), by requiring that the D7-brane closes at 
u = Umin, namely 9{umin) = 7'"/2. Therefore we have that for large g 3> uh the bare 
quark mass is related to the parameter Umim via: 



£ oc pVl + L'2 . 



(325) 



The solution regular at p = is simply: 



L = v^tiM-g^ sin 9 = const . 



(326) 




(327) 
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Equation (327) is particularly useful because it gives us an estimate of the critical 
mass rricr at which the meson melting phase transition takes place (for large q). More 
precisely, it gives us an estimate of the bare quark mass m* of the critical embedding 
separating the Minkowski and ergosphere classes of solutions (the one having a conical 
singularity at the ergosphere). Indeed, one can show that for the critical embedding 
'^min = "f^H, then we obtain: 



where we have defined the dimensionless parameters rh = m/uH and q = q/uii- It 
is interesting to verify this result numerically. To this end we solve the equation of 
motion of the critical embedding numerically for wide range of the parameter q. From 
the asymptotic of the solution at large u = u/uh can extract the parameter m*. 
The resulting plot is presented in Figure 50. The solid line in Figure 50 corresponds 
to equation (328), one can see the good agreement even at small q. Note that the 
critical embedding is unstable and is by-passed by the first order phase transition of 
the fundamental matter. Nevertheless, the parameter m* is a good approximation of 
the actual critical mass rhcr at which the phase transition happens. 

Now let us study the effect of the R-charge chemical potential on the phase transi- 
tion. To this end we derive the equations of motion for the Minkowski and ergosphere 
classes of embeddings, using equations (307) and (310). Note that the lagrangian in 
equation (310) gives us two equations of motion: for 9(u) and 0i('u). The second one 
can be integrated and substituted into the equation of motion for 9{u) using equa- 
tion (311). Next we define the dimensionless parameters: 



m* = ^/l + q'^ , 



(328) 



u = u/ub_; m = m/«H; c = c/u'^; q 



q/uH . 



(329) 



159 
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Figure 50: Plot of versus q. The solid line in the figure correspond to equation 
(328), one can see the good agreement even at small q. 

After that we solve numerically the equation of motion for 9(u) and using equations 
(305), (306) and (323), we calculate the calculate the parameters rh and c', correspond- 
ing to the bare quark mass and the quark condensate of the dual gauge theory. The 
resulting plot of the equation of state in — c' versus fh coordinates, for different values 
of the parameter q, is presented in Figure 51. 

One can see that the overall structure remains the same as in the q = case. 
However, the critical mass is increasing with q as can be expected from the behavior of 

that we have studied. Furthermore, for sufficiently large q, the ergosphere type of 
embeddings develop a negative condensate — c' which approaches zero from bellow in 
the limit m — > 0. The plot in Figure 52 presents a zoomed in version of the equation 
of state near the phase transition. One can see the typical first order phase transition 
pattern. The corresponding value of the critical mass, calculated using the equal area 
law, is represented by the vertical solid line. It is interesting to study the phase diagram 
of the theory and in particular how does the critical temperature of the theory depend 
on the R-charge chemical potential. 
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Figure 51: A plot of the equation of state for values of the parameter q = 0, 0.5, 1.0. 
One can see that the critical mass is increasing with q. 

6.3.3 Single charge: The phase diagram 

Let us study the phase diagram of the theory. To this end let us obtain expressions 
for the temperature T and the R-charge chemical potential /Ir as functions of the 
parameters uyl and m. Using equations (281), (292) and (309), we obtain the following 
expressions for the temperature: 




(330) 



and for the R-charge chemical potential: 




(331) 



Note that in this set up the temperature T and the R-charge chemical potential are 
not completely independent. Indeed, for the ratio T/fin we get: 

>^. (332) 
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Figure 52: Magnified plot of tlie first order transition pattern of the equation of state 
for q = 0.5. The solid vertical line corresponds to the bare quark mass of the critical 
embedding rhcr ~ 1.126. 

Therefore, in the (/iR,r) plane only states above the line T = ^fJ-R correspond to 
physical states. This suggests that it is more convenient to study the phase diagram 
in the (^ur, T — ^/xr) plane. Furthermore, the change of variables from uh and q to 
/XR and T is invertible in two separate intervals of values of q , namely q"^ <2/ \/3 and 
q"^ > 2/\/3- Of those, only the first interval contains the point q = which corresponds 
to the pure meson melting phase transition considered in Chapter 3. Therefore, we 
shall consider only values of q satisfying q"^ < 2/\/3- Next note that: 



where A is the t'Hooft coupling and niq is the bare quark mass of the dual gauge theory. 
This suggests introducing the following dimensionless physical parameters: 



= ^T^g^a'^N = Xa' 



ruq = m/(27ra') = muH/(27rQ') 



(333) 




(334) 
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The corresponding phase diagram is presented in Figure 53. Note that The curve 
separates meson states from melted meson states. Note that we have used as 
an approximation for the critical mass. We have also assumed that the D7-brane 
embeddings which reach the ergosphere fall into the black hole. In this way we ignore 
effects as the one described at the end of Chapter 5, where we showed that ergosphere 
type of embeddings that are sufficiently close to the critical embedding have a conical 
singularity above the event horizon. We assume that stringy corrections which become 
important near the point of singularity would resolve this singularity and enable the 
D7-brane to connect to the black hole, via a thin tube. 




Figure 53: The area enclosed by the curve corresponds to a confined phase of the 
fundamental matter. The area outside corresponds to melted mesons. 
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6.3.4 Three equal charges: The general set up 

When we set all three charges to be equal [18], {i.e., qi = q2 = q3 = qin equation (277), 

the metric is given by: 

0,2 

R 

3 

+ ^\^B?dn'i + i4{Rd(t)i- A^idtf"^ > (335) 



i=l 



where 



From a computational point of view, it is again more convenient to use 



L{u) = usm.e = m+-^ + 0\^ \ . (337) 



However, just like in the single charge case of last section, the Minkowski solutions 
exhibit a positive slope of q^ /2 for the condensate which is again an artifact of the 
coordinates. To show that, we introduce a new radial coordinate r{u) satisfying: 

j->2 

Hf-^du^ + R^ = -^{dr"^ + r'^dO'^) , (338) 
and the corresponding radial coordinates are: 

L = rsin^, p = r cos 9. (339) 
Using that r{u) from equation (338) has the asymptotic behavior: 



r{u) = u+j^ + 0[^ ) , (340) 
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one can show that: 



L = m+(c + m^)^ + 0[^) . (341) 



Note also that one can write: 



o2 

Hf-^du'^ + R^= ^„ , J dp'^ + dL^) , (342) 
L -\- p 



and following the standard argument presented in ref. [60] one can show that: 



c - my . (343) 



Therefore, to explore the functional dependence of the quark condensate on the bare 
quark mass, we again choose 9 and 0i to be our transverse coordinates to the D7- 
brane. And again we work with the gauge choice A'\ from equation (280) for the 02 
and (^3 directions. With this choice, the world-volume of the D7-brane is given by: 



Co^u^ cos^ eHi l- 4^ sin2 e\ I + ^O'"^ si^^) cos (344) 
y u^R'^j V H 

Again, just like in the single charge case, one can see that the action vanishes at the 
ergosphere given by equation (302). For the case in consideration the equation (302) 
reduces to: 

{u'^ + q^f - u^{u^ + sm^ e) = Q (345) 

For Minkowski types of embeddings we can still use the lagrangian in (344), but for 
ergosphere embeddings we need to generalize the ansatz and let the D7-branes ex- 
tend along the 0i direction, thus enabling them to remain non-space like beyond the 
ergosphere. The resulting lagrangian is given by: 



C oc cos^ 9hJ (1 - sin2 ^)(1 + ^9'^) + ^ sin^ 9<p'^ sin cos V (346) 

y u R J H H 



165 



Now, just like we did for the single charge case, we can integrate the equation of motion 
for derived from (346). This amounts to the equation: 



cos^ OR^fsin^ Ocj)'^ 



^(1 - <^ sin^ e){l + ^6'^) + ^ sin^ e<Pf 



= Ki 



(347) 



For the on-shell action we get: 



>Coc«^cos^0sin^i?A 



H{u^m - 4(1 + 5 sin2 e)) (l + ^9'^ 



sin 2-0 



\ ((u2 + q2)3 - uj^u'^){u'^ + g2) sin2 ^cos6 e-Kf 2 



. (348) 



Now it is straightforward to determine the constant of integration Ki: 



Kf = uy^{u^ + q^) sin* 9o cos^ Oq , 



(349) 



where uq and 6*0 determine the point at which the D7-brane enters the ergosphere. 
Note that Oq and uq satisfy equation (345). 

6.3.5 Three equal charges: Properties of the solution 

The analysis of the three charge case is analogous to the on for the single charge case. 
A crucial difference is that there is an upper bound for the parameter if < 2\/3/3 at 
which the temperature drops to zero. Following the same procedure as in the single 
charge case we can derive the equation of motion for 0{u) and use it to calculate the 
equation of state in the (m, — c')-plane. The resulting plot for various values of the 
parameter q is presented in Figure 54. One can see that the critical mass is again 
increasing as a we increase q. However, since the for the three charge case we have 
extremal horizon at f = 2/(3\/3) and we consider only the case q^ < 2/(3\/3) (no 
naked singularity), the critical mass has an upper limit (unlike the one charge case). 
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Figure 54: A plot of the equation of state for values of the parameter q = 0, 0.3, 0.6. 
One can see that the critical mass is increasing with q. 

6.3.6 Three equal charges: The phase diagram 

Now let us focus on the study of the phase diagram. As we saw earlier in the case 
of three equal charges the temperature of the background is given by equation (293) 
which we replicate below: 



The radius of the event horizon is the biggest root of f{u) given in equation (336). 
In analogy to the one-charge case it is convenient to express that temperature as a 
function of q and uu- 




(350) 




(351) 



where 



3I/32 _^ 2^/^{-9f + v^81^4 _ 12)2/3 



(352) 



62/3 (-9^2 + ^81g4 - 12)1/3 
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Similarly just like in the one charge case we can express the R-charge chemical potential 
as a function of q and un- 



un 6^/^q{-9f + V81g4 _ 12) V3 

itR = — r , . (353) 

i?2 31/32 + 2V3(-9g2 + ^81g4 - 12)2/3 ^ ^ 

One can verify that unlike the one charge case the temperature and R-charge chemical 
potential are completely independent and as un and q vary take all possible values in 
the (r, /XR)-pane. It is once again convenient to introduce the following dimensionless 
parameters: 



^/XR ^ 27T 6^/3^(-99^+ v^Slq^- 12)1/3 
rriq m 31/32 + 2V3(-9g2 + ^8J~43l2)2/3 



mq m \ UE{qrJ 



Now using (the bare quark mass of the critical embedding) as an approximation 
for fhcr, we obtain the phase diagram presented in Figure 55. One can see that there 
is a critical value of the chemical potential above which even at zero temperature there 
are no bound meson states. Note that, just like in the one charge case, we assume 
that all of the embeddings which are type ergosphere eventually fall into the black 
hole and correspond to melted mesons. This is equivalent of assuming that the conical 
singularity of the embeddings close to the critical one is being repaired by a stringy 
corrections. In this way those embeddings would smoothly connect to the black hole. 

6.3.7 The physical meaning of Ki 

Let us focus on the physical interpretation of the constant of integration Ki defined 
in equations (313) and (349), for the single and three equal charges cases respectively. 
The main point is to interpret the physical meaning of the coordinate 01 in the dual 
gauge theory. The correct interpretation is that 0i is contributing to the phase of the 
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Figure 55: The area enclosed by the curve corresponds to a confined phase of the 
fundamental matter. The area outside corresponds to melted mesons. There are no 
bound meson states above the critical value of the parameter \/A^ ~ 6.47. 

bare quark mass. Indeed, in a T-dual description the mass separation between the 
D3- and D7-branes corresponds to the expectation value of the U{1) gauge field on 
the D7-brane [75], and rotation in the transverse plane parameterized by (pi, rotates 
the components of the D7-brane's gauge field which correspond to the components of 
the bare quark mass. Let us make the discussion more explicit. The mass term of the 
theory can be written in general as: 

n = ^{mgOg + m*0;) + no , (355) 

Here TCq is the term of the hamiltonian that doesn't depend on the bare quark mass. 
Usually we can choose coordinates such that rriq is real. In that case {Oq) is the 
corresponding condensate of the theory. Here we assume that: 

mq = \mq\e''^-; {0)g = \{0)q\e-'^^ , (356) 
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then the expectation value of the term (355) can be written as: 



\mg\\{0)g\cosA(f); Ac/) = 0^ - . (357) 

Now we identify the expectation value of the hamiltonian with the regularized action 
via: 

(H) = -Sreg[L,</.i] (358) 

On gravity side, as we learned in Chapter 1, the value of the bare quark mass 
is controlled by the asymptotic separation of the D3- and D7-branes 171 — J-'ooj Via 
\mq\ = m/(27ra'). Therefore we can relate the variations of \mq\ and to the 
variations of and 0i|oo- Then using a standard variational argument one can show 
that: 

|(0),|cosA,^ = -(27ra')^|oo; K||(0),| sin A</. = (359) 



Therefore we see that: 



tan A0 = - ( ) / ( ) = ^ . (360) 



Therefore we sec that the constant of integration Ki is directly related to the phase 
difference between the bare quark mass and the condensate of the dual gauge theory 
A(j). Furthermore, using equation (357) one can write a modified version of equation 
(168) for the absolute value of the condensate of the theory: 



1(^)^1 = T^^^i^ V + ■ (361) 

Note that A^ = for Minkowksi embeddings and A^ ^ for ergosphere type 
of embeddings. Therefore, we conclude that the phase difference A4> is a natural 
order parameter for the first order phase transition subject of our study. Note also 
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the striking similarity to the case of external electric field, considered in the previous 
chapter, when the order parameter of the phase transition is the global electric current 
along the direction of the electric field. In the following chapter, we shall shed more 
light on the discussed similarity. 
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Chapter 7: Universality in the large Nc dynamics of 
flavor 



7.1 Thermal vs. quantum induced phase transitions 

In this chapter we discuss common properties of the phase transitions studied in Chap- 
ters 3 ,5 and 6. In particular we focus on discrete self-similarity exhibited by the theory 
near criticality. To begin with, let us summarize the properties of the set up. In all 
cases considered so far, the background corresponded to the near horizon limit of a 
non-extremal D3-brane background. In addition we introduce D7-branes in the probe 
limit. The D7-brane wraps a, C and extends in the radial direction of the asymp- 
totically AdS5-space. The size of the varies as a function of the radial coordinate. 
The D7-brane embeddings then naturally form two classes: embeddings that reach 
the horizon and hence fall into the black hole, and embeddings for which the wrapped 
5^ shrinks to zero size at some radial position. For these, the D7-brane world-volume 
simply closes smoothly before the horizon. In an Euclidean presentation, the com- 
pact, unbounded parts of the D7-brane have the topology x since the Euclidean 
time has a periodicity set by the inverse temperature of the system. The classes are 
then distinguished by one or the other compact space shrinking away. The authors 
of ref. [12] proposed that the (topology changing) transition of the D7-brane embed- 
dings corresponds to a type of confinement/deconfinement phase transition, now in 
the meson sector of the theory. This system has been extensively studied and it was 
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shown that it is a first order phase transition providing a holographic description of 
the meson melting phase transition of the fundamental matter. 

There is also a unique critical embedding separating those two classes. This solution 
reaches the horizon and has a shrinking S^. It has a conical singularity. Solutions of 
this type will occupy much of our attention in this chapter. Many of these features 
generalize to the general Dp/Dq system. In ref. [67] the Dp/Dq system was considered 
and some universal properties, associated with this critical solution separating the two 
classes of embedding, were uncovered. In particular it was shown that for a certain 
temperature the theory exhibits a discrete self-similar behavior, manifested by a double 
logarithmic spiral in the solution space. This space of solutions is parameterized by 
the bare quark mass and the quark condensate. (Geometrically these correspond, 
respectively, to the asymptotic separation of the D7- and D3- branes and the degree 
of bending of the D7-branes away from the D3-branes.) 

The region of solution space where the self-similar spiral is located is unstable, in 
fact: There is a first order phase transition associated with the physics of the system 
jumping between branches of solutions and bypassing it entirely. Nevertheless, it seems 
that important features of the full physical story can be captured by examining the 
neighborhood of this critical solution. It is remarkable that the critical exponents (or 
better "scaling exponents", so as not to confuse the physics with the nomenclature 
of second order phase transitions) characterizing this logarithmic structure exhibit 
universal properties and depend only on the dimension of the internal S" wrapped 
by the Dg-brane. The precise value of the critical temperature is irrelevant. The 
structure is determined by focusing on the local geometry near the conical singularity 
of the critical Dg-brane embedding, and the exponents are then naturally determined 
by the study of possible embeddings in a Rindler space [32, 67]. 

The studies described above concern a thermally driven phase transition. As the 
temperature passes a certain threshold, thermal fluctuations seek out the new global 



173 



minimum that appears and the system undergoes a transition to a new phase. In 
this chapter we study transitions of the system under the effect of two different types 
of control parameters: an external electric field and an R-charge chemical potential, 
revisiting work done on these systems in Chapter 5 and Chapter 6 of this work. We 
show [30] corresponding scaling exponents are again universal and depend only on the 
dimension of the internal sphere wrapped by the Dq^-brane. We find that the key 
properties of the critical solution can be determined from the local properties of the 
geometry, and we find that this geometry arises naturally by working in a rotating 
frame, arrived at using T-duality. The resulting physics is not controlled by thermal 
dynamics, the local geometry is not Rindler, and so the exponents are different. The 
phase transition is driven by the quantum (as opposed to thermal) fluctuations of 
the system, as can be seen from the fact that they persist at zero temperature. It is 
satisfying that we can cast these different types of transition into the same classifying 
framework. 

7.2 Thermal phase transition 

Let us begin by reviewing the result of refs. [32, 67]. We will be using the notations 
of ref. [67]. Consider the near-horizon black D/j-brane given by: 



where H{u) = (R/u) f{u) = 1 — (uh/u) and i? is a length scale (the AdS radius 
in the p = 3 case). According to the gauge/gravity correspondence, string theory on 
this background is dual to a (p + l)-dimensional gauge theory at finite temperature. 
Now if we introduce Dg-brane probe having d common space-like directions with the 




(362) 
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Dp-brane, wrapping an internal S"^ C 5"^"^ and extended along the holographic coordi- 
nate u, we will introduce fundamental matter to the dual gauge theory that propagates 
along a (d + l)-dimensional defect. 
If we parameterize <S'^~*' by: 

dnl^p = de^ + sin^ ednl + cos^ 0dn'^_p_„ , (363) 

where dfi^ is the metric on a round unit radius m-sphere, the DBI part of the La- 
grangian governing the classical embedding of the probe is given by (We consider only 
systems T-dual to the D3/D7 one, which imposes the constraint p — d + n+ l = 4.): 

£ (X e-*y-|5a/3| = —u"" sin'* Q^YVf^ (364) 

The embeddings split to two classes of different topologies: 'Minkowski' embeddings 
which have a shrinking S*" above the vanishing locus (the horizon) and yield the physics 
of meson states and 'black hole' embeddings that reach the vanishing locus, corre- 
sponding to a melted/deconfined phase of the fundamental matter. These classes are 
separated by a critical embedding with a conical singularity at the vanishing locus, as 
depicted in Figure 56. 

It is convenient to introduce the following coordinates: 

= 1 (^«^ + ^ju-'-v - v^-"^^ , (365) 
L = r cos 6 , and p = r sin 6 . 
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Figure 56: Schematic diagram depicting the Minkowski (left) and black hole (right) 
embedding solutions that are separated by a "critical" embedding (centre) which has 
a conical singularity at the event horizon. 

Then one can show [50, 60] that the asymptotic behavior of the embedding at p — ^ oo 
encodes the bare quark mass iriq = m/2'!ra' and the quark condensate {'ipip) oc — c of 
the dual gauge theory via the expansion: 

^ip) = ^+-^ + --- (366) 

After solving numerically for each embedding of the Dg-brane, the parameters m and 
c can be read off at infinity. Prom the full family of embeddings, a plot of the equation 
of state of the system c(m) can be generated. The resulting plot for the D3/D7 system 
[6] is presented in Figure 57. The two different colors (and line types) correspond to the 
two different classes of embeddings. The equation of state is a multi-valued function, 
and there is a first order phase transition when the free energies of the uppermost and 
lowermost branches match. 

The main subject of our discussion is the spiral structure in the solution space near 
the critical embedding [67, 68]. In the enlarged portion on the right in Figure 57 it is 
located to the lower left, roughly at m = 0.9185, — c = 0.0225. The spiral structure 
that is hidden near this point is a signal of the discrete self-similarity of the theory 
near the critical solution. 
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Figure 57: Plot of the equation of state c(m). The zoomed region shows the location 
of the first order phase transition. There is a spiral structure hidden near the "critical" 
solution in the neighborhood of m = 0.9185, — c = 0.0225. 

In order to understand the origins of the spiral, we zoom into the space-time region 
near the tip of the cone of the critical embedding [32, 67] using the change of variables: 



3-p 

y_ /'^\^ 

R \ r) 



u = U]i + irTz ; 
Here T is the temperature of the background given by 



7-p 



(367) 



AttR \ r) 



(368) 



Leaving only the leading terms in z results in the following metric: 

"2 - -{27rTfz'^dt'^ + dz'^ + dy'^ + y'^d^l + dxj + ... 



ds^ 



(369) 
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The metric (369) corresponds to flat space in Rindler coordinates. The embed- 
dings of the Dg-branes in the background (369) again spht into two different classes: 
Minkowski embeddings characterized by shrinking 5" {y = 0) at some finite zq, and 
black hole embeddings which reach the horizon at 2; = for some finite y = jjo (the 
radius of the induced horizon). The equation of motion is derived from the Dirac- 
Born-Infeld action of the Dg-branes which has the following Lagrangian: 

jC (X + ■ (370) 

The equation of motion derived from this reads: 

zyy" + {yy' -nz){l + y'^) = . (371) 

Solutions of this equation enjoy the scaling property y{z) — j^y{jjLz), in the sense that 

if y{z) is a solution to the equation (371) so is j^y{fiz). Under such a re-scaling the 
initial conditions {zo,yo) for the two classes of embeddings scale as: 

zo zq/h; yo yo/ix] (372) 

This suggests the existence of a critical solution characterized by = yo = 0. One can 
check that y = \fnz is the critical solution. It has a conical singularity at y = 2; = 0. 

To analyze the parameter space of the solutions we can linearize the equation of 
motion (371) near the critical solution by substituting ?/(z) = ^Jnz + ^{z)^ for small 
^(z). The resulting equation of motion is: 

z^^\z) + (n + \){,z^{,z) + i{z)) = , (373) 
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which has a general solution of the form: 



with 



-(Acos(a„lnz) + 5sin(a„lnz)) , 



(374) 



n 



rn = ^ ; o;„ = ^ A/4(n + 1) - . 



Note that q„ are real only for n < 4 which are the cases naturally realized in string 
theory [68]. Now the scaling property of equation (371), combined with the form of the 
solutions (374) suggests the following transformation of the parameters {A, B) under 
the re-scaling of the initial conditions given in equation (372): 



cos (a„ In ^) sin (a„ In ^u) | i A 
— sin (a„ In jj) cos(a„ Itll/j.) I \B 



(375) 



For a fixed choice of the parameters A and B, the parameters {A',B') describe a 
double spiral, whose step and periodicity are set by the real and imaginary parts of 
the critical/scaling exponents r„ and a„. 

Equation (371) has a Z2 symmetry [32] relating the two classes of solutions (black 
hole and Minkowski embeddings) . If the parameters {A, B) describe one class of embed- 
dings, then the parameters {—A, —B) describe the other. In this way the full parameter 
space near the critical solution (given by A = 0, = 0) is a double logarithmic spiral. 

This self-similar structure of the embeddings near the critical solution in our 
Rindler space is transferred by a linear transformation to the structure of the solutions 
in the (m, c) parameter space. If we call (m*, c*) the parameters corresponding to the 
critical embedding from Figure 56, then sufficiently close to the critical embedding we 
can expand: 



m — 


m* 




(a\ 


= M 











+ 0{A^) + 0{B^) + 0{A,B) . 



(376) 
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' m — m* I 






I c-c* 1 







^1 *l\ 


Mm' — m* 




|to — m \ 


\ |c'-c*| ^ 







The constant matrix M cannot be determined analytically and depends on the proper- 
ties of the system. Generically it should be invertible (numerically we have verified that 
it is) and therefore in the vicinity of the parameter space close to the critical embedding 
(m*,c*) there is a discrete self-similar structure determined by the transformation: 

m' — m* \ 1 / cos (a„ In ii) sin (a„ In n) 

d — c* j ^ \— sin {an In /x) cos(a„ In /x) 

Let us define two solutions to be "similar" if: 



(378) 



Then one can see from equation (377) that this is possible only for a discrete set of /xs 
given by: 

gW«". k = l,2,... . (379) 

Note that in general the matrix M in equation (377) will deform the spiral structure 
given by the transformation (375). However, the scaling properties of the theory 
remain the same as they are completely determined by the scaling exponents: r„,a„. 
Furthermore, one can see that the scaling exponents depend only on the dimension 
of the internal sphere S"" wrapped by the Dg-brane and are thus universal, in the 
sense that the detailed value of the critical temperature is irrelevant. It is the spiral 
structure that ultimately seeds the multi-valuedness of the space of solutions, twisting 
the (m, — c) curve back on itself as in Figure 57. Therefore, it is the spiral — and the 
neighborhood of the critical solution from where it emanates — that is responsible 
for the presence of a first order phase transition in the system. Whether there is a 
spiral or not can be read off from the scaling parameters (rn, Q;„), and since [68] for all 
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consistent Dp/Dq systems the condition n < 4 is satisfied the corresponding thermal 
phase transition (meson melting at large Nc) is a first order one. 

7.3 Quantum— induced phase transitions 

In this section we will consider a different class of phase transitions. These are arise in 
the presence of external fields, and can happen even at zero temperature, and so since 
the fluctuations driving the transition are no longer thermal, they might be expected 
to be in a different class. Naively, the broad features of the equation of state — multi- 
valuedness and so forth — have similarities with the thermal case, and so it is natural 
to attempt to trace the extent to which these similarities persist. We will find that 
once we cast these systems in the language of the previous section, the similarities and 
differences will be quite clear. 

We will first concentrate on the case of an external electric field. The flavored 
system, at large enough electric field, has an insulator/ conductor phase transition, as 
studied in Chapter 5. As with the thermal transition of the last section, the mesons 
dissolve into their constituent quarks, but this time it is due to the electric field over- 
coming their binding energy. The transition is of first order. 

As we saw in the previous section the scaling properties of the thermally driven 
phase transition are naturally studied in a Rindler frame with a temperature set by the 
temperature of the background. In Chapter 5 it was shown that in analogy to the ther- 
mally driven phase transition there is a nice geometrical description of the electrically 
driven phase transition, and the structure of the system can be again characterized by 
an unstable critical embedding with a conical singularity at an appropriate vanishing 
locus (analogous to the event horizon). Here, we will generalize this description to the 
case of the Dp/Dq system. 
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Furthermore, after an appropriate T-duality transformation we will show that 
the vanishing locus corresponds to an effective "ergosphere" due to a rotation of the 
coordinate frame along the compact directions of the background. The instability near 
criticality is then naturally interpreted as an instability due to the over-spinning of 
the D(g — 1) brane probes (in the T-dual background) as they reach the ergosphere. 
We then study the structure of the theory near criticality by zooming in on the space- 
time region in the vicinity of the conical singularity. Once again, we will find that the 
structure is entirely controlled by the dimension of the internal sphere, S*", wrapped 
by the T){q — l)-branes (in the T-dual background) — details such as the value of the 
electric field and the temperature of the system, are irrelevant. 

7.3.1 Criticality and scaling in an external electric field 

Let us consider the near-horizon black Dp-brane given by the background in equa- 
tion (362). Following a similar idea [29] for producing a background magnetic field, 

if we turn on a pure gauge i?-field in the {t,Xp) plane [5, 25, 52], in the dual gauge 
theory this will correspond to an external electric field, oriented along the Xp direction: 



B = Edt A dx, 



(380) 



The resulting Lagrangian is: 



C (X e 



-\9al3 + 




1 f-E^H 



«"sin"eVl + /M2^'2 . 



(381) 



This leads to the existence of a vanishing locus at = n* given by: 



,7-p _ 



= u 



J-p 

'H 



(382) 
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at which the action (381) vanishes. Notice that this is distinct from the horizon, and 
even at zero temperature will be present. A study of the local physics near this locus 
will therefore pertain to non-thermal physics. 

The embeddings split into two different classes: Minkowski embeddings which have 
a shrinking 5" above the vanishing locus and correspond to meson states and embed- 
dings reaching the vanishing locus, corresponding to a deconfined phase of the fun- 
damental matter. These classes are separated by a critical embedding with a conical 
singularity at the vanishing locus. Our goal is to explore the self-similar behavior of the 
theory near this critical embedding and calculate the corresponding scaling exponents. 

In order to make the analysis closer to the one performed in refs. [32, 67], for 
the thermal phase transition (described in the last section) , we T-dualize along the Xp 
direction. This is equivalent to a trading of the pure gauge i?-field for a rotating frame 
in the T-dual background. Indeed, the geometry T-dual to equation (362), with the 
5-field given by equation (380), is given by: 

1 ^""^ 1 1 /d ^ \ 

ds^ = H-^-fdt^ + dxj) + 2mEdtdxp + m + u^dnl_p + dxl] , 

i=l ^ •' ^ 

e^ = 9sH^-h /=l-(v)'"'- (383) 

The background given by equation (383) corresponds to the near-horizon limit of a 
stack of Nc D(p— l)-branes smeared along the coordinate otp. Now if we place a probe 
D{q — l)-brane having (d—l) spatial directions shared with the D(p — 1) — branes, 
filling the radial direction u and wrapping an internal inside the S^~^ sphere of the 
background, we will recover the action (381), as we should. 

Note that in these coordinates we have an effective "ergosphere" coinciding with 
the vanishing locus given by equation (382). Now the critical embedding is the one 
touching the ergosphere and having a conical singularity at « = tt*. In the (m, c)-plane 
this embedding corresponds to the center of the spiral structure (m*,c*). 
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Despite the analogy with the analysis of the thermal phase transition, in this case 
there is a crucial difference because of the necessity (from charge conservation) for the 
D{q — l)-brane to extend beyond the ergosphere. Indeed, since the D(g — l)-brane 
is an extended object one can find static solutions that extend beyond the ergosphere 
and are non-superluminal. To this end one should allow the D(g — l)-brane to extend 
along the direction of rotation otp. In the original coordinates (before T-dualization) 
this is equivalent to a non-trivial profile for the Ap component of the gauge field, 
which corresponds to the appearance of a global electric current along the Xp-direction 
[5, 52]. This is the reason why we refer to the corresponding phase transition as an 
insulator /conductor phase transition. After the transition, the quarks are free to flow 
under the influence of the electric field, forming a current. 

Let us describe how this procedure works in the case of a general D(p — 1)/D(g — 1)- 
intersection. Again we will work in the T-dual background (383). Let us consider an 
ansatz for the D(g — l)-brane embedding of the form: 



Now after integrating the equation of motion for Xp and plugging the result in the 
original Lagrangian, we get the following on-shell Lagrangian: 



9 = 9{u) ; 



(384) 



this leads to the action: 




(385) 




(386) 
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It is easy to verify that if we choose the integration constant in equation (386) to 
satisfy: 

= E'^H^uf' sin^" 00 , (387) 

then the action (386) is regular at the ergosphere {u = u*). Note that at the critical 
embedding Oq = 9^ = and the constant in equation (387) is zero. As we learned 
in Chapter 5, this constant is proportional to the global electric current along the Xp 
direction of the original Dp/Dg-brane system. 

We are interested in the scaling properties of the theory, near the critical embed- 
ding solution. Despite the fact that the Lagrangians (381) and (386) describing the 
Minkowski and ergosphere classes of embeddings are different, the fact that at the criti- 
cal embedding they coincide (if^=0) shows that the corresponding equations of motion 
share the same critical solution. Furthermore, as we will see, the critical exponents are 
the same for both types of embedding. 

Let us introduce dimensionless coordinates by the transformation: 

Du^ y /u^\^ fu^\^ ^(u^\^ 

H^Exp^Xp, (388) 

1 

where = (7 — p)^/*/i?* u^, iJ* = H\u=u^ and /* = f\u=u»- To leading order in z 
and y the metric (383) is given by: 

ds^ = -Dzdf + dz^ + dy'^ + y^dnl + H^uldn^_p_^ (389) 

1 

+2dtdxp + -^^^j^dxp + ^2 ■ 

* i=l 
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First consider the case of Minkowski embeddings, characterized by a distance zq 
above the ergosphere at which they close {y = y{zo) = 0). The Lagrangian describing 
the D(g — l)-brane embedding is: 



A oc y"'z 



,/2 



(390) 



The corresponding equation of motion is given by: 




,nl/2 



,/2 



- ny 



n— 



^Z^/'^^/l + ^ = Q . 



(391) 



Equation (391) possesses the scahng symmetry: 



y^y/iJ'\ 



z -> z/ii; 



(392) 



in the sense that if y = y{z) is a solution to equation (391) so is the function j^y{fiz). 
Now under the scaling (392) the boundary condition for the Minkowski embedding 

scales as zq — > zo/fi- This suggests the existence of a limiting critical embedding with 
zo = 0, and indeed: 



is a solution to the equation of motion in equation (391). The corresponding D(g — 1)- 
brane has a conical singularity at y = z = 0. Now before we linearize equation (391) 
and calculate the critical exponents let us consider the case of the ergosphere class of 
solutions characterized by the radius of the ergosphere induced on their world-volume. 
Because of the possibility to extend beyond the ergosphere we should consider the 
analog of the ansatz from equation (384): 



y{z) = V^nz 



(393) 



y = y{z); 



Xp — Xp(z) . 



(394) 
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The corresponding Lagrangian is: 



L oc y'^^jDz{l + + (^Pz + l)x'^ , (395) 

where F = D/E^H^. After integrating the equation of motion for Xp and substituting 
it into the Lagrangian (395) , we obtain the following on-shell Lagrangian: 

(X —j= „ V 1 + ■ (396) 

It is easy to see that the Lagrangian (396) is regular at z = 0, y = yo- The equation of 
motion for y(-z), derived from the Lagrangian (395) and with the substituted solution 
for Xpi^z) is: 



(397) 

It is easy to check that equation (393) is a solution to equation (397). Furthermore, 
for z ^ 1/F one can see that equation (397) has the scaling symmetry (392) (note 
that equation (392) suggests yo yo/ lA- Linearizing equations (391) and (397) near 
the critical solution (393) by substituting: 



y{z) = ^/2nz + i{z) (398) 

results in the same equation: 

z'eiz) + (n + l/2){zi'{z) + i{z)) = . (399) 
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The general solution of equation (399) is given by: 



^{z) = — (^cos(Q!„ln2;) + 5sin(a„lnz)) , 



(400) 



where the scaling exponents are given by: 



2n- 1 




(401) 



Note that the scaling exponents again, while quite different from those of the thermal 
case (see equation (374)) depend only on the dimension of the internal S"' wrapped by 
the Dg-brane and are thus universal for all Dp/Dq systems. Furthermore, the discrete 
self-similarity holds for n < 5. By similar reasoning to the thermal case [68], since for 
all consistent systems realized in string theory we have that n < 4, for such systems we 
may expect that the electrically driven confinement/deconfinement phase transition is 
first order and has the described discrete self similar behavior near the solution that 
seeds the multi-valuedness of the equation of state. 

The rest of the analysis is completely analogous to the thermal case considered in 
the previous section. Therefore we come to the conclusion that close to the critical 
embedding (specified by and c*) the theory has the following scaling property: 



with r„ and a„ given by equation (401). 

It is interesting to compare the analytic results some numerical studies. Let us 
consider the D3/D7 system. From equation (392) on can see that the variation of the 
scaling parameter n in equation (402) can be traded for the variation of the boundary 
conditions of the probe, namely zq for Minkowski and yo for ergosphere embeddings. 



m' — m' 




(402) 



v 



c — c 
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On the other hand, close to the critical embedding, the change of coordinates in equa- 
tion (388) suggests that: 



^0 oc yo and uq — oc zq , 



(403) 



where uq and a-re the boundary conditions for the embeddings in the original (not 
zoomed in) background. Note that the parameter uq is related to the constituent quark 
mass Mc [68] (in the absence of an electric field) via Mc = (uq — un)/{2Tra'). 
Close to the critical embedding we have that: 



for some fixed boundary conditions ucin and 0o,in- Now equation (402) suggests that for 
Minkowski embeddings the plot of (m — m*)/('Uo — 'f^*)^""'"^ versus q;„ ln('Uo — ■"*) should 
be an harmonic function of ln(uo — u=f) with a period 2tt. Similarly for ergosphcrc 
embeddings the plot of (m — m*)/9Q"-~^^ versus a^ln^o should be a harmonic function 
of ccnln^o with a period 27r. Note that the physical meaning of 9o can be related to 
the value of the global electric current (see equation (387) and the comment below). 

As can be seen in Figure 58, for both types of embeddings the numerical results are 
in accord with equation (402) and the analytic results improve deeper into the spiral 
(large negative values on the horizontal axis). Our numerical results confirm that the 
critical exponents are indeed ra = 5/4 and 03 = \^/4, as the general analytic results 
yield. 

7.3.2 Criticality and scaling with R— chetrge chemical potential 

Now we study the case when the external parameter is an R-chargc chemical potential 
in the dual gauge theory. Introducing flavor to this background was the subject of our 



u*)/{uo - II*) and fi = 2;o,m/2o , 



(404) 
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Figure 58: The solid line is a fit with trigonometric functions of period 27r. The plots 
confirm that the scaling exponents of the theory are rs = 5/4 and 03 = \/3T/4. 

study in Chapter 6. Let us for convenience replicate some of the important equations 
here. 

The relevant geometry is given by: 

3 



where 



/=^WiW2W3-^, W. = l + 4' Al = t^^2 (406) 
3 2 

A = 'Hi'H2W3 with //i=sin0, /i2 = cos^sin'i/', /i3 = cos0cos'^/'• 



^=l 
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Here the parameter ttn would be the radius of the event horizon if the angular momen- 
tum of the geometry was set to zero {qi = 0) . The radius ue of the actual event horizon 
is determined by the largest root of f{u) = 0. The temperature of the background 
that we calculated in the previous chapter is given by: 

^ - ^ - - - - ^) - ("-"'X"-"^) • (^°^) 

where u\ and U2 are the other two roots of f(u) = 0. 

The background (405) has an ergosphere determined by the expression: 

3 

A{HiH2H3)-^f - HifijiA^^ = . (408) 
1=1 

Since the background (405) is asymptotically AdSs x S*^, we can "remove" the ergo- 
sphere (408), by going to a rotating frame. This is equivalent to gauge shifting A] 
from (407) such that A[^ = —R/J-li + A]. The parameters are set by the condition 
^riwE = hence: 

Prom the behavior at infinity (u — > oo), it is clear that fx^^ correspond to the angular 
velocities of the frame along . In the dual gauge theory these correspond to having 
time dependent phases of the adjoint complex scalars or equivalently to R-charge 
chemical potentials for the corresponding scalars [18]. 

In order to restore some of the symmetry of the metric (405), we will consider 
the case when q2 = Qs- This corresponds to having an (parameterized by V', '/'2, 1^3) 
inside the deformed S^. Now if we introduce D7-branes filling the AdS-like part of the 
geometry and wrapping the S^, we will add fundamental matter to the gauge theory. 
Furthermore, we are free to rotate the D7-branes along 0i and the corresponding 
angular velocity is interpreted as a time dependent phase of the bare quark mass. 
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(Recall that in introducing D7-branes to the D3-brane system we actually add flavors 
as chiral superfields into the J\f = 2 gauge theory). If that phase is the same as the 
phase of the complex adjoint scalar, fj,^t, it is equivalent to an R-charge chemical 
potential for both the adjoint scalar and the chiral field. 

On the gravity side of the description this is equivalent to letting the D7-branes 
have the same angular velocity /x^ as the rotating frame of the background. Moving 
to the frame corotating with the D7-brane corresponds to moving back to the gauge 
choice for Aj from equation (407). The price that we pay is that we again have an 
ergosphere in the bulk of the background. It will be convenient to gauge shift A''^ and 
A'^ to A"f and A"^, correspondingly, so that the ergosphere is at: 

A(WiWi)- V - Til sin^ e{Alf = 0. (410) 

The shifted forms, A"^ and A"'^, vanish at the locus given by equation (410). 

The possible D7-brane embeddings then naturally split into two classes: Minkowski 
embeddings that have a shrinking above the ergosphere and ergosphere embeddings 
which reach the ergosphere. These classes are again separated by a critical embedding 
which has a conical singularity at the ergosphere. In analogy to the T-dual description 
of the previous subsection for the external electric field case, the ergosphere embeddings 
will have to be extended along (f>i (as we saw in Chapter 6) so that they can stay 
non-space-like beyond the ergosphere. However, in this chapter we are interested 
in the scaling properties of the theory for parameters (m, c) in the vicinity of the 
critical parameters (m*,c*), corresponding to the critical embedding. As we saw in 
the previous section, modifying the ergosphere class of embeddings so as to be regular 
at the ergosphere does not alter the properties of the theory near the critical solution. 
In particular the scaling exponents characterizing the discrete self-similar behavior of 
the theory remain the same. So henceforth we will focus on the study of the Minkowski 
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type of embeddings. The analysis is completely analogous to the one performed in the 
previous subsection. 

In order to focus on the space-time region close to the conical singularity of the 
critical embedding, we consider the change of coordinates: 



where 

^^erg = ul- qI (412) 

is the radial coordinate u of tip of the critical embedding or equivalently the 6 = ir/i 
point of the ergosphere. It can be shown that for the values of q2 for which the geometry 
is not over spun (and so has an horizon) the corresponding value of Uerg is real. 
After leaving only the leading terms in z and y, we get: 

V? 

ds'^/a' = -Dizdt^ + dz^ + dy^ + y'^d^ll - 2qidtd(pi + ^dx'^ + Ridcpj , (413) 

R 

2 2 

j<~>3 j/2 , -2 ; J i2 , 2 / 7 /2 7-> 4(/i«H j-,2 ^erg + ^1 r)2 

aO = dtp + sm ipd(p2 + cos V"03 ! -^i = ^3 ; Ri = — ^2 ^ ■ 



u 



The metric in equation (413) is of the same type as that in equation (389), namely 
flat space with some compact directions in a rotating frame. Therefore the analysis is 
completely analogous to the one for the electric case and hence the scaling exponents 
are again given by equation (401) with n = 3, because the D7-branes are wrapping an 
internal S^: 

ra = 5/4 ; = Vsi/A . (414) 

We can again verify this numerically. It is convenient to do this for the single charge 
case, namely qi 7^ 0, 52 = 93 = 0. The plot analogous to Figure 58 for the electric case, 
is presented in Figure 59. The plot represents the variation of the bare quark mass 



193 



parameter m as a function of the initial boundary condition uq — uu, for Minkowski- 
type embeddings. The parameter m=K corresponds again to the bare quark mass for 
the state corresponding to the critical embedding. The good agreement with the result 
for the critical exponents in equation (401) is clear, and the accuracy of the analytic 
description improves as we go deeper into the spiral (to the left). 

I m — m» I 
(«o-«h)"* 




Figure 59: Plot of the relation between the bare quark mass parameter m and the 
distance above the ergosphere (uq — uh)- The plot is for q = 0.5 in units in which 
uh = 1. The solid line is a fit with trigonometric functions of a period 27r. The plot 
confirms that the scaling exponents of the theory are = 5/4 and 03 = \/3T/4. 

An important observation is that our result does not depend on the values of the 
R-charges, nor the temperature. In fact, this physics persists at zero temperature, 
such as at extremality with all three charges equal qi = q2 = Q3 = Q, or more generally. 
(Extremality is when ue = ui or U2 for which T = 0. See equation (407).) The 
fact that we have the same structure at zero temperature (extremal horizon) further 
confirms that the key properties of the corresponding phase transition is indeed driven 
by the quantum (rather than thermal) fluctuations of the system. 

7.4 Criticality and scaling: Some generalizations 

In this section we generalize the procedure for the study of the critical behavior em- 
ployed in all three different systems of phase transition (thermal, or in the presence 
of electric field or R-charge chemical potential). This may lay the groundwork for 
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other types of phase transitions that may arise in future studies, seeded by spirals 
with different universal behavior. 

Note that in all cases there is some vanishing locus. The different classes of Dq- 
brane's embeddings are being classified with respect to whether they fall into that 
vanishing locus, or wrap an internal 5" sphere that is contracting to zero size above 
the vanishing locus signaling the end of the Dq'-brane. 

In all cases there is a critical embedding that separates the two classes of embed- 
dings. The critical embedding reaches the vanishing locus and has a conical singularity 
there at some finite radius (u* = uh or Uerg for the thermal and R-charge cases). 

The main point of the analysis is that after zooming into the space-time region 
near the conical singularity we obtain the metric: 



where D is a non-essential constant. The Dirac-Born-Infeld Lagrangian of the brane 
is then: 



Note that to extract the key behavior (that we are studying) of this critical embedding 
(and its neighborhood) there is no need to modify the embeddings which reach the 
vanishing locus (as we did for the ergosphere class of embeddings). The critical solution 
and the linearized equation of motion is the same for both classes. Therefore it is 
sufficient to consider the Minkowski type of embeddings and analyze the Lagrangian 
(416). The resulting equation of motion is: 



ds^ = - 



(415) 



(416) 




(417) 
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It is easy to check that equation (417) has the scahng property (392) and the hmiting 
critical solution is given by: 

y*(^) = ■ (418) 



Now after the substitution: 



In 

y{z) = J-z + az), (419) 



we obtain the following linearized equation: 

z^fiz) + (n + k/2){ze{z) + a^)) = . (420) 
The general solution of equation (420) can be written as: 

e(^) = ^(^cos(a(f))lnz) + Ssin(aWlnz) , (421) 



where 



r(f) = {n + k/2 - l)/2; = ^ V4(n + ^/2) - (n + ^/2 - 1)2; (422) 



are the scaling exponents characterizing the self-similar behavior of the theory. Both 
being real, they control the shape of the spiral which emanates from the critical solu- 
tion. The oscillatory behavior is present for n < 3 + 2^/2 — k/2. For these values of 
n the theory exhibits a discrete self-similarity and the equation of state c = c(m) is 
a multi-valued function suggesting that the corresponding phase transition is a first 
order one. 

While there is the possibility of complex scaling exponents and hence possibly 
second order phase transitions (if the multi-valuedness goes away when the spiral 
does) , this is not realized in the examples that we know from string theory. 
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Note that we have k = 2 for a thermal induced phase transition and k = 1 for 
the quantum induced phase transitions that we studied (external electric field and 
R-charge chemical potential), arising from the two most natural types of a vanishing 
locus that one may have: an horizon, and an ergosphere. Perhaps other systems will 
yield different values of k. 

7.5 Concluding remarks 

We have succeeded in casting two important types of phase transition (in large Nc 
gauge theory with fundamental flavors) into the same classifying framework as the 
meson-melting phase transition. These quantum fluctuation induced transitions (so- 
called since they persist at zero temperature) , resulting in the liberation of quarks from 
being bound into mesons as a result of the application of an external electric field, or a 
chemical potential for R-charge, turn out to have the same underlying structure. It is 
distinct from that found for thermal fluctuation induced transitions. The structures are 
controlled by the local geometry of the spacetime seen by a critical D-brane embedding 
(it is the borderline case between two physically distinct classes of embedding), and 
while it is Rindler for the thermal case with an horizon at the origin, it is (after a 
T-duality in order to geometrize the discussion as much as possible) a rotating space 
with a simple "ergosphere" type locus. The technique of characterizing the physics in 
terms of this underlying classifying space [32, 67] is rather pleasing in its utility, and 
we extended our analysis to the natural generalization of this space, extracting the 
scaling exponents that might pertain to physics from future studies. 

Of course, there is much interest in how much we can learn about finite Nc physics 
(for applications to systems such as QCD) by studying universal features of large Nc- 
Unfortunately, it is almost certain that much of this is far from robust against 1/Nc 
corrections. The spiral structure is rather delicate, and the stringy corrections arising 
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in going away from the large Nc limit would generically severely modify the classifying 
spacetimes we've been studying, erasing the spiral and its self-similarity. The absence 
of the spiral is necessary for there to be (at best) a second-order transition at finite 
Nc, since it results in multi-valuedness of the solution space, requiring the system to 
perform a first order jump. 

However, It is tempting to speculate that the nature by which the spiral is de- 
stroyed by 1/Nc corrections might (especially since the setting is so geometrical) be 
charactcrizable in a way that allows universal properties of the second (or higher) order 
phase transitions to be deduced from the properties of the spiral at large Nc- 
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Chapter 8: Conclusion 



In this work we employed holographic techniques to study dynamics of flavored strongly 
coupled large Nc non-abelian gauge theories in the quenched approximation. In par- 
ticular we focused on the gauge theories dual to the Dp/D{p + 4) systems and most 
of our study was dedicated to the D3/D7 system. One of the main properties that we 
studied, was the phase structure of the theory under the influence of various control- 
ling parameters such as: external magnetic field, external electric field and R-charge 
chemical potential. We were also interested in describing non-perturbative phenom- 
ena such as chiral symmetry breaking. Another aspect of our study was the meson 
spectroscopy of the dual gauge theory. 

In our study of the zero temperature flavored large N(. N = 2 Yang-Mills theory 
coupled to an external magnetic fleld we showed that, as expected, the supersymmetry 
is completely broken by the external fleld. Furthermore, in order to lower its energy, 
the theory develops a quark condensate that scales with the appropriate power of the 
energy scale which is set by the external magnetic fleld. Our study thus provided 
a holographic description of the phenomena of magnetic catalysis of chiral symmetry 
breaking. An intriguing property of the theory is the spiral structure at the origin of the 
quark condensate versus bare quark mass plane. In particular it is interesting how this 
pattern governs the spontaneous breaking of the chiral symmetry. The detailed anal- 
ysis of this structure, performed in Chapter 2, revealed a discrete self-similar behavior 
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of the dual gauge theory near criticahty. Wc calculated the corresponding critical ex- 
ponents for the bare quark mass, the quark condensate, and the meson spectrum. Our 
study of the meson spectrum confirmed the expectations based on thermodynamic con- 
siderations that the lowest positive branch of the spiral corresponds to a stable phase of 
the theory and that the inner branches are real instabilities that are characterized by a 
tachyonic ground state and cannot be reached by a supercooling. The lowest negative 
m branch of the spiral is tachyon free and thus could be metastable. Further study of 
the meson spectrum revealed a Zeeman splitting of the energy levels and characteristic 
Gell-Mann-Oakes-Renner relation [34] between the mass of the pion associated to the 
softly broken chiral symmetry and the bare quark mass. A very interesting direction 
of future study would be the generalization to the case of the Dp/D(p + 4) system. In 
particular, it would be interesting to verify the universality of the pattern in which the 
spontaneous chiral symmetry breaking takes place. 

Another important issue that we address in this work is the thermal properties of 
the gauge theory dual to the D3/D7 brane intersection. Following the pioneer paper 
ref. [12], in Chapter 3 we introduced flavor D7-brane to the near horizon limit of the 
black 3-brane geometry. Our numerical study revealed a first order phase transition 
pattern in the equation of state in the quark condensate versus bare quark mass plane. 
This phase transition has been extensively studied in the literature and is now believed 
to be describing the melting of mesons into the A/" = 4 super Yang-Mills plasma. The 
study of this phase transition and in particular, the generalization of the corresponding 
phase diagram to include the influence of additional controlling parameters is one of 
the main subjects of our study. 

Encouraged by the rich physical picture accessible through the holographic tools 
employed in Chapter 2 and Chapter 3, in Chapter 4, we considered the influence of 
an external magnetic field on the meson melting phase transition. Inversely, one could 
say that we studied the effect that the temperature has on the magnetic catalysis of 
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chiral symmetry breaking. Our study of the corresponding phase diagram supported 
the expectations from earher field theory results that the temperature has a restoration 
effect and increase the value of the magnetic field at which spontaneous chiral symmetry 
breaking happens. The meson melting phase transition exists only below a critical 
value of the applied field. This is the critical value above which spontaneous chiral 
symmetry breaking is triggered (in the case of zero mass). Above this value, regardless 
of the quark mass (or for fixed quark mass, regardless of the temperature), the system 
remains in a phase with a discrete spectrum of stable masses. Evidently, for these 
values of the field, it is magnetically favorable for the quarks and anti-quarks to bind 
together, reducing the degrees of freedom of the system, as can be seen from our 
computation of the entropy. Meanwhile, the magnetization is greater in this un- 
melted phase. Our study of the meson spectrum revealed the existence of a chiral 
symmetry breaking metastable phase existing bellow a critical ratio of the magnetic 
field and the temperature square. This metastable phase eventually wins at sufficiently 
strong magnetic field and the melted phase seizes to exist. We also verified that the 
meson spectrum of the pion, associated to the spontaneous chiral symmetry breaking, 
exhibit Zeeman splitting of the energy levels at weak magnetic field. Furthermore, for 
sufficiently strong magnetic field, when the chiral symmetry is spontaneously broken, 
a characteristic Gell-Mann-Oakes-Renner relation [34] was observed. 

A natural extension of the phase diagram of the meson melting phase transition 
is the addition of an external electric field. In Chapter 5 of this work we addressed 
this problem. A key technical feature that makes the electric case different from the 
magnetic case is the existence of a vanishing locus (above the horizon) at which the D7- 
brane action vanishes. In order to be able to extend the D7-branes beyond that locus 
one needs to excite additional U{1) gauge field on the D7-brane. This additional gauge 
field sources holographically a global electric current along the direction of the external 
electric field [52]. The D7-branes then naturally fall into two different classes, namely 
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D7~branes which reach the vanishing locus and eventually fall into the black hole, and 
D7-branes which close above that locus. The scenario is quite similar to the case of 
finite temperature studied in Chapter 3. The transition from one type of embeddings 
to another type of embeddings corresponds to transition from melted to un-melted 
phase of the fundamental matter and correspondingly from zero global current to 
non-zero one. This is why we refer to this transition as to an insulator /conductor 
phase transition. The bound quarks dissociate under the separating electric force. It 
is somewhat intriguing that this process takes place also at zero temperature. We 
addressed this issue in more details in Chapter 7. Another interesting feature (again 
technical) is that even after we have regularized the D7-branes at the vanishing locus, 
by introducing additional U{1) gauge field on the D7-brane, the embeddings, entering 
the vanishing locus above some critical polar angle, develop a conical singularity before 
falling into the black hole. Unlike the conical singularity of the critical embedding 
separating the two classes of embeddings, this conical singularity does not lead to a 
discrete self-similar behavior. We proposed that perhaps this singularity is fixed by 
stringy corrections. Of immediate interest is the calculation of the meson spectrum 
for this class of embeddings as it could shed light on the question of their stability. 
We leave this task for the future. It would also be interesting to study the meson 
spectrum at zero temperature in the dissociated phase and check for the existence of 
quasi-normal modes, such a study could verify the stability of that phase. 

Another controlling parameter that one can consider is an R-charge chemical po- 
tential. In order to address this problem, in Chapter 6, we introduced flavor D7-brane 
to the geometry of spinning D3-branes. We provided a detailed calculation of the 
Hawking temperature of the background and elaborated on the extremal case for three 
equal charges when we have vanishing temperature. We also calculated the R-charge 
chemical potential of the theory at zero temperature. Then we focused on the proper- 
ties of the D7-brane embeddings in two special cases, namely one R-charge and three 
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equal R-charges. Wc obtained the equations of state in a numerical form for both cases 
and demonstrated that the critical mass is increasing with the R-charge of the geome- 
try. We also obtained the phase diagrams of the meson melting phase transition in the 
temperature versus R-charge chemical potential plane. While studying the properties 
of the probes, in both cases, we found a subtlety closely related to the one that we 
encountered in Chapter 5, for the study of the electric field case. This is the existence 
of a vanishing locus, this time due to the ergosphere of the background. Letting the 
probes extend along the direction of rotation of the geometry initially transverse to the 
D7-brane, solved the problem of extending the solutions beyond the ergosphere. Since 
this direction is isometry, there is a corresponding constant of integration (analogous 
to the global electric current from Chapter 5). It was shown that the physical inter- 
pretation of this constant of integration is related to the phase difference between the 
bare quark mass and the quark condensate of the dual gauge theory. An interesting 
direction of future study is the analysis of the meson spectrum of theory, especially in 
the deconfined phase. 

Finally we would like to comment on the universal properties of the gauge theories 
holographically dual to the Dp/Dg systems T-dual to the D3/D7 intersection. One 
such property was described in ref . [67] , where the authors studied the meson melting 
phase transition for the general class of Dp/Dq systems and revealed a remarkable 
self-similar structure of the theory. The structure is controlled by the local geome- 
try of the space-time seen by a critical D-brane embedding (it is the borderline case 
between two physically distinct classes of embedding) which is a Rindler space. This 
critical behavior is characterized by scaling exponents that can be calculated analyt- 
ically. They are the same for all Dp/Dq systems and depend only on the dimension 
of the internal sphere wrapped by the probe branes. In Chapter 7, we reviewed this 
work and showed that similar structure exists for the quantum fluctuation induced 
insulator/conductor phase transition studied in Chapter 5. The theory again exhibits 
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discrete self-similar behavior, characterized by scaling exponents which are distinct 
from the one for the thermal phase transition but are again universal for all Dp/Dq 
systems. After a T-duality, in order to geometrize the discussion, we showed that the 
local geometry controlling the structure is a rotating space with a simple "crgosphere" 
type locus. We also considered the case of R-charge chemical potential driven phase 
transition, studied in Chapter 6, for the case of the D3/D7 system. We found the same 
structure, as for the insulator /conductor phase transition, characterized by the same 
scaling exponents. From technical point of view this is not surprising since in both 
cases the vanishing loci are of "ergosphere" type. It would be interesting to extend 
this study to the case of spinning Dp-brane geometries as this would provide a further 
check of the proposed universal behavior. Let us also clarify that these studies are in 
the context of the large Nc limit and most likely the observed self-similar behavior is 
an artifact. Furthermore, at finite Nc the phase transition is expected to be of second 
order. There is, of course, the possibility that some universal properties of the Dp/Dq 
systems will still be present at finite Nc- 
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